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Der Resonanzbegriff bei linearen 
gewohnlichen Differentialgleichungen zweiter Ordnung 
RUDOLF IGLISCH | 
Vorgelegt von L. COLLATZ » 
§ 1. Problemstellung, Hauptergebnisse 
In der Differentialgleichung 
(1) y+ alt) y’ + b(t) y = f(A) 
mogen alle drei Koeffizientenfunktionen die gleiche — nicht notwendig kleinste — ; 
Periode P besitzen: 
(2) Gaeta ai BOE) bt), fh FP) = f(z). 
Die zugehérige homogene Gleichung ist 
(3) gp + alt) p + b(t) p =0. ve 
Man setze 7 
t \ ‘ 
(4) p (i) =exp (Jali) ai), also p'—ap=O. ‘ 
Dann verifiziert man unmittelbar durch Rechnung, daB die Funktion - 
(5) 2() = 2) 9 
der sog. zu (3) adjungierten Gleichung 
(6) 2’ —a(t)z' + (0) —a'(t))z=0 
geniigt. Obwohl in der vorliegenden Note unmittelbar weiter nichts davon ge- 
braucht wird, soll in § 2 der Zusammenhang der mit P periodischen Lésungen y 
von (3) und (6) naher untersucht werden aus folgenden beiden Griinden: Erstens : 


werden die dabei zu benutzenden SchluBweisen auch spaterhin wiederkehren ; 
zweitens werden die Ergebnisse in einer anschlieBenden Arbeit bendtigt, die dem 
Resonanzbegriff der zu (1) entsprechenden nichtlinearen Differentialgleichung 
gewidmet ist. Der kundige Leser wird die Resultate dieses Paragraphen unschwer 
in einen umfassenderen Rahmen einordnen k6énnen. 


Definition 1. Bei der inhomogenen Differentialgleichung (1) legt der Reso- 
nanzfall vor, wenn die adjungierte homogene Gleichung (6) eine mit P periodische 
Lésung 
(7) z(t) =P) op) =h(E+P) eC +P) =20 + P) 
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besitzt, fiir die . 


(8) SeOQeOtOQdt=C+0 


0 

ausfallt. 
Dann ist nach (2) auch a(é) z(t) mit P periodisch, also infolge (4) auch 

p'(t) p(t), von selbst natiirlich 2’(t) = (p (4) (t))’. Daher gilt wegen p y= (fp @)'— p' 


(9) pt+ P)g'(t+ P) =P) ol). 


In §3 wird dieser Resonanzfall untersucht. Das Ergebnis laBt sich zusammen- 
fassen in 

Satz 1. Im Resonanzfall nimmt jede Lésung y(t) von (1) — unabhdngig von den 
Anfangswerten — mit unbeschrankt wachsendem t beliebig groke Werte an, u.zw. 
wachsen diese Ausschlage im wesentlichen mindestens linear mut t tiber alle Grenzen 
(vgl. dazu Satz 5). 


Definition 2. Bei (1) liegt der Hawptfall vor, wenn (6) keine mit P periodische 
Lésung z(¢) besitzt. 
In § 4 erhalt man 


Satz 2. Im Haupftfall gibt es fiir alle Werte von t beschrankt bleibende Losungen 
y(t) von (1), z.B. die eindeutig existrerende mit P perrodische Losung. 


Definition 3. Bei (1) liegt der Ausnahmefall vor, wenn zwar (6) mit P perio- 
dische Lésungen (5) besitzt, wenn jedoch fiir alle diese periodischen Lésungen 


P 
(10) SPO we ial 0 
gilt. 

In §5 wird gezeigt: 

Satz 3. Auch im Ausnahmefall gibt es fiir alle Werte von t beschrankt bleibende 
Lésungen, z.B. die mit P periodischen, die aber jetzt eine ein- oder zweiparametrige 
Schar bilden. 

Im Falle konstanter Koeffizienten auf der linken Seite von (1) laBt sich diese 
Gleichung — unter Heranziehung der Fourier-Entwicklung von /(‘) — explizit 
losen. Die hier aufgestellten Satze finden ihre Bestatigung. 


Es sei noch bemerkt, daB bei den nachfolgend durchgefiihrten Abschatzungen 
kein Wert darauf gelegt wurde, jeweils die besten Schranken aufzusuchen. 


§ 2. Die periodischen Lésungen der homogenen 
und der adjungierten Gleichung 
Bewiesen werden soll folgender 


Satz 4. Besitzt die homogene Gleichung (3) eine mit P periodische Loésung q, (t) 


und ist p(t) davon unabhangige Lésung von (3), so gibt es ein bis auf einen konstanten 
Faktor eindeutig bestimmtes 


(14) P(t) =x Q(t) +A q2(2) 


Resonanz bei linearen Differentialgleichungen 184 


derart, dag 
(12) z(t-+P) =p(t+P) e¢+P) =P) vo) =z) 


mit P periodische Lésung der adjungierten Gleichung (6) ist, falls p(t) nicht die 
Periode P besitzt. Ist y(t) auch mit P periodisch, so sind alle Funktionen z(t) 
aus (5) mit P periodisch. 


Zum Beweise multiplizieren wir die aus (3) mit g(t) bzw. q,(¢) gebildete 
Lagrangesche Identitat mit # (¢): 


(YP — PH1) +40 (Y' GL — HHI) =9, 
integrieren sie zwischen ¢ und t+ P, wenden auf das erste Integral partielle 
Integration an und beachten die zweite Beziehung (4). Dann erscheint 


(13) (Pre — Pept? =0 

oder mit Benutzung der Periodizitat von q, (é): 

(14) gi) [A (E+P) o'E+P)—P o'O]—H [A E+P) pE+P)—2 9 (OJ=0. 
Nun sei /, eine — sicher existierende — Stelle, an der 

(15) p (4) =0 

erfiillt ist, demnach @ (¢,) ==0. Aus (14) wird 

(16) P(A +P) vo (4 +P) =f (4) va). 


Dieser Beziehung geniigen samtliche Losungen q(t) von (3). Jetzt soll durch 
passende Wahl von x und / in (11) zusatzlich die Gleichung 


(17) Pla +P) pa +P) =p (4) vA) 

erfiillt werden. Das liefert fiir x und A nach Eintragen von . die Beziehung 

(18) [PG +P) —P(4)]aG)~+[PG+P) o2(4+P) — (4) 2(4)] 4 =0. 
Fall I. Ist 


Pla +P) — b(4) +0, 


so ]4Bt sich aus (18) zu vorgegebenem A eindeutig x berechnen, es gibt also eine 
bis auf einen Faktor eindeutig bestimmte Funktion q(t), die (17) befriedigt. — 
G2 (t) ist tibrigens in diesem Fall nicht mit P periodisch; denn bei gegenteiliger 
Annahme wiirde aus (14) mit g,(¢) statt m(é) fiir /=2, wegen (15) folgen 


P1 (4) G2 (4) PG SR IZ) = p(4)] =0, 
wahrend alle drei Faktoren von Null verschieden sind. 
Fall II. Ist 


p(t, +P) —p(4) =0, damit aus (4) allgemein p(¢+P) =/(é), 
so vereinfacht sich (18) unter Beachtung von f(t) =O zu 
A[2(t + P) — %2(4)] =0. 


ee 
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a) Ist g(t, +P) + M2(4), mithin p,(/) nicht mit P periodisch, so erhalt man 
14=0, x beliebig; g(t) ist wieder bis auf einen Faktor eindeutig bestimmt. 

b) Ist 
(19) Galt, +P) = G2(h), 


so sind x und 4 beide willkiirlich. Aus (14) mit y(t) statt y(t) folgt in diesem Falle 
(20) p> G+ eS 2 (4). 


Daher besitzt nach (1), (2), (19) und (20) auch @,(t) die Periode P; damit ist 
nach (11) jede Lésung p(¢) von (3) mit P periodisch. Es gilt also folgendes 
Zwischenresultat. Ist nur y(t) mit P periodisch, so ist p(t) durch (17) bis auf 
einen Faktor eindeutig bestimmt; sind ,(¢) und q,(t) beide mit P periodisch, 
so gentigt jedes p(t) aus (11) der Beziehung (17). 
Jedes der Gleichung (17) geniigende q(¢) liefert nun nach (5) ein 2(¢) mit 


(21) z(, +P) =z). 
Aus 
2) =PO PO +PH oO =20 ['O +40 9] 
kommt unter Beachtung von (16) 
(22) 2h + P) = 2'(h). 
Nach (2), (6), (21) und (22) ist mithin dieses z(¢) mit P periodisch. — Jetzt folgt 
Satz 4 unmittelbar aus dem Zwischenresultat. 


Zusdtze. 1. Die mit P periodischen Lésungen von (3) und (6) haben die gleiche 
Zahl von Parametern. — Geht man namlich von einer mit P periodischen Lésung 
2 (¢) von (6) aus, so liefert die gleiche Uberlegung eine mit P periodische Funktion 


%(t)- exp (—f a(t) dt) = p(t) p() exp (— J a(?) dt) = y(t). 


II. Stets ist a(t) =f'(d): p(t) mit £(t)>0, also (3) selbstadjungiert, wenn man 
mit p(t) hinaufmultipliziert. Aber: Dann und nur dann hat diese selbstadjungierte 
Gleichung lauter mit P periodische Koeffizienten, wenn #(t) mit P periodisch ist. 


III. (3) und (6) kénnen héchstens dann eine zweiparametrige mit P periodische 
Lésungsschar besitzen, wenn p(t) mit P periodisch ist. 


Zwei Beispiele mégen den Sachverhalt erldutern. 
Erstes Beispiel. Die Gleichung 


[(4 + bcos t) p’]’ + esint gy’ + [a+ (2b —c) cost] ~ =0 
hat eine Losung 
G,(t) = sing. 


Damit der Faktor von "’ positiv ist, sorge man fiir a> |b|. Man errechnet in 
tblicher Weise 


@o(t) = C- sind - 7a (a+ boost)? ay 


sin? ¢t 
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Das Integral laBt sich bei ganzzahligem Zahlerexponenten elementar berechnen. 
_ Man erhalt z.B. mit willkiirlichem « 


Putye == Ot) == a cost, 
Miri c=420 P(t) =alacost+ d], 
fiir c=3b q(t) =a[(a?+ 0?) cost + 2a + O71]. 


In den ersten beiden Fallen ist g(t) wie q(t) mit 27 periodisch, im dritten nicht, 
falls b=-O ist. In jedem Fall besitzt dagegen f(t) die Periode 22; denn man 


errechnet 
—b) sint 1 
pt) =exp| a at] = Ae tee 
b t cme 
i nights (a+b cost)? 


Zweites Beispiel. Jetzt ein Beispiel fiir den Fall einer nichtperiodischen Funk- 
tion p(t). Die Gleichung 


{[> | Ce 1p cos #| ob + 
t [27 er 1 — BA) sin t+ 75") e+ Ly — 1 + 6%) cos] p =O, 


wo die erste eckige Klammer dauernd positiv gemacht werden kann, hat die mit 
2z periodische Lésung 


y—a 


Q, (t)-= Bcost+ sint 4 B 


Man rechnet 
p(t) =exp All sintdt exp ak dt lk 
B eee te cost P at(2 Be x y—1—r) cos c| 


Der erste Faktor ist mit 27 periodisch, der zweite nicht. Beide Integrale lassen 
sich elementar auswerten. Man erhalt 


y (y—a) 


pl) =|at (ty I B?) cos ar PR=VI=w x 
Toate (yam : 
exp es —— arc tg t ( pie nd) tg — 
|e a 1 p) a+ (73 32 = y—1—B?) 


Man beachte, daB man an den Unendlichkeitsstellen des Arguments des arc tg 
immer weiterwachsend auf den folgenden Ast des arc tg tibergehen muB. — 
(2 (t) kann nicht mit 2a periodisch sein, man braucht es nicht erst zu rechnen. 


§ 3. Der Resonanzfall 
Die adjungierte homogene Gleichung (6) habe eine mit P periodische Lésung 
(5), fiir die (8) gilt. Im Gegensatz zur Aussage des Satzes 1 sei die Existenz 
einer Lésung y(¢) von (1) angenommen, die fiir alle Werte von ¢ beschrankt 
bleibt, fiir die also etwa 


(23) 


y(t\| SE 
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gilt. Wendet man auf die aus (1) und (3) gebildete Lagrangesche Identitat die 
gleichen Schliisse an wie in § 2, so erhalt man analog zu (13), wenn diesmal iiber 
ein Intervall der Lange n- P integriert wird, die Beziehung 

t+nP 


(24) LP! 9 pes J P(t) p(t) f(t) ae 

oder unter Beachtung der Periodizitat von p(t) y(t) und der Gleichungen (8) 
und (9) 

(25) pelt) [ye+n"P) —y'O]-—POeO [ye+nP) —y@] =n. 


Wegen der Periodizitat sind die Faktoren vor den eckigen Klammern beschrankt, 
also etwa 


(26) lpQeOQ SP, [POH WS”. 


Aus (25) kommt dann folgende er fiir den ersten Summanden: 
(27) D\y'(¢+nP) —y'(t)| 2n|C| —2ME, 


wobei die Annahme (23) beachtet wurde. Integriert man (27) zwischen ¢ und 
t+41 und benutzt noch einmal (23), so erscheint 


4DE = n| Cl 2M E, 
(28) 2E[2D +M]=n|C|, 


See lCipes 
= 22 D5 


Da n beliebig groB angenommen werden kann, ist dies ein Widerspruch gegen 
(23). Man kann das Ergebnis auch so aussprechen: 


Satz 5. Im Resonanzfall gibt es fiir jede Lésung y(t) von (1) in gedem Intervall 


(29) tStsi+nP+1 
mindestens eine Stelle t*, fiir die 
*)| > Sad 


ausfallt. — Damit ist Satz 1 in vollem Umfang bewiesen. 


§ 4. Der Hauptfall 

Die homogene adjungierte Gleichung (6) habe keine mit P periodische Lésung 
z(t). Entsprechend Satz 2 ist zu beweisen, daB (1) eine mit P periodische Parti- 
kularlosung hat. Dies ist ein einfaches und allgemein bekanntes Resultat; es sei 
aber doch erlaubt, hier einen Beweis im Rahmen der bisher benutzten SchluB- 
weisen anzufiihren. 

Zunachst sei einmal die Existenz einer mit P periodischen Lésung y(#) von 
(1) angenommen. Aus (24) folgt dann mit 7 =1 


y') [bE +P) +P) — 2) oO] —y@ [PE+P) +P) -— 20 0] 
t+P 


= JPW) PFW at 


(31) 
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fir y=1,2. Die Koeffizientendeterminante D(t) dieses Gleichungssystems fiir 
y(t) und —y(é) ist unter Beriicksichtigung der zweiten Beziechung (4) und 
bY, =(PY,)'—?' 9, = (Py,)'—4p¢,, ferner der ersten Gleichung (2) sowie (5) 
_la€+P)—a® a¢+P)—40 

ag(t +P) —2o(t) a(t +P) — a(t] 
Spater wird gezeigt werden, daB 
(33) D(t)=+0 _ fir alle ¢ 
gilt. 


(32) D(t) 


Jetzt soll die mit P periodische Lésung y(t) von (1) konstruiert werden. Man 
rechne dazu aus (31) fiir einen speziellen Wert ¢* die Anfangswerte y(¢*) und 
y’(t*) aus, die eine etwa vorhandene mit P periodische Lésung y(t) haben miiBte. 
Mit diesen Anfangswerten bestimme man y(t) aus (1). Es ist zu zeigen: Dieses 
y(t) ist von selbst mit P periodisch, d.h. es gilt von selbst 


(34) SET chee e) UE oe and «ES Spey (t*) 
Aus (1) und (2) folgt dann ja die behauptete Periodizitat. 


Zum Beweise subtrahiere man (31) fiir /=?¢* von (24) mit m=1 und gleich- 
falls ¢==7*: 
b (+P) 9, (8 +P) [y+ P) — ye) — 


35 
pt — p(i* + P) o,(* + P) [y (* + P) — y@*)] =0. 
Setzt man 
6 Avy —|@ (t) —21(¢) 
on ea 


so ist (35) fiir »=1, 2 ein lineares homogenes Gleichungssystem fiir die beiden 
eckigen Klammern mit der Koeffizientendeterminante A(f*-+ P), wie man mit 
Hilfe der gleichen Umformung feststellt, die (32) als Koeffizientendeterminante 
von (31) lieferte. Bekanntlich ist aber A(t) +0 fiir alle ¢. (Der Beweis der Tat- 
sache lauft iibrigens ganz analog dem gleich durchzufiihrenden Beweis von (33).) 
Daher ist aus dem Verschwinden der eckigen Klammern in (35) unmittelbar (34) 
abzulesen; die eindeutige Existenz einer mit P periodischen Lésung y(¢) von (1) 
ist somit sichergestellt. — Natiirlich kann man nun dieses periodische y(¢) aus 
dem Gleichungssystem (31) fiir alle ¢ unmittelbar berechnen, ohne auf die Dif- 
ferentialgleichung (1) zuriickzugehen. 

Um den Beweis fiir (33) nachzutragen, differenziere man (32) und driicke die 
zweiten Ableitungen 2, nach (6) aus. Man erhalt so 


D'(t) =a(é) D(t), also D(t) =c-exp(fa(t) dt) =cp(t). 
Es ist zu zeigen, daB c+-0 sein muB. Bei gegenteiliger Annahme c=O ware 
D(t) =0, d.h. 
[a(t+P) —2()]- [a+ P) — 40] =[2¢+P) — 2] [a+ P) — 40). 
Daraus erhalt man 


Ah) a(n we 2 (t+ P) —2, (#) 
37) 7 UES fo rrr er 


wobei jede Gleichung so lange gilt, als der Nenner von Null verschieden ist. 


= CONSt, 
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Fall I. Es gibt einen Wert é), fiir den beide Nenner in (37) gleichzeitig ver- 
schwinden. Dann gilt bei jeder Wahl von « und # in 


(38) a(t) =a%(t) + B22(?) 
(39) 2(to + P) = 2(t)- 
Man kann nun « und f so wahlen, da auch 

(40) 2" (to EP) = 2’ (t) 


wird. Aus (39), (40), (6) und (2) folgt dann aber, da dieses z(t) die Periode P 
besitzt, was gegen die Voraussetzung des Hauptfalles bei Gleichung (1) verstéBt. 


Fall II. Wenn keine gemeinsame Nullstelle der beiden Nenner in (37) existiert, 
so gibt es ein Wertepaar A, uw derart, daB fiir alle Werte von ¢ 


Ag (t+ P) + w2e(t +P) =Ay(t) + 1 20(2) 


ist. Das bedeutet aber wieder die Existenz einer mit P periodischen Lésung z(t) 
von (6) und fiihrt zum gleichen Widerspruch wie in Fall I. 


SchluBbemerkung. Da man die allgemeine Lésung von (1) erhalt, indem man 
zu der mit P periodischen Lésung y(#) eine beliebige Lésung «q@, (¢) +f @.2(t) der 
homogenen Gleichung (3) hinzufiigt, bekommt man weitere fiir alle Werte von ¢ 
endlich bleibende Lésungen, falls ,(¢) oder q,.(¢) fiir alle ¢ beschrankt bleibt. 


§5. Der Ausnahmefall 


Hat die adjungierte homogene Gleichung (6) eine einparametrige Schar mit 
P periodischer Lésungen z (¢), fiir die jedoch (10) gilt, so ist die erste Gleichung (31) 
von selbst erfiillt; zur Bestimmung der Anfangswerte y(t*) und y’(¢*) wie im 
vorigen Paragraphen bleibt also nur die zweite Gleichung (31) iibrig. Wie vorhin 
folgt, daB jede durch die so berechneten Anfangsbedingungen bestimmte Lésung 
von (1) von selbst die Periode P hat. Es gibt also jetzt eine einparametrige Schar 
mit P periodischer Lésungen y(#) von (1). 

Hat (6) eine zweiparametrige Schar mit P periodischer Lésungen z(t), fiir 
die (10) gilt, so sind die beiden Gleichungen (31) von selbst erfiillt. In diesem 
Fall sind also alle Lésungen y(t) von (1) mit P periodisch. 
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Der Resonanz fall bei nichtlinearen 
gewohnlichen Differentialgleichungen zweiter Ordnung 


RUDOLF IGLISCH 


Vorgelegt von L. COLLATZ 


§ 1. Problemstellung, Hauptergebnisse 
In der Differentialgleichung 


(1) y+ gly, 9’, t) =f(e) 
sei fiir alle ¢ 
(2) gly, v,¢+P) =el(y,y’,t) und f¢+P) =f). 


Der Einfachheit halber werde etwa vorausgesetzt, daB g(y, v’, t) stetige partielle 
Ableitungen bis zur zweiten Ordnung besitzt und f(t) stetig ist. Es sei eine mit 
P periodische Lésung y(¢) von (1) bekannt. Betrachtet wird dann bei geringer 
Abanderung der rechten Seite die Gleichung 


(3) zg" + g(z,2',t) =f(t) + BF (ie) 
mit 
(4) e+ P) = Fi) 


bei geniigend kleinem |8|; es sollen Aussagen tiber die zu y(¢) benachbarten 
Lésungen z(¢) von (3) gemacht werden. Mit 


(5) z(t) =y(Z) + u(¢) 


liefert Subtraktion von (1) von (3) 


(6) BPS U.N bye 895i, 1) = PE) 
oder nach Anwendung der Taylor-Formel 
(7) Wl Ley ly. yt) W + By (9, V's 0) @ = BF) — FG yy 0? + By yu! + By yo’; 


dabei soll durch das Zeichen ~ angedeutet werden, da Zwischenwerte als Argu- 
mente in die betreffenden Funktionen einzusetzen sind. Betrachtet werden nur 
Lésungen w(t) mit 

(8) |u(t)| Se und |wu'(t)|<e 
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bei geniigend kleinem ¢. Die Rolle der homogenen linearen Differentialglei- 
chung™® spielt jetzt 


(9) P+ By (VY, EP + by (VV ft) P = 9- 
Die Funktion # (t) wird analog I, (4) 

t \ 
(10) p(t) =exp(S ey (yO, 9,4) at); 


die zu (9) adjungierte homogene Gleichung (entsprechend I, (6)) braucht hier 
nicht extra angeschrieben zu werden. 

Satz 1. (Resonanzfall.) Die zu (9) adjungierte homogene Gleichung habe eine 
mit P periodische Losung 


(11) z(t) =P) - pd), 
fiir die 
P 
(12) PRAT AC RACE sete! 
ausfallt; dann nimmt jede Lésung u(t) von (6) — unabhangig von den Anfangs- 
werten — mit wachsendem t Werte von der Grifpenordnung VIB | an. — Evine rechte 


Seite der GréBenordnung |B| in (6) erzeugt also Auslenkungen der gréferen GroBen- 
ordnung \/{B|. (Vgl. hierzu Satz 4.) — Der Beweis wird in § 2 erbracht. 

Satz 2. (Hauptfall.) Hat die homogene Gleichung (9) keine nut P perrodische 
Lésung y(t), so hat die inhomogene Gleichung (6) Lésungen u(t), die fiir alle Werte 
von t von der Grépenordnung | B| bleiben, z.B. die in diesem Fall eindeutig bestimmte 
mit P periodische Lisung. (Vgl. die Satze 5 und 6.) — Den Beweis hierfiir findet 
man in § 3. 

Im sog. Ausnahmefall, wenn also die zu (9) adjungierte homogene Gleichung 
zwar mit P periodische Lésungen (11) besitzt, fiir die jedoch 


[FOO pat =0 


ist, lassen sich ohne weitere Voraussetzungen iiber g(y, y’, t) keine allgemeinen 
Aussagen machen. 
§ 2. Der Resonanzfall 

Sei also eine (12) gentigende mit P periodische Lésung (11) der zu (9) adjun- 
gierten Gleichung vorhanden. Verfahrt man mit der aus (7) und (9) gebildeten 
Lagrangeschen Identitat wie in I §§ 2 oder 3, so gelangt man statt zu I, (25) zu 
13) b(t) p(t) [w(t-+nP) —w')| —p) 9’) [ue +nP) —u(t)] =nBC— 

t+nP 
a af {yy u? “= Zee" uu’ - ey u'?} p (t) Y (t) dt. 


Es soll zunachst gezeigt werden, daB die beiden Ungleichungen 
(14) |u| SN|Blt und |u'()|< NIB 


* Die Lektiire dieser Arbeit setzt die Bekanntschaft mit der Note des Verf. ,,Der 
Resonanzbegriff bei linearen gewohnlichen Differentialgleichungen‘‘ [Diese Z. 3, 
179—186 (1959)] voraus. Zitate von Gleichungs- oder Paragraphennummern aus 
dieser friiheren Note werden durch eine vorgesetzte I gekennzeichnet. 
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mit einer spater zu bestimmenden Konstanten N nicht fiir alle Werte von ¢ gelten 
kénnen. Mit den Bezeichnungen I, (26) gelangt man wie von I, (25) zu I, (27) 
jetzt von (13) zu 


(15) D|u'(t +n P) — u'(t)|=|B| nC —2MN|B|}— S|B| N2n, 
wo “aus der Beschrankung des Integrals in (13) kommt. Daraus erhalt man unter 
nochmaliger Benutzung von (14) 
2N|B|}(D + M) 2 n|B| (|C| — £N?). 
Wahlt man jetzt N so, daB 


(16) |G] — ZN? = B> 0 
wird, so bekommt man einen Widerspruch fiir 
p> 2ND+M) 


B\pi 
Das Ergebnis 1aBt sich auch so aussprechen: 
Satz 3. Mindestens in einem Intervall der Lange 
seal Elie el 

B |p| 
mug entweder einmal | u(t)| >N|B|? oder | u'(t)| >N|B 


Um gemaB Satz 1 eine analoge Aussage fiir |w(é)| allein machen zu k6énnen, 
benotigt man den 


3 geworden sein mit N aus (16). 


Hilfssatz 1. Bet den beiden Differentialgleichungen 


(17) wu’ + a(t) wu’ = git) 
und 

(18) v’ +av’ =G(t) 
sev 


= 00), 2) — Gl). =H, 


(19) V' (to) = U' (lo) = U0 > 0, 
v (fp) =u (to) = to: 
Man setze 
(20) m = Min (e*, 1). 
Dann gilt im Intervall 
(21) REMESTES TS 
(22) u'(t) = v'(t) = ugm — = (e*— 1), 


falls die rechte Seite dieser Abschatzung posit ist. 
Zum Beweise setze man 
(23) w(t) =v(/) — u(4) 
und subtrahiere (17) von (18): 
w +aw’ = — (x —a(t)) uw’ + (GQ) — gd). 
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AuBerdem ist 
w (ts) =O = w' (tp), w'' (ty) > 0. 
Fiir ¢< fp ist 


t 
w'(t) = eb) f exh) [— (a, — af)) w + (GH — g(t))|dt<0, 
to 


solange u’>0 ist fiir <fy, also erst recht, solange dort v’>0 ist, wie ein Blick 
auf (23) zeigt. Aus (18) rechnet man 


t 
w(t) = elit) f G (t) eX —*) dt + up eet), 
ty 


also nach der zweiten Beziehung (19) und Berechnung des Integrals 


v(t) = ug et 9) — an 


= [et Ho) is 1] 5 


Betrachtet man nur das Intervall (21) und benutzt (20), so erhalt man aus der 
letzten Abschatzung die Behauptung (22) des Hilfssatzes 1. 


Die Anwendung dieses Hilfssatzes geschieht so: 
Wir wollen die Annahme 
(24) |u(é)|<T|p|* fiir allet 


mit einer spater zu bestimmenden Konstanten T7<N zum Widerspruch fihren. 
Dazu schreiben wir (7) so um: 


(25) a” Se {gy-(¥, NG t) ie By yr U a 5 Byy u'} Ui = BF) ¥G Sy (y, a, t) Caos a byy uw. 


Als Stelle 4) des Hilfssatzes wahlen wir eine Stelle, wo w’|(¢))|>N|6|? ist, was 
nach Satz 3 geht, wenn (24) angenommen wird. Wir wollen zunachst 


(26) u' (to) > N|B|? 
voraussetzen. 


In (25) laBt sich wegen (8) die geschweifte Klammer nach unten durch eine 
Schranke « abschatzen, der Betrag der rechten Seite nach oben durch eine 
Schranke der Form 


H =(p|[F+UT%) +|B[ART, 
wenn man (24) beachtet. Verwendet man zu (25) als Vergleichsgleichung 
Oak Ud 


und beachtet (26), so liefert der Hilfssatz gemaB (22) 


(27) u'(t) 2 mN [Blt — |p| =+ 2 ce" 4) —RT|pjbS—* 


fiir das Intervall (21), falls die rechte Seite positiv ist. Wir wollen sie sogar 
groBer machen als 27|6|!. Dazu braucht man nur fiir 


HN ae Tee ae O (V/\81) See 


od 
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zu sorgen, d.h. fiir 


ee west e~ —4 
mN —O(fp|) >|2+R=—*] 7 
oder 


(28) Tez 


Integration von (27) zwischen ty—1 und ¢, liefert dann die zu (24) widerspruchs- 
volle Beziehung 


(29) | u (to) — “(fp — 1)| > 2T ||P. 


Diese Ungleichung gilt sogar ohne die Absolutstriche. Die gleiche Uberlegung 
fiihrt aber auch zu (29), falls statt (26) 


u(t) << — N 


ae 
angenommen wird; in diesem Falle wiirde aber linker Hand in (29) w (¢;—1) — u (ty) 
zu stehen haben, so daB (29) mit den Absolutstrichen in beiden Fallen richtig ist. 
Aus (24) und Satz 3 liest man ab den 
Satz 4. Mindestens in einem Intervall der Lange 
2N(D+M) 


ay pray 
B\B 


mup einmal |u(t)|>T|B 
Damit ist aber Satz 1 bewiesen. 


3 geworden sein mit T aus (28). 


§ 3. Der Hauptfall 

In diesem Paragraphen wird vorausgesetzt, daB die homogene Gleichung (9) 
keine mit P periodische Lésung g(t) besitzt. Zunachst soll bewiesen werden der 

Satz 5. Zu geniigend kleinem |B| besitat (1) genau eine kleine Losung u(t) mit 
der Periode P. 

Zum Beweise benétigt man einen bekannten Hauptsatz aus der Theorie der 
impliziten Funktionen, der hier formuliert sei als 

Hilfssatz 2. p(A, B, B) und p(A, B, B) mogen in einer gewissen Umgebung des 
Wertetripels (0, 0, 0) stetige erste Ableitungen besitzen; auBerdem sea 


(30) y(0,0,0) =O und w(0,0,0) =0, 


ferner 
(31) ¢a(0,0,0) @p(0, 0, 0) 


+0 
y4(0,0,0) Pp (0, 0, 0) 


Dann gibt es zu jedem geniigend kleinen B eindeutig ein kleines Losungssystem A, B 
der Gleichungen 
(32) g(A, B,B)=0 und (4, B,f) =0. 

Die Anwendung des Hilfssatzes geschieht so: Man charakterisiere jede Lésung 
w(t) von (6) durch ihre Anfangswerte 


(33) u(o.) =A, w(0) =B, 
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also u=u(t; A, B, B). Dieses uw ist dann und nur dann mit P periodisch, wenn 


gleichzeitig 


=u(P; —u(0; A, BB) = 
(34) p (A, B, B) (PsA, B, B) u(0; : , B) O, 


y(A, B, B) =w'(P; A, B, B) — u'(0; A, B, B) =0 


gilt. Aus (6) sieht man unmittelbar, daB w(¢; 0, 0, 0) =0 gilt, so daB (30) erfullt 
ist. Die Ableitungsvoraussetzungen des Hilfssatzes sind nach bekannten Satzen 
iiber die stetige Abhangigkeit der Lésungen von (6) von dem Parameter 6 und 
den Anfangswerten A, B gleichfalls erfiillt. Differentiation von (6) nach A oder B 
liefert, daB 

gf =ua(t; 4, Bp) und gf) =up(t; A, B,8) 
Loésungen von 


pt’ + ey(y+u,y +u',t) oe +e,(y +4, ¥ +w,t) p* =0 


sind. Also erweisen sich 
gy, (t) = ug (t;0,0,0) und p(t) = ug(é; 0,0, 0) 


als Lésungen von (9), u.zw. sind sie linear unabhangig, da 


gi(0)=1, (0) =0 


y 


P2(0)=0, 92(0) =1 
ist. Die Determinate (31) wird D*(0), wenn 


Milé+P)—gill) G¢2(@+P) — galt) 


Dee, ; i 
gilé+P)— gilt) g2(t +P) — ga(2) 


gesetzt wird. DaB D*(t) +0 ist ftir alle Werte von ¢, wird genauso bewiesen wie 
I, (33), nur hat man statt I, (6) jetzt unsere Gleichung (9) zu verwenden. Damit 
sind alle Voraussetzungen des Hilfssatzes 2 erfiillt. Es existiert mithin zu jedem 
gentigend kleinen 6 genau eine mit P periodische kleine Lésung w(t) von (6). 
SchheBlich brauchen wir noch 
Satz 6. Die mit P periodische Lisung u(t) von (6) geniigt fiir alle t einer Un- 
gleichung der Form 


(Gell | «(¢)| < const |B]. 


Beweis. Genau wie I, (34) erhalt man jetzt fiir dieses mit P periodische u(t) 
die Beziehungen (vy = 1, 2) 


u'(t) [b(t +P), ¢ + P) — PM) y )] — ul) [PE +P) gE +P) —2O G0] 
t+P t+P 


as B EM) p (2) Py (2) dt -s a Woy u> ote 285° uU u' Se by u’?} p(t) Py (t) dt. 


Die Koeffizientendeterminante ist wieder D(?) aus I, (32), also fiir alle Werte 
von t von Null verschieden. Die homogene Gleichung (9) und die dazu adjun- 
gierte (vgl. I, (6)) haben ja nach I, § 2, Zusatz I genau die gleiche Anzahl mit P 
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periodischer Lésungen, bei uns also keine. Mit Hilfe der Cramerschen Auflésungs- 
formel errechnet sich nun 


{+P t+P 
rend ae A oa ee pea Se (fp q.dt 
u(t) =— ( t+P t+P ; 
P(t +P) pot +P) — Dé) Go (?) ges NaS - bb py dt 
Setzt man 
Z = Max |u(d)|, 
0<st<P 


so folgt, wenn man fiir ¢ eine Stelle mit w(t) =Z einsetzt, aus der vorigen Glei- 
chung sofort eine Ungleichung der Form 


CK pie Zs 


B| 


also 
L(V LZ) 


und damit (35), wenn man (8) beachtet. 
Satz 5 und Satz 6 enthalten Satz 2. 
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Singular Integral Equations of Volterra Lype 
and the Finite Part of Divergent Integrals 


PL, BUTZER 


Communicated by A. ERDELYI 


1. Introduction 
In an earlier paper by the writer [2] the calculus of MrkUSINSKI was applied 
to non-singular integral equations of the convolution type. The object of the 
present essay is to discuss and solve by operational methods singular integral 
equations of the types 


(I) =€%) (¢> 0), 


0 
(11) FP { fdut+ei)=i) (#>0), 


of the first and second kind respectively for integral and non-integral values of 
a=11. In the more common case «<1, both integrals exist in the ordinary sense 
provided for example f is a bounded function. For «=1 an interpretation of 
the divergent integrals may be sought by Caucuy’s principal values for «=14 
and HapaMmarb’s finite parts for «>1. Equations of this type occur in aero- 
dynamics, as for example the Prandtl equation in airfoil theory. 


According to HADAMARD [7] we have the definition 


ts : 
Iv) =. [ _f(u) ep 
(1.1) EP ics Tau) ae | (ees Q(e)], 
0 

where Q(¢) is the unique linear combination of negative powers of ¢ and of log ¢ 
which must be subtracted in order that the limit may exist. To find the explicit 
expression of Q(¢) for a particular value of «, the integral and non-integral cases 
of « must be treated separately. Thus, if f is absolutely continuous, 


(1.2) re fg i) Mi 
and ; 
FP Pee) SEO 
firs (t—u)3 ta Vt i Vi—u ou 


0 


+ This problem was first stated in a colloquium lecture held by the author at the 
Technische Hogeschool, Eindhoven, Holland, May 23, 1958. 
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In the first formula the “finite part’’ FP can be replaced by the Cauchy principal 
value. The evaluation of the finite part? becomes increasingly laborious for 
higher values of «, and even after such steps have been performed the equation 
in question remains to be solved. Operational methods such as the Laplace 
transform may be suggested for example. 


In the following we proceed along different lines and define an operator 
{F P 1/t*} for integral and non-integral values of «=1 for the half-line t>0. 
This will be an operator in the sense of MixusiINsxtr [10], [5], having properties 
not unlike the Laplace transform of the distribution [4/f*] of ScHwartz [12], 
defined for the half-line 20. A generalized version of the operational calculus 
of MIKUSINSKI is then introduced and applied to solve the integral equations. 
Later in Section 3 a solution of the general equation (I) is given for integral and 
non-integral «=1, while particular equations of type (II) are solved in Section 4. 

Proceeding to summarize the results required from MIkusINsKI’s calculus, 
let C be the class of functions of ¢ which are continuous for ¢=20. Throughout 
this paper f or {f(t)} denotes the function per se and f(t) the value at ¢t. The 
addition of two functions / and g in C is defined in the usual manner and the 
product fg as the convolution 


(1.3) je =O} {60} ={ J He —w) e(w) dup. 


This defines a commutative ring with no unit element. According to a theorem 
of T1tcHMARCH this ring has no divisors of zero. It can therefore be extended 
to a quotient field, the field K of operators of MIKUSINSKI whose elements //g, 
g=-i0y, f. gee, are called operators. This field contains the complex numbers, 
the class C, the class L of functions which are Lebesgue integrable over (0, #) 
for every t>0 and also contains a unit element 1 which is the element correspond- 
ing to the Dirac delta function 6(t). The integral operator /={1} with the 


property 
HHO} ={ St) au} 


is also contained in K. More generally 


(1.4) ro=| [ Co" juan, WOM om 


n! 
0 


Furthermore K contains the differential operators s =1/l. If the (7 — 1)" deriva- 
tive ("— is absolutely continuous for ¢= 0 and f/”) = Wie (t)}, it can be shown that 


n—1 . . 
Sie = { ae > st a) (0) : 
4=0 


For details to the above the reader is referred to MIkusINsKI [10] or 
ERDELYI [5]. 


2 G. TempLe [15] considered the finite part of divergent integrals from the point 
of view of “weak” functions. F. J. BurEAu [1] has given a detailed discussion of the 
notion of finite part of divergent integrals, also emphasizing the historical aspect. 
See also LiGHTHILL [J8]. 
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2. The Operator {FP 1/t*}, a2 1 
Proceeding from these known facts, we consider functions which are integrable 
over every interval (e, ¢), e>0, but not over (0, ¢). For special functions of this 
kind, e.g. for 1/f*, «1, we define the operator {F P 1/f}. Let [«] denote the 
greatest integer <a. 
Definition 2.1. For «=1 we define? {F P 1/t*} by means of the convolution 
product 


t 
[a] hie, (tu) du 
(2.1) l (rp se}—{FP a at 
0 


provided the latter integral exists as HADAMARD’s finite part defined in (1.1). 
I’({«]) is the usual Gamma function with argument [«]. 

On employing the binomial theorem, we may write the right-hand side of 
(2.1) in the form 


a | aie ah : (= \F P / — | 


0 


The case when « takes integral values and the case for non-integral values 
are treated separately. In the first of the two cases we have: 


Lemma 2.1. If «=k, k any positive integer, then 


{F P1j}=s*heK, 


where 

h = (1/P(k)) {C, + (— 1)* * log } EL, 
and 

h—2 
ay (ea pre ahs 
ce Far ear | i ) fame 
Note that (.=0, 6,.-=—1, G;=3, G,=—3, 6(,=28,..., where the C,’s alternate 
mM Sign. 
To establish Lemma 2.1 one may write (2.2) as 
pee | ean ene he Res k-1 
Davee een oy aes ead | (4) du 
09) Hay QC a ee fas easyer fi 
0 

and since : 

t t 

du d pfi-« 
2.4 FP {* =1 Eis 
(2.4) ae eet | ea OE, 


the expression (2.3) is readily seen to equal h. Obviously {F P 4/t*t=s*h does 
not represent a function but is an operator in the sense of MikusINsk1. If s*h 
had been a function, it would have been 


ve 6” (t — u) h(u) du =h” (2). 


: This notation may perhaps be a poor choice as {F P 1/t*} is an operator and not 
a function, but we prefer its use here in comparison with the function {1/*} for « <1. 
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For the evaluations necessary in Section 3 we proceed to identify the operator 
{F P 1/t*} not with the quotient h//* but with an equivalent quotient of two dif- 
ferent functions of C. 


For this purpose we employ the function‘ ([6] p. 21 7) 


(2.5) RO ream 


which is a finite and continuous function of ¢. 
Then again use shall be made of the identity ([16] p. 179) 


(2.6) Flog + —y]o[(¢—u) edu —1 (¢> 0) 


where y = —I’(4) =0.577215 67... is EULER’s constant and & is an arbitrary 
constant. 

Multiplying both numerator and denominator of the quotient h/l* by {v(te*)} 
with 


(2.7) Gi (—1) ly 

we obtain the required equivalent form for {fF P 4/¢*, namely 

(2.8) UE Ger a ogee); (1) e 
D(k) Ik fy (te®)} I’ (k) lF{v (te8)} 


using the identity (2.6). 

We now give another definition of {/ P 1/t*} for integral «. 

Definition 2.2. 

LE BA = set 
—s k (=2)* Be 

fF Pap ty) — eee ee (EEN bg oh ta 
It is readily seen that this definition is equivalent to that of Lemma 2.1, often 
being of wider application however. It shows the close connection with the 
distributional calculus of SCcHwARTZ® defined over the half-line ‘= 0 ([12], voles 
p..42). 

Next we consider the second case where « is any positive real number +1, 
Ga 

Lemma 2.2. If «=k+ 6, where 0<B<1 and k ts integral valued, then 

= 1)kP(1—f) sh+8-4 

Fe (B+ 2)... (B+R—1) 
Proof. By hypothesis, we may write (2.2) in the form 


aan {tae (FSSA) F? fotal ate S(T] 


4 This function was first used by setneleehe in his theory of convolutions, a theory 
which may perhaps be regarded as a predecessor to the theory of MIKUSINsKI (see 
TEMPLE [14)). 

5 Also treated in SAUER [11] p. 242. For the actual connection between the theory 
of Scuwartz and that of MIKUSINSKI, see FENYO [17]. 


(2.9) 


(2.10) AE CK. 


14* 
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Making use of the identity 


seen ie Ltt (—yEr(A) POLO 
eerste i) nee 


and the property ([{10] p. 101) 
(2.12) {h=ru—pe* <1), 


the product of the operator (2.11) with ik takes on the form (2.10), establishing 
the lemma. 
Definition 2.3. We define by induction for | Rs WE, Perey 


B(B+1) (B+2)...(B+h—1) 


Pip oe ror 


This definition is equivalent to that given by Lemma 2.2, provided of course 
the property (2.12) is used. 
It is now possible to summarize the results obtained in the following table. 


—{F Pj te}, G84. 


Table 1 
1) {mnt ea, ta 
2) (= r(a4a) ps F OE) A>0 
3) (}=t=— 
4) at =I(i—_py? = — Cet han 
5° {pt} —stlog 4} 
6) fe Pa} = pags (i+ (— 1)! tog 2), R152, 3, 
7) {FP ag} = USS O<B<1 
8) {EP arg = | 1) apy ay b= 12,3, a Spee 


The first three elements belong to the space C, the fourth to LZ, while the last 
four belong to the field K of operators. Entries 1) to 4) have already been given 
by MixusiXsk1 ([//] p. 328) and are stated for the sake of comparison. 

Having introduced the operator {F IP {t*\; “«=1 into the Mikusirski calculus, 
it remains to investigate the rules of operation of this operator with other func- 
tions in C. Just as the convolution product of /"*1 with / was discussed in (1.4), 
the product of {F P 1/d*} with {f(é)} can be treated in its turn. 


® Compare also Lavorne [9] concerning a similar result in the distributional 
calculus; also O’KEEFFE [19, 20]. 
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Theorem 2.1. If «21 and f\@1-)) is absolutely continuous for t=O and 
f(0.) =/'(0) =--- = f1-*)(0) = 0, then 
an. 


EPA HO} =| 
Proof. Firstly we establish the result for non-integral values of c<=1. Suppose 
a=k+B, 0<fP<1, k any integer =1. If /*~” is absolutely continuous, then 
it is known (cf. [8] p. 665) 


t 
Wap galat Lie weld (0) speak 
sesh re aime (FP {du 
0 


REB—1 tbl kt B—1 (t—u)e+B-1 
(2.13) Me 
a #(0) oe (0) 
EB e PUP a (bP 1) (Reb — 2) ete 
feel) (—1)* 1.0) ay, 


(b+ B—1)...(B+1) BP ” RFB—1)...B+)B J (uy? 


On the other hand, according to Definition 2.3, if / is absolutely continuous, it 
follows that 
faa pagar 4s serene aap eeneginh opr 


Repeating this operation 
Eners \°) f(0) 
{F Pag} app Past t apn es 
Pilon tiAR 2) it (a  ( 
: wero ob ah eae OG 


on the assumption that /~” is absolutely continuous. A comparison ®* of (2.13) 
with (2.14) completes the ae for the non-integral case. 


(2.14) 


A comparison of formula (1.2) with 


(2.15) {F P 1/t} f = (0) {log ¢} + f’ {log 2} 


in case / is absolutely continuous extends the result for the integral value k = 1. 
If #*-” is absolutely continuous, we proceed similarly as in the non-integral 
case. Comparing the expressions obtained by k-times repeated application of 
formula (2.15) with that of (1.2) likewise repeated k times, the desired result 
follows. 
3. Solutions of the singular Integral Equation (1) 

Solutions of (I) in case « is integral and in case « is non-integral can be given. 

Theorem 3.1. Let k be any integer =1 and let f*—” be absolutely continuous 
for t=0 and gE. In the case k=, and {(0)=0, 10774 S101 nen = 23) the 


6a Here the vanishing of the derivatives of / at zero is needed as the right term of 
the equation in Theorem 2.1 is only defined (as follows from (2.13)) if the given deriva- 
tives vanish at zero. 
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solution of 


(3.4) BM bees iat (t> 0) 
is given by : 
(3.2) i (0) = (=r fase | fete) ele a é]dw|du 


where v(t) is defined in (2.5) and E im (2.7). 


Proof. According to Theorem 2.1, equation (3.1) takes on the operational 
form 


{FP AEH EOS = {85 
By Lemma 2.1 and the relation (2.8) it follows that 
k 


(3.3) fas% = (1h Pa)? folte)} g. 


The solution then is 
f(t) =(—1)*L B) fat fdty. ay. 1F glw )v[(t—- 1— wv) & | dw, 


which is equivalent to (3.2). This solution satisfies the given equation (3.1). 


Remark. From equation (3.3) it is obvious that for k =1, the solution of (3.1) 
is given by 


fj =-# Jet )v[(t —w) e’]du 


(since €=y here), provided the indefinite integral in the square brackets is 
absolutely continuous, which is the case if g is absolutely continuous. Under 
this hypothesis and applying a theorem on the differentiation of a convolution 
product (cf. Dortscx [3] p. 118) we may write the solution in the form 


f(t) = — g(0) v(t e’) =f gu )v[(¢—w) e’| du. 
For k =2, since & =1+-y, the solution is 
t 
=f g(u) v[(t¢—u) e+” ] du, 
0 
provided g€E and /(0) = 
Theorem 3.2. Leta =k+f,0<B<1, k any integer =1. I} f?-» is absolutely 


continuous for t=0 and gE L, and f(0)=f'(0)= ++: =f*—)(0) =0, the solution of 
t 
4 BP | oes 
ie h (t—u)k+B ee ae 
as given by ‘ 


t 


1) = (1) b+ B14) SABA f wh" g(a) du. 


A 
0 
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Proof. Since /“~» is absolutely continuous, we apply Theorem 2.1 and it 
remains to solve 
{EP 4/PtP\ p= g. 
By Lemma 2.2 
[= se 1)* '(k+B) Pre l y 
P(1—B) P(B) 


and since |J'(1—) I’(8)] += (1/z) sin Ba, the theorem follows. 

The different integral equations with the corresponding solutions obtained 
from Theorem 3.1 and 3.2 for particular choices of k and « are listed below. The 
first integral with solution is the well-known Abel equation (see also [2]), given 
for the sake of comparison. The table can easily be extended for different and 
larger, integral and non-integral values of «. 


Table 2 
Integral Equation Solution 
t t 
iC tal fig, 
") ie du= 3) H(@) I ak (t—u)3 ‘ 
t 
f(u) = Mis v —wu) ev 
2) ee) Lo au = e(e) i io (Uw edu 
: eg =1 few) 
3) Eel, (a, du = g(t) HO) == ai i (any du 
0 
t t 
f(u) 
4 Ie yg GU = g(t) f() = | g(v) v[(t¢—w ety] du 
Jw Z J 
t t 
f(u) 3 4 
5) FP { 1 au—g() He) = | ew) g(u) de 
fy x! 
6) FP [ fo du=el f(t) = Se) du fate v{(t—u) e*7] du 
: (t—u)3 ' i 


4. Special singular Integral Equations of the Second Kind 


Let us consider a Prandtl-type equation used in airfoil theory, one of the 
Volterra type instead of the usual one with range of integration over a fixed 
finite interval (—y, 7) (e.g. ELtiott [4], SOHNGEN [13]). 


Theorem 4.1. Let f be absolutely continuous for t=0 and gC L. The solution of 


(4.1) aaah MO) dutg=1 (>0), 
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is given by 
t t 
(4.2) f(0) =f glu) gi —w) du + f(0) |1—f au) au, 
0 
where q(t) is defined by the integral along the line Rp =c 
c+10o 
1 ept 


c being a suitable, large positive quantity. 
Proof. We may write (4.1) in the operational form as 

(FP ah f +e=t, 
and multiplying both sides of this equation by {v(te”)} we have 
(4.4) {o(te”)} {log t} sf’ = (f — 8) {u(te”)}. 
Making use of the identity (2.6) and the fact that 

sf’ =s[sf—f(0)] =s?f — sf(0) 
which holds since / is absolutely continuous, we can write (4.4) as 
[(t-r ol = geer 270), 


where for simplicity we set {v(te”)} =v. Hence 


Be 1 f(0) 
(4.5) LSgeeg a Ogee 
But for v€C, it is a known fact ([9], []) that 
ae ~ a len te Le ~ __ 4\" an 
(4.6) ee: = 21 4)"-1y", eu 2! tu" 


The problem now is to evaluate? these sums explicitly. With v(t e”) as defined 
in (2.5) it can readily be established by induction that the 2“ convolution power 
v" (te”) is given by 

fo, 0) 
ee tide: eve yn—-l jutn—l1 
(4.7) v" (te”) Tay ea du 


for nm =1, 2,3,.... (This can also be established indirectly using Laplace trans- 
forms and formulas (18) and (19), p. 222 of [6].) Now in the calculus of residues 
it is a known fact that (e.g. Dorrscu [3] p. 213) 


c+t0o 
8 Onics ept = ia ie. 
a) OX AINA a Dak ap P(u+n) ’ 


d The writer is greatly indebted to A. Erp&éLy1 for generously providing the key 
to this argument and the contour integral representation for q (t) 
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where c is a suitable positive quantity. Formally then we proceed as follows. 
Inserting (4.8) in (4.7) we obtain 


eo) eS ‘ c+i oo sf 
(—1)"-1 0" (t e”) Cis Pet a dust [ d 
2! = 5 P'(n) Bi ah alae Piet? p 
oO c+ico te Se 
(=4) aes ee era a) du 
n=l c—ico 
= pt me a 
a Fe io [ge 5 dade du 
Cc— 100 
c+ico 
Lert eet ¥p\ Ub 
271 fe ‘fp e p oa 
cC—10 


e. ert 
one iz (y— log p) — ie: 


The interchange of the order of integration here effected is justified; also term- 
by-term integration as 


Sf 


is convergent for ~ sufficiently large. Hence if c is sufficiently large we have 
obtained a contour integral representation for the first sum in (4.6). Regarding 
the second sum in (4.6) we have 


co 


1 
du— >’ y” (log y—y)” (|p| ="), 


n=1 


en yj Ne! 
ean | 


co 


3B (<1) 1 — SF (110 1 fa}, 


n=0 n=1 
where 1 is the unit element in K. In view of equation (4.5) our theorem now 
follows. 
The method is not restricted to this equation but can be extended to cover 
singular integral equations which reduce to algebraic equations of the form 


(FP 137 O}+e =F, 


k integral, provided’ /“~» is absolutely continuous, and /(0)= --- = f®-)(0) = 0. 
In the non-integral case we consider only one special case: 


(4.9) siete Le! A) snes agian (t>0). 


In the operational form 


with solution 
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Hence ynle 


1 1 SN i) ee late 
NO sgt 21> Unter 


To determine this sum-function the sum can be broken up into two parts cor- 
responding to terms with odd and even powers of |/¢ respectively. Those with 
even powers give 

(V3 os (2 Oa ae 

ray + re) t PG) 


Those with odd powers 


a CRAG 28 W'S th(Vé 
Te re Tw ray eV A, 


where 


Now the derivative h’(t) satisfies h'(t)=2-+2th(t), whence it follows that 
t 
h(t) =2e [e-” du. Then 


® 2o 2 Fi t : x 
ye (— 1)" a ee pag ote ies du. 
n=1 


L'(n/2) Va Vx 
The solution f(¢) of (4.9) then takes the form 
f(t) =e (1 —erf yt), 


where erf¢ is the usual error function. Indeed, this function is a solution to 
the original equation (4.9). 
The more general problem® of treating singular integral equations of the type 


FP { 4M 1(¢—wdu=g(d) (¢>0), 


say, with kernel of the form 7 (¢)/é, «21, where 7 and / are bounded functions 
for +20 for example, instead of the simpler form 1/f*, «21, remains open. 
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Introduction 
Considered are asymptotic solutions for differential equations of the forms 


yy" — (x — 1) 7a (%, AM) y =, (1a) 
yy’ — Mx (% — 1)? gi (*, a7) y = 0. (1b) 


Here x and y are the dependent and independent variables respectively, primes 
indicate differentiations with respect to x, A is a parameter which is allowed to 
assume large values, g,(x,A7) is a one-valued analytic function of x and 47 
which, at least for large values of 4, is free of singularities and zeros within the 
region of the complex x-plane under consideration. 
The point x=1 is a singular point; in its vicinity the solutions can be re- 
resented in the form 
i yp= (x — 1)" B(x 1), (2a) 
vine (a) a7) (2b) 


except if @,— 0. is an integer. P,, A and later similar expressions stand for 
functions which are regular at x =1 and can therefore be represented by power 
series. The exponents @, and 0, may be defined so that 


Re 0, = Re 0. 


The purpose of the present paper is to derive asymptotic representations (for 
large values of A) for the solutions (2a) and (2b). The classical formulae for 
asymptotic representations always give two linearly independent expressions for 
the exact solutions of (1). At x=1 they even have a form corresponding to (2). 
However, it is not justified to assume that the asymptotic expression possessing 
at x1 a behavior like (2b) represents the exact solution defined by (2b). In 
fact, for a differential equation with transition point, (1b), such a coordination 
would be wrong. 

The expressions (2) give more than merely the behavior at x=1. They 
imply that, if y; and y,, are continued around a path in the complex x-plane 
which encircles x =1 in the positive sense, their original values will be multiplied 
by factors e**'* and e?*'% respectively. This property makes it possible to 
characterize the solutions y, and y,; at a distance from the singular point, and 
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on this basis the desired representation for y,, will be obtained. The present 
paper gives the main features of the argument leading to this result. More details 
can be found in [/]. 


_ 1. Region without Transition Point 
8a(x, 47) = [gr(*, 4) + oA]? 
Of the two branches of g,, that one will be used for which 


Let 


9 ~—<are[Ag.(1,A4%)| Sa. 


The argument 47 of g, and other expressions will be omitted in subsequent 
formulae. The exponents of the regular singular point at x=1 are given by 


(3) 
Let jag 
h(x) = ee — 1) g2(€) dé. (4) 


To define the argument of the function h(x), we require a condition for the 
lower limit of the integral, vzz. 


Lint [arg(& —1)] =72. (4a) 
Then two linearly independent asymptotic solutions of (1a) are given by 
Yi2=h' (x) texp(F Ah(x)). (5) 


They are asymptotic in the following sense: Assume that at a point x =0 one 
has an exact solution y for which one pair of the following relations holds 


y (x) = ¥,(*) (14+0(4%)), 
y' (x) = Yy (x) (1 +0(A)) 
y (x) = ¥,(x) (1+0(4%)), 
y' (x) = Ye (x) (1+0(A)). 


Then these relations hold for all points of the x-plane which can be reached 
along paths for which the real parts of the exponential function in (5) does not 
decrease. 06 will subsequently be called the “point of matching’. Equation (4) 
may be written in the form 


(6a) 


Or 


(6b) 


h(x) = h*(x) + h(x), (7) 
where 4 
We(a) =e (4) PE 1) tae 
i (7a) 
= ¢'* g,(4) [In(x —1) —4z] 
and 3 
a), = oe I (es(é) ee) 1) GAs. (7b) 


Then for the vicinity of x =1 
i .= e—4** (0 (1)—# eX? (—all)iath())] (x — 4)8=482(1) (1 + O(x% = 1)) ‘ (8) 
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Thus at points near to x =1 the asymptotic expressions have the form (2) 
The function h, and accordingly the regions of validity of the asymptotic expres- 
sions (5), depend upon the function g,. It can however be shown that h has 
the same topological character as the function /*. Therefore future discussion 
of the regions of validity can be restricted to h*. According to (7a), h* has a 
logarithmic branch point at x =1. To obtain a survey of all Riemann sheets, 
we consider h* as a function of 
a variable 7 defined by 


imaginary 
ax/s 


| | 
| | | | 
ee | 
x A_,(h=2119, (1) et 
This relation maps the entire x- ae | | | — 
plane into a strip of the 7-plane } | | | 
of width 2z in the imaginary Wi | | 
direction. Fig.1 shows the re- VE 
gion considered, in the x-plane. map of the v7 PE 5. 
boundary of Tcariee. 
the region C | Seat —OXIS. 
S ~ AW Gt ; 
Ss / 
Ne § Geil Ke Wey: 
SS complex tL 1 ea ae Bk oe 
< x-plane SE eae 
C ‘ ~ aes 
\ : 
/ T| 
real axis Va = 
a = | 
a ean — 
. Ap (h=-4I0 lG> (1) | | 
NX 
@ | 
Fig. 1. Region C in the complex x-plane Fig. 2. Region C in the complex 7-plane with lines Re Ah = const; 


A go(1) real 


The corresponding representation in the 7-plane is shown in Fig. 2. (The straight 
dash dot lines in Fig. 2 and in the subsequent figures are the boundaries of such 
a strip; they correspond to a branch cut in the x-plane extending along the real 
axis from x =1 to negative infinity.) The point x =1 is mapped onto + o in 
the 7-plane; x= oo, to — oo. Travelling along a closed path around x =1 in 
the positive sense, one arrives at a point 7 which arises from the starting point 
by substracting 227; e.g., the origin of the x-plane is mapped onto the points 
A_,, A_,, Ag, Ay, etc. When Ay is the starting point, A, is reached by travelling 
around x =1 once in the positive sense. The region C in which the asymptotic 
solutions are considered may be bounded (Fig. 1); the boundary in the 7-plane 
is shown in Fig. 2 but omitted in all subsequent Figures. Thus the shaded part 
of Fig. 1 will map into the shaded part of Fig. 2. 


The solid vertical lines in Fig. 2 are lines of constant real part of Ak* drawn 
under the assumption that Ag.(1) is real. Fig. 3 is a corresponding sketch for 
a complex value of Ago (1). 

The definitions of g. and / are such that Re AA tends to + coas one approaches 
the point x =1, or as Rey tends to + co for a bounded imaginary part of 7. 
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Therefore in the expression Y, in (5), the real part of the argument of the ex- 
ponential function increases if the real part of 7 decreases. Consequently the 
point of matching will be chosen at x =b =1, i.e., at oo to the right in Figs. 2 
and 3. Then all points of the region C in Figs. 2 and 3 can be reached by ad- 
missible paths, and in the entire region C there is an exact solution yg for which 


1S valid the relation Vo (x) a vi (x) (1 - O (A+)) ' 


| imaginary 


real N - axis 


Fig. 3. Complex 7-plane with lines Re 4 = const; Fig. 4. Region of validity for Y, 
A g.(1) complex 


Since the exact solution y, becomes small in comparison to y;; as one approaches 
x =1, it can be shown without difficulty that, except for a factor, y, is identical 
with y,. This factor can be found immediately from (8). Thus 


y= eh ™ (1) —4 e772 (—82(1)in+h,(1)) vane ans 0 (A>) : (10) 


In order to obtain a large region of validity for the asymptotic solution Y,, 
the point of matching }, must be taken far to the left in the complex 7-plane. 
But then it is impossible to decide on the basis of the appearance of the expression 
Y, at the point x =1 whether the exact solution represented is a multiple of yy 
or a linear combination of y, and y,,. 

Here the idea described in the introduction is helpful. Let us denote by y, 
the exact solution which matches Y, at the point x =0. Then there exists a 
region extending over several sheets of the x-plane, for which 


ya(x) = Yo(x) (1 +0(4%)). (11) 


The shaded part of Fig. 4 is the region of validity of this expression. Let us 
ascribe the numbers —1,0,1,... to different sheets of the complex x-plane in 
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such a way that one arrives upon the sheet with the next higher number by a 
complete circuit around « =4 in the positive sense. Let us differentiate between 
values of x lying in different sheets by superscripts in parenthesis. Then the 
form of Y, given in (8) shows that 


¥, (0) — Ermer, (4M), (12) 
Furthermore, because of the definition of y, and y;,;, one has 
yy (XD) = oP FA y, (4)), (12a) 
zp (XO) = eB *4& grey (249?) . (12b) 
y, can be represented as a linear combination of y, and yy 
Ya(%) = Cy Vy (%) + Or Vir (*)- 


This relation is valid at points x as well as at point «7. Using (11) and (12) 
one thus obtains 


Cy Vy (4) + Ory Nyy (x) = Y, (x4) (4 +O (A>)) ) 
Cy Yq (HO) CP FAL Coy Yyy (x4)) €? 4402 — Ye (x10) €?™*%2 (4 + O(A4)). 


Using the asymptotic representation for y,, (10), and the explicit expressions 
for 0, and 05, (3), one finds from the above equations 


Y= em? 25 (41)? et4(—8(L) tz +h,(1)) x 


, : 4 Ki) ¥ilb)  e-Bataety 
Yn {i +O0t pale, Y,(x) Y,(b) aiastista 


2. Differential Equations with a Regular Singular Point 
and a Transition Point 

For the differential equation with a transition point (1b), the procedure is 
in principle the same. One first finds overlapping Riemann sheets in which the 
asymptotic solutions are valid. Then the relations connecting the solutions in 
these two sheets will give the desired result. A slight complication arises from 
the more complicated structure of the Riemann sheets and the fact that the 
asymptotic solutions appear here in a form which is less tractable than in the 
previous case. The main points of the proofs will be sketched and then the 
results will be given. For further details see [7]. Let 


h(x) =e f (1) ga (@) ae. (13) 
For the lower limit x =0 of the integral, let 
lim [arg(f —1)] = iz, 
lim [arg (é)] =0. 
Two linearly independent asymptotic solutions are now defined by 


V, = eiti2+2in/s g3 eho h'(x)~3 Ho) (Ah (x) ein) (14a) 
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and 


Y, = em yh lie h(x)* h'(x) 4 H® (A h(x) ef”). (14b) 
The following asymptotic representations for the Hankel functions will be needed: 
ea) le eae Fe AOR) =m eater 20 (15a) 
and 


Bie ao \2 eI z~4[ cos (z—a/4) +O (z) cosh (Imz)], m<argz<3a. (15b) 
Using (15a) one finds 

Y,=h' (x)rFe""(1 + O(A)), —3n<argih<in, (16a) 

Y,=h'(x)-te-*"(1+ O(A4)), —8a<argah<2a. (16b) 
The last expressions are identical with the expressions (5); but they are less 
general than (14) since their region of validity does not include the origin and 
the lines argAh=}a and arg Ah =z respectively. The expressions (14) are 
approximations of exact solution in the following sense. Let us assume that 


exact solutions y, and yg, are chosen such that at a point «=0 the following 
relations are fulfilled: 


Va (b) = ¥.(0); — ¥p(b) = ¥, (2), 


9), (Oy ¥,, (0); “ya (0) = Y; (6): 
Then one has 


Meg = oe mm O(A-8) | 2(x)|* 


h'(x) 


“4M, (My + |Al?|A(x)|2)~het44* (47) 


for all points that can be reached from x =0 along a path for which Re -Ah 
does not decrease; the positive sign is to be used for y, the negative sign for yg, 
M, and M, are two positive constants, independent of A and x. In the error term 
the following restriction in the argument of AA must be observed 


—iasargd<hsia TOYO (17a) 
—3aSarghihsia for y,. (17b) 


Moreover, admissible paths will not go through a point for which h’=0, unless 
simultaneously 4 =0. 


The function / is written in the form 


h=h*+h,, (18) 
where 
h*(x) = ef g, (4) fet(E — 1) 28, (19a) 
jn x) =e f €4(€ — 1) [gp (6) — go(t)]a€. (19) 
0 


There is a wide class of functions # having regions of validity with the same 
topological structure as that of the function h*. To problems of this kind our 
future considerations are restricted. 
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To show the behavior of the asymptotic expressions at the point *=1 we 
write / in the form 


as e*|go(t 1) f( E—1)1d& + ha(x)| = e'*[g.(1) [In(x — 1) — 7] + he(x)] (20a) 
where 


hy (x) = fe-ay 1 [54 go(€) — &o(1)] 46. (20b) 
hy is regular at x =1. Then one finds from (16) for the vicinity of point x =1 


Vg ei gy (1) BoA Leath mA] (4p — 4) FA6 [4 +O(A4)] [1+ O(x—1)]. (24) 


a, 


imaginary axis 


Fig. 5. Lines Re 7 = const in the complex x-plane for a differential equation with transition point; / g,(1) real 


The expression h*, (19a), can be integrated 


h*(x) =e gy(1) [2 Vx +n ies ial. 


On the basis of this formula, or by qualitative considerations, one can carry out 
the discussion of the regions of validity. Assume first that Ag, (4) is real. A sketch 
of Re Ah* in the x-plane is given in Fig. 5. In the vicinity of « =1 the dominant 
term in h* is In(x—1), the lines of constant real part of AA* are circles, at a 
greater distance they are closed curves. In the vicinity of x =0 the expression 
(19a) for h* can be integrated to give 


W(x) = g9(1) @) ab + 
Thus h* maps itself into two sheets of the complex x-plane. The lines of constant 
real part have a branch point at x =0, the branches form angles +42 and a 
with the positive real axis. For very large values of x the term x? prevails in 
the expression for h*; the curves of constant real part become parabolic. These 
considerations lead to Fig. 5. As in the previous section a representation in 
terms of the variable , (9), is chosen in order to give a better survey of all 
Riemann sheets. However, this variable does not remove the branching of the 
Riemann surface at x =0; therefore branch cuts, starting at the maps of point 
«x =0, must be introduced. They will be chosen to follow one of the branches of 
the lines Re A =const which run through the maps of x =0; no need will arise 
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for admissible paths to cross these barriers. The mapping from the x-plane to 
the 7-plane is conformal for *=0; therefore the character of the branching 
remains unchanged. Expressing h* in term of 7, one finds: for x—> oe, i.e., for 


* ht ~ — 2g9(1) 7” 


for x1, 1.e., for Rey oo 


bro en (hyn: 


a) 
= 
= 
2 
= 


Fig. 6. Lines Re Ah=const in the complex 7-plane for a Fig. 7. Lines Re 2h =const in the complex 7-plane for a 
differential equation with transition point; Ag.(1) real differential equation with transition point; Ag,(1) complex 


On the basis of these remarks the equivalent to Fig. 5 in the 7-plane has been 
drawn in Fig. 6. The points A_,... A, are maps of the origin of the x-plane. 
The line connecting them is the closed line Re A =const of Fig. 5 which starts 
and ends at the origin. The wavy line to the right of the points A in Fig. 6 is 
the map of a closed curve around x =1. As one goes to the right such lines will 
approach vertical straight lines. This gives the connection with Fig. 2. The 
curves to the left of the row of points A represent curves Re Ah =const of the 
x-plane which lie outside of the closed curve Re Ah =const going through the 
origin. 

Similar discussion can be carried out if Ag,(1) has a complex value. The 
result is shown in Fig. 7. 

In studying the regions of validity one must choose a starting point for the 


various paths, and then determine which regions of the 7-plane can be reached 
Gy 
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by paths for which Re(+-A/) (the positive sign for Y,, the negative sign for Y;) 
never decreases. If one branch of a certain line Re Ai =const can be reached 
by an admissible path, all of its points can be reached. One may, therefore, 
consider the lines Re 4A const as contour lines, and determine the region in 
which an asymptotic solution is valid, by proceeding from one line Re Ah =const 
to the next one, as a rising flood would proceed from contour line to contour 
line. In such a procedure the origin of the x-plane needs special attention. In 
its vicinity # =const + const x}, and this character is carried over into the 7-plane. 
The first constant is unessential as far as the character of the contour lines is 
concerned. The function # extends through two sheets because of the expression 
2 x. In different sheets the direction of 
2 increasing Re Ah is different. 
The vicinity of such a point is shown 
E in Fig. 8, and it is assumed that the 
y direction of increasing Re Ah is known 
at a certain point P of one sheet. Then 
it can be determined at other points 
by consideration of continuity. In this 
process one finds quite automatically 
that a branch cut must be introduced 
along Ask. One notices that once the 
level of AB or AC is exceeded, the region 
ie in of validity extends immediately along 
CE AK. Points of the line A& can be reached 
Fig. 8. Contour lines in the vicinity of a transition point Only by means of the path AK, 1.e., along 
a path going through point A. It was 
mentioned before that an admissible path must not go through a point for which 
h’=0 unless h=0 simultaneously. In our example, paths going through point 
Ay in Figs. 6 or 7 are allowed, but all other points “‘A’’ are forbidden; they are 
marked by double circles. For a branch cut starting at one of the ‘‘forbidden”’ 
points “A” the asymptotic solutions are not valid; actually even the vicinity 
of such branch cuts must be excluded. These branch cuts are indicated in the 
subsequent figures by double lines. A branch cut for which an asymptotic 
solution is valid is indicated by a single heavy line. On the basis of these remarks, 
the reader will recognize how the limits of validity shown for Y, in Fig. 9 and for 
Y, in Fig. 10 arise. In each case the point x= indicates the point of matching. 
In (16) limits for arg 2h were imposed. The branch cuts starting in Figs. 9 and 10 
at point A, correspond to these limits. 


On account of its behavior at x =1, (21), the exact solution represented by 
Y; can be immediately identified with a multiple of y,. One thus obtains 


yy weit g, (4)4 et ladin—natty] (22) 


except for the vicinity of the double lines in Fig. 10. To characterize these lines 
in the x-plane one first determines the curve for which 
Bue 


arg Ah = aE 
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using paths which do not encircle x =1. This line starts at x =0 and extends to 
infinity. In the Riemann sheet for which h=0 at x =O this line represents a 
branch cut for # but not for the function Y,. In all other Riemann sheets this 
line represents the limit of validity of the asymptotic solutions. Naturally, for 
practical purposes the asymptotic solution is needed in only one Riemann sheet. 

Replacing the Bessel function in (22) by its asymptotic representation (15a), 
one finds 


yer gy (4) ete ea —())] h'(x)~#e7* (4 +0 (A) —ia<argih<n. (22a) 


oe a a Ap 

= =O) Sa eee a Ca 
Ap 

Fig. 9. Region of validity for Yy in the complex 7-plane Fig. 10, Region of validity for Yg in the complex 7-plane 


Here the curve arg Ah =3a and its counterpart in the other Riemann sheets 
represents the limit of validity in all sheets of the complex x-plane. The re- 
presentation thus obtained is identical with (8). Fig. 11 shows the limit of the 
validity of this equation in the complex x-plane. The middle part represents 
the sheet which contain Ay. The heavy line represents the limits of validity. 
By crossing the dash dot lines one arrives in other sheets, which are shown above 
and below. 
A representation of the same kind which is applicable for the vicinity of the 
line arg Ah =-+ 3a is found by means of (15b): 
yy = eb? 9 (1) e4 Leal ia—tald)] yy 
x 2e'7/4 h’(x)—4 [cos (Ah (x) i — $2) + O(A4) cosh (Re |A/|) ] (22b) 
an <arghh< 3a. 
The region of validity (as expressed by the condition for arg Ah) in the 7-plane 
can easily be seen in Fig. 10. The corresponding region in the complex x-plane 
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is the shaded part of Fig. 12. A representation similar to Fig. 11 showing three 
Riemann sheets separately has been adopted. 

The region of validity for Y, shown in Fig. 9 contains parts which originate 
from each other by a complete circuit around x =0 and x =1, @g., the points 
x—b and x=b.,. Their behavior can be used to determine a representation 
for y,,. For a path which encircles the points x =0 and the x =1 in a positive 
sense, one finds (see [1] for details) 


h’ (xe?) = eh (x), (23a) 


h(x e") = 2migy(1) —h(x). (23D) 


CA 


VL 
=/ arg Ah = 


a Ce a 


Az 
Fig. 11. Region of validity for Eq. (22a) Fig. 12. Region of validity for Eq. (22b) and (28c) 


Let us now consider the asymptotic solution for y,, (17), and insert into it the 
asymptotic expression (15a); then 


va (x) =h'(x)bet*)(140(04)) —Suc<argdh<}n. (24) 
This representation is valid in the region CA_, Ay E and in the region EA, A, F 
of Fig. 9. If point x is chosen in the region CA_,A,)E the point xe?” will lie in 
the region E A, AE . If one considers (24) as a representation for the first region, 
insertion of (23) yields the following expression for the second region: 

Vie eee Veen eam eset) 7! (enema) (1 +0 (A>) é (25) 


yz 18 now considered as a linear combination of y, and y,;: 


Vu = Cy Vy + Cy Vyy- (26) 
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Applying this equation to a point x and to a point x e?** reached by a path 
encircling x =0 and x =1, and inserting into the equations thus arising the above 
asymptotic representation, (24) and (25) for y,, one obtains 


Cy V(X) + Oy Vy (%) =’ (x) “Fe? (4 + O(A4)), 
CV; (x) ame eit Ge. (x) e202 pmi/2 ert Ags(1) pr (x)-3 ena +0(A>)). 


Fig. 13. Region of validity for Eq. (28a) Fig. 14. Region of validity for Eq. (28b) 


From this system of equations an expression for c; yy can be obtained and since 
an asymptotic expression for y,; is available in (22) one then obtains, after some 
manipulation, an asymptotic formula for y;; 


' a3 sie entailed 
sa imma SB Tbice ne + 2simh@aidgny Ost 
+ 0(A-8) |h(x)|* | h'(x)| “4 My [MS + | ]> | 2 [8] -# x (27) 
Dhog( 2) Cae TERY) hee 
x {le | 4 ant 


where /i,(%) and g(x) are one-valued functions of x. One has 
h,(*)=h(x) for —$aSargihsgn 
and 3 


hg(x) =h(x) for — $uSargihsen 


h,, and hg occur in the error terms only. 
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The coefficient of Y, in (27) can be expressed in terms of 0, and 0g. Equa- 
tion (27) is valid in the common region of validity of the expressions Y, and Y,, 
Figs. 9 and 10. 

Asymptotic representations for y,,; in terms of exponential functions may also 
be desirable. They are obtained by inserting into (27) different asymptotic 
representation of the forms (15) for Y, and Y,. Naturally, one must observe 
the regions of validity of these representations. In order to obtain overlapping 
regions of validity, the following three expressions are needed: 


yor a — oie gy (1) amt eile RGU 


x W(x) Her (4 +O(A%)) 4 


e274 gi(1)+21/2 


gh) (4 +0(47))} (28a) 


2 sinh (271 1g, (1)) 
—3a<arghih< 5a. 
The region of validity is shown in Fig. 13. 
Vry (x) = et m2 25 (1)2 e? (g2(1) i7—h,(1)) x 


xh (ayHe*)(4 +0(0)) 4 


e271 4g.(1)+71/2 
2 sinh (227 / g,(1)) 


e**)(1-40(24))) (28b) 
—$u<argah<3na. 
The region of validity is shown in Fig. 14. 
Vay (X) =e? ge (1)F eA eal) Fx —A) ((x))—4 [2sinh (277-0 g, (1) |x 
X hi (x) Be? 4 Aealt) +d g—AR() (4 + O(AA)) — eB treat) CAM) (4 + O(A4))} (28c) 
ou <argah<3n. 


The region of validity is the same as in Fig. 12. 
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Vom Defekt ausgehende Fehblerabschatzungen 
bei Differentialgleichungen 


JOHANN SCHRODER 


Vorgelegt von L. COLLATZ 


In einer vor kurzem erschienenen Arbeit [4]! wurden Fehlerabschatzungen 
fiir Funktionen u*(x) hergeleitet, welche auf einem #-dimensionalen Gebiet 8 
einer Differentialgleichun 

6 6 M{[u] =f (x, u) (0.1) 


und auf dessen Rand J” gegebenen Randbedingungen (oder Anfangsbedingungen 
oder ahnlichen Forderungen) 


Uful=y, G=1,20.) (0.2) 


gentigen. Zu diesen Abschatzungen bendtigt man zwei Funktionen w,(x) und 
u(x), die durch die Gleichungen 


M(m|=f(x,%), Ulal=v @=1,2,...) 


verkniipft sind. In vielen Fallen ist es médglich, solche Funktionen zu finden 
(siehe die Beispiele in I), in anderen jedoch kann dies schwierig sein, insbesondere 
bei komplizierterem Gebiet 8. Hier werden nun aus den in I bewiesenen Ergeb- 
nissen Fehlerabschatzungen hergeleitet, fiir die man nur eine Naherungsfunktion 
u (x) braucht und die von einer Abschatzung der Defekte 


—M[uo] + f(%,%), Gilt] — v1 
ausgehen. Dariiber hinaus enthalt diese Mitteilung auch Eindeutigkeitsaussagen. 
In Nr. 4 wird ein (als Vorschlag aufzufassender) Rechengang zur Ermittelung 
einer Naherung u(x) und zur anschheBenden Fehlerabschatzung angegeben und 
eine Randwertaufgabe bei einer gewodhnlichen Differentialgleichung danach 
gerechnet. 
1. Voraussetzungen und Bezeichnungen 


a) Wir verwenden die Bezeichnungen aus I, ohne sie noch einmal zu erlautern. 
Die in Nr. I, 2.1 genannten Voraussetzungen seien erftillt. AuBerdem gelte auf 


BA Gr eye Ve {Ul ee et (1.4) 


% bedeute die Menge der Funktionen u(x), welche (auf $) stetig sind und die 
von der Lésung w*(x) verlangten Differenzierbarkeitseigenschaften besitzen. 


1 [4] wird im folgenden — bei Formelnummern und dgl. — als I zitiert. Die Lek- 
tire von I wird vorausgesetzt. 
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Fiir jede feste Funktion w€% existiere eine Lésung h =h{u| der Randwert- 
ania Mih]=0, U[A]=U[w] (=1,2,...). (1.2) 

Der Einfachheit halber wollen wir ferner voraussetzen, daB die auf dem 
Gebietes) @: «eB, oli SyvSvi(o) 
erklarte Funktion f(x, y) und die auf (I, 1.4) 

: xeB, OSyGA*(z) 

definierte Majorante } (x, y) stetige erste Ableitungen nach y besitzen und (bei 
geeigneter Bezugsfunktion w (x)) 


~ 


0Sh(%,»)Si(xy') fir OSySy’, (1.4) 
ty (x, w(x) +2)| <4, (x, z|) (1.5) 
gilt, soweit die vorkommenden Argumentepaare in & bzw. © liegen (vgl. Hilfs- 
satz I, 1). 
b) u(x) sei eine im Gebiet & verlaufende Naherungsfunktion aus Y fiir die 
Lésung w*(x) mit stetigem M[uo]. 
d [uy] = — M[uo] + f(%, Mo) 


ist der zu u(x) gehorige Defekt der Differentialgleichung. 


09(x) bedeute eine stetige Funktion mit 


| % (x) — w(x)| S ao(*) S O(a). (1.6) 


(Ihre Existenz wird vorausgesetzt.) 
c) Im Eindeutigkeitssatz 2 sind ferner gewisse Voraussetzungen erforderlich, 


welche bewirken, daB die Integralgleichung (3.6) der Aufgabe (3.7), (3.8) aqui- 
valent ist: 


Die Aufgabe M[uj=r(x), U,[u] =o (ead ee) 


sei lésbar, wenn 7(x) einer bestimmten (von der speziellen Aufgabenstellung ab- 


hangenden) linearen Menge stetiger Funktionen 8 angehért. 8 enthalte die 
Funktionen? 


2=|d[mol]| ; 

z= 09(%), 

z=h[u] bei beliebigem w € Y, (127) 
B=] (a, Ule))) bei behebigem w€ 8, (1.8) 
z = G(x, &) u(&) dé bei beliebiger stetiger Funktion w(x). (1.9) 


* Enthalt 8 auBerdem die Lésung u=g(x) der Aufgabe M[u]=0, U;[uJ=y,;, 
so folgt aus den Voraussetzungen dieser Nr. 1c, daB die Integralgleichung (I, 2.6) der 


natin 2.1), (I, 2.2) Aquivalent, die in Nr. I, 2.1 genannte Voraussetzung c) also 
eriu 1St. 
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Bei gewohnlichen Differentialgleichungen sind die genannten Voraussetzungen 
im allgemeinen fiir 8 als Menge der stetigen Funktionen erfiillt. Bei Aufgaben 
der in Nr. I, 3.2 behandelten Art verwendet man als 2 die Menge der auf 8 stetigen, 
auf 8 mit stetigen ersten und zweiten partiellen Ableitungen versehenen, langs 
I stetig differenzierbaren Funktionen u(x)=u(x,, x.) und als 8 die Menge der 
auf % stetigen Funktionen, welche auf 8 einer Hélder-Bedingung geniigen. Dann 
sind die Funktionen (1.7) und (4.9) aus 8 [/], und es liegt auch die Funktion (1.8) 
in 8, falls f(x, y) stetige erste Ableitungen besitzt. 


2. Ergebnisse 


Satz 1. Die in Nr. 1, a und b genannten Voraussetzungen seien erfiillt, und es 
gebe eine Funktion n(x) € 2X mit stetigem M{[n] und folgenden weiteren Eigenschaften : 


0Sn(x) SO(%) — @(%), (2.1) 
| 4[mo]] S M[q] — f(x, + 0) + f(%, @), (2.2) 
|2[%] — ¢(x)| SAly], (2.3) 
p(x) +(x) S u(x) Sy(x) — (x). (2.4) 


Dann besitzt die gegebene Aufgabe (0.1), (0.2) (d.h. (I, 2.4), (I, 2.2)) eine Lésung 
u*(x), und es gilt die Fehlerabschatzung 
| w(x) — uo (x)| S (x). (2.5) 

Dariiber hinaus kann man aussagen, daB unter den Voraussetzungen dieses 
Satzes die nach (I, 1.6) gebildete Folge wu, (x) gleichmaBig gegen u*(x) konver- 
giert und die Fehlerabschatzung | u*(x) — uw, (x)|<7(«*) + @9(*) — n(x) (1 =0, 1, 
2,...) mit den durch (I, 1.15) erklarten Funktionen @, (x) gilt. 

Die Voraussetzung (2.3) ist z.B. erfiillt, falls u(x) den gegebenen Rand- 
bedingungen (0.2) und 7(x) den zugehérigen homogenen Randbedingungen ge- 
niigt, denn dann ist [u)|=g(x) und h[{y]=0. Bei partiellen Differentialglei- 
chungen und komplizierterem Grundbereich §% ist es oft praktisch nicht méglich, 
U(x) und (x) so zu wahlen. In vielen Fallen ist (2.3) aber auch erfiillt, wenn 
Uo(x) und 7(x) auf dem Rand gewissen Ungleichungen geniigen: 


Zusatz 1. Aus 


Vip eaee Ban Zz) bel eeUiese Nummern 1 =1,,1,-.- (2.6) 
=0 fiir 1 > 1; 
(u(x) EU) folge w(x)20. Dann kann man (2.3) durch 
[U.[%] — yi | = Gln] fiir 1 = 11, %5,..- (2.7) 
Ulmol=ri, Uly]=0 fir i+4; 


erseizen®. 

Nimmt /,(x, y) in der Umgebung der Lésungsfunktion auch negative Werte 
an, so wirkt sich bei der Fehlerabschatzung oft ungiinstig aus, daB in (1.5) der 
Betrag |/,(x, w+2)| abzuschatzen ist. Man kommt jedoch sehr oft auch mit 
einseitigen Abschatzungen aus: 


8 Bei der in Nr. I, 3.2 behandelten Aufgabenklasse geniigt es also, auf I” |u)—y| <7 
zu fordern. 
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Zusatz 2. Fiiy jede Konstante C=O erfiille N[w] =M[u|+Cwu die in Nr. la 
und die in Zusatz 1 von M{u] geforderten Voraussetzungen. Dann bleibt der Satz 
richtig, wenn man Statt (1.5) nur 

Ty (hs Oe) 2) i, eaaral) (x,w+zEG; x, |2/€ S) (2.8) 
verlangt und (2.3) durch (2.7) ersetzt?. 

Zusatz 3°. Statt (2.4) braucht man nur zu fordern, dap 

Max [®(x), (x) —n(@]Z es), Min[P(x), (x) +n] <¥(%) (2.9) 


gilt, wenn D(x) und P(x) irgendwelche Funktionen bedeuten, mut denen 


Plasreleh se | Gleb lb yee se 14) (2.10) 


fir x, yvE® st. 
Satz 2. Die Voraussetzungen des Satzes1 (bzw. die abgednderten Voraus- 
setzungen der Zusdtze 1 und 3) und die der Nr. Ic seien erfiillt. Besitzt dann die 


Aufgabe M[C] =| [moll 4 f(x, L 09) f(x, Dale 


U[e| = Ul) @=1,2,.--) 
nur eine Lésung €(x) mit OSC (x) <n (x), so ist u*(x) auch die einzige 1m Gebiet © 
verlaufende Lisung des gegebenen Problems (0.1), (0.2), fiir welche (2.5) gilt. 
Satz 2 bleibt auch richtig, falls N[w] =M[w}+Cu fir C&O alle in Satz 2 


und Zusatz 1 von M{[w] geforderten Eigenschaften besitzt und man (1.5) durch 
(2.8) sowie (2.3) durch (2.7) ersetzt. 


3. Beweise 
Beweis des Satzes 1 und des Zusatzes 3. Wir zeigen, daB mit den Funktionen 


ols), a(x) = ela) +f G(x, &) H(E,wol@)) a8 6.1) 
ee a(x) =Heln] +f G(x, §)|deluo]] a8 * (3.2) 


die Voraussetzungen des Satzes I, 1 erfiillt sind. 
Fir jede Funktion wE€% mit stetigem M{[w] gilt 


4 G(x, é) M.[u] dé = u(x) —h[u] ®. (3.3) 
Damit ist wegen (2.3) 


|, (*) — to(x)| =| g(x) — luo] + f G(x, &) de [org] a | 


% 


<h[n] Ta G(x, &) |d, [|| dé = 0, (x). 


* Beider gewohnlichen Randwertaufgabe (M[w] =) —u’”=a (~) + b(*) u(=f (4%, w)) 
u(0) =a, w(1) = kann man damit / = (|6| +) y statt j= |b| y verwenden, und der 
ieee te sin wx fihrt dann auf die notwendige Bedingung — «0 <b <2? statt 

tae 
Verwendet man die Zusdtze 2 und 3 gleichzeitig, so bedeuten hier g(*) und 

(% §) die zu N[w] = M[u] + Cu statt M[u] gehérigen Funktionen, wobei C wie im 
Beweis von Zusatz 2 zu wihlen ist. 
° Der Index & deutet an, da® die Differential 1 i 
an, da operation d (bzw. WM , 
als Funktion von € anzuwenden ist. 2 ae 
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Aus (2.2) und (3.3) folgt weiter 


0; (x) <h[n] +f G(x, €) [Mz [n] — #(E,n(€) + 00(6)) FALE, 09 (€))] dé 


= 1) (x) ah G(x, ) [F(E,m() + @0() —F(E, 00 (8) 48, 
d.h. i 
1 (%) 2 04 (x) zi G (x, &) [#(E,n() + @o(€)) — FE, eo (€))] 4. (3.4) 


o (x) =H (x) —o,(x) geniigt damit der Ungleichung (I, 1.16). 
Man hat nach (3.4) auch o(x) 20 und wegen (2.1) a(x) <0 (x) — o9(x) —o,(x), 
so daB ferner 09(x) +0,(*%) <@(x) ist. SchlieBlich gilt 
ty (x) — 0 (x) = M(x) + (14 (x) — Ug (x)) — 97 (x) + 04 (x) = ug (x) — (2), 
te (2) +(x) = to (x) + (o4 (2) — t4o(x)) g(x) — 04 (x) S g(x) q(x), 


so daB infolge (2.4) auch (I,1.17) bzw. unter der Voraussetzung (2.10) auch 
(I, 1.20) erfiillt ist. 


Beweis des Zusatzes 1. (2.7) bedeutet wegen U,[g]=y, (¢=1, 2,...) 


SAO NAA ci Bea See aes. 


nasal 


O 

2M), Bie ae = a, oe 
U, Uy—egli Ae 
Masts heater a 


Daraus folgt nach Voraussetzung 


h[n—uw+g]20 und h[yn+u—g]=0 


d.h. wegen h[g]=g: 
|} [49] — ¢(x)| =| Fela — g]| SATn. 


Bewets des Zusatzes 2. Es gebe eine Funktion 7 (x)€ YU mit stetigem M[7]|, welche 
den Bedingungen (2.1), (2.2), (2.4) und (2.7) gentigt, wobei fiir f(x, y) nur (2.8) statt 
(1.5) gefordert wird, und es bedeute dann @ (x) =u (x) — (x), p(x) =u (x) +7 (x) 
und @ das in ® enthaltene beschrinkte Gebiet der Punkte x, y mit «C8, 
G(x)Sy<p(x). Wir setzen ferner p(x, y) =Cy+f(x, y) und (x, y)=Cy+ 
f(x, y) mit einer Konstanten C>0, die so groB gewahlt ist, daB in (s py (%, Vy) 20 
gilt. Dann ist |p, (x, w(x) +2)|S,(x, |2|) fiir x, w(x) +2€ @ und x, |2| EW. 

Mit N[u], p(x, y), p(x, y) und & statt M{[u), f(x, ¥), f(a, y) und ®& sind dann 
alle Voraussetzungen des Satzes 1 erfiillt, und es folgt die Existenz einer Losung 
u*(x) mit (2.5). 

Bewets des Satzes 2. Aus (2.2), (2.3) folgt — vgl. (3.4) — t=Tt mit T=7-F Oo 
und dem durch 


T 0 =o, (x) + 09(x) oe G(x, 8 [f(E, 0) —/(E, oo) ] dé 
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definierten Operator T und daraus nach Nr. I, 1.5 weiter die auf 8 gleichmaBige 
Konvergenz der durch @,;;= 70, (7 =0, 1, 2,...) erklarten Folge 0, gegen eine 
stetige Funktion 9* = 7 0* und 
on (2) Sensa(%) Sor) St(x) — (w=0,1,2,...). 6.5) 
Fiir € =0*— Q, gilt also 
C(x) =h[) +S G(x, ) |de[uo]] 2E +S C(x, &) [76,5 +008) FE ol6))] 48. (3-6) 
B B 
Da C(x) eine Summe von Funktionen der Arten (1.7), (4.9) ist, gehort C(x) zu ee 
Daraus folgt nach Nr. 1c weiter, daB die rechte Seite der folgenden Differential- 
gleichung (3.7) eine Funktion aus 8 und €(x) daher Losung der Aufgabe 
M{C] =| 4 [o™a]] + 7(%,5(2) + g0(4)) — F(%, @0(#)), (3-7) 
Ul] =U[n] (4 =1,2,...) (3.8) 
ist. AuBerdem hat man infolge (3.5) 
0S¢0(x) Sy (x). (3.9) 
»i15-2,. v2) istreine. Folge vp -dehiniera 


Durch Uj E =NH+ Oo, Up Tt, (2 a 


fiir welche G80) Steg AeA Te) (3.10) 


gilt. Man beweist dies mit vollstandiger Induktion, indem man t=2TJ rT und die 
Monotonie des Operators T benutzt. Entsprechend wie in Nr. I, 1.5 die gleich- 
maBige Konvergenz der Folge o,,(%) gegen o* = T o* gezeigt wurde, beweist man 
hier, daB die Folge t,,(x) gleichmaBig gegen eine Funktion o** = T o** konver- 
giert, fiir welche wegen (3.10) 0*(*) S0**(x) Sy (x)+ 09(x) gilt. Auch € = 0**— 9, 
geniigt den Beziehungen (3.7), (3.8) und (3.9). Daraus ergibt sich nach Voraus- 
setzung 0*— 0¢7=0**— 0). Es konvergiert also auch die Folge 1, gegen 0*. Mit 
Satz 1b aus [3] folgt dann die Behauptung. 


4. Rechengang und numerisches Beispiel 


Zur Ermittelung einer Naherungslésung w(x) und anschlieBenden Fehler- 
abschaétzung wird (als Beispiel) der im folgenden beschriebene Rechengang vor- 
geschlagen. Die Rechenschritte 2 und 3 kénnen dabei auch durch andere Metho- 
den zur Berechnung von u, ersetzt werden. 


1. Man priife, ob die in Nr. la genannten Voraussetzungen erfiillt sind (f(x, y) 
wird jedoch erst im Rechenschritt 5 ermittelt) und terle die Randoperatoren U, in 
zwer Klassen (1 =1, bzw. 1+-1,) ein, derart daB aus (2.6) u(x) =0 folgt. Nimmt man 
alle U; zur zweiten Klasse (i +:1,), so ergibt (2.6) w(x) =0. Die erste Klasse (i =1;) 
sollte aber im allgemeinen méglichst viele U, umfassen. — Diesen Schritt 1 wird 
man nicht bei jeder einzelnen Aufgabe, sondern fiir gewisse Klassen von Problemen 
durchfiihren (siehe z.B. Nr. I, 3 und I, 4). 


ESTE Cae y) im y nichtlinear, so berechne man mit dem Differenzenverfahren 
im Punkten XpEB Naherungen v(x,) fiir die Funktionswerte u*(x,) der Lésung. 
Bei linearen Aufgaben entfallt dieser Rechenschritt. 
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3. Man mache fiir die Lésungsfunktion u*(x) einen Ndherungsansatz 
U(x) = ¥9(%) + >) a, v, (x) (v, (x) EU mut stetigem M{2,]), 
v=1 


welcher fir alle Werte. der Parameter x, mindestens die Randbedingungen U, [uy] =y; 
mit 1=1, erfillt. Die Parameter «, bestimme man so, daB die in uy linearen Nahe- 
rungsdefekte 


d [9] = — M[uo] + f(x, 0(2)) + fy (x, 2(2)) (uo(2) — 0(a)) 
im den Punkten x, und die Randdefekte 
d;=U,|u|—y; fir 1 =4, 
in passend gewahlien Randpunkten dem Betrage nach méglichst klein werden (soweit 


nicht auch die Randbedingungen mit 7 =7, fiir alle «, erfiillt sind). Dazu benutze 
man etwa die Fehlerquadratmethode oder Tschebyscheffsche Approximation. 


4. Man berechne die Defekte 
4 [uo] = —M[uo} + f(x, %) und d;[u%]=Ul[u)—y, fir 1=7%,. 


5. Man ermittele in einem geeignet gewahlten Gebiet & eine Majorante f(x, y), 
welche die in Nr. 1a geforderten Eigenschaften besitat, und bestimme eine (1.6) 
gentigende Funktion @9(x). Im allgemeinen wird man dabei w(x) =u (x) und 
dann @y (x) =O setzen. Sind die Voraussetzungen des Zusatzes 2 erfiillt, so braucht 
statt (4.5) nur (2.8) gefordert zu werden. (Zur Wahl des Gebietes ® und der 
Majorante f(x, y) vgl. auch Nr. I, 2.3.) 

6. Man berechne eine Lésung n(x) =0 (CA mit stetigem M[n]) der Aufgabe (2.2), 
(2.7) und priife, ob die Voraussetzungen (2.1) wnd (2.4) (bzw. (2.9)) erfiillt sind. 
Ist dies der Fall, so existiert eine Lésung u*(x) der gegebenen Aufgabe (0.1), (0.2), 
und es gilt die Fehlerabschatzung (2.5). 

Numerisches Beispiel. Als einfaches Beispiel rechnen wir die schon in Nr. I, 3.4 
behandelte Aufgabe? 


(M[u] =) —a” =3x%0—u8(=f(x,4)), u(0)= 3, “(1) =14. 


1. Die in Nr. 1a genannten Voraussetzungen sind erfiillt. Man kann hier — 
wie meistens bei gewohnlichen Differentialgleichungen — leicht erreichen, daB 
die auftretenden Funktionen die entsprechenden Randbedingungen erfiillen, kann 
also beide Randbedingungen zur zweiten Klasse nehmen. 


2. Das iibliche Mehrstellen-Differenzen-Verfahren zur Schrittweite h = 4 ergab 
die Werte v(x,) (x,=-h) in der zweiten Spalte der Tabelle 1 (vgl. Nr. I, 3.1). 


3. Bestimmt man in dem die Randbedingungen erfiillenden Ansatz 


6 
(x) =$(1 +) +2(1 —%) Da, a 
y=1 
die Parameter «, so, daB die Summe 
8 ~ 
D (4(%%))? 
k=0 


? Es sei bemerkt, daB sich die in Nr. I, 3.1 hergeleitete Fehlerabschatzung verbessern 
1aBt. Mit o,=0,0022 erhalt man dort als mégliche Schranke o(%) = 0,0055¥(1 — *) 
<0,0014. 
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mit Fi (x,) = ul (x4) — (6x4 0%) — 3.07%) Mo (Hp) + 3% 0? (He) — 20°(%) 
méglichst klein wird (Fehlerquadratmethode), so erhalt man 

ot; =0,134285, 4 = 0,036253 
&)=0,196777,  % = 0,018504 
&, = 0,152856, a = 0,037323 


und die Werte u(x,) in Spalte 3 der Tabelle 1. 


Tabelle 1 
eee ie ee ee ee ee Se 

8x v(x) Uy (x) Uo (x) 

6) 0,5 O55 0,5 

1 0,580145 0,580147 0,580146 
2 0,661 302 0,661 307 0,661 299 
3) 0,741 769 0,741775 OAL TT 
4 0,818 845 0,818 852 0,818 845 
5 0,888 686 0,888 694 0,888 683 
6 0,946 269 0,946 277 0,946271 
7 0,985 564 0,985 569 0,985 563 
8 1,0 1,0 1,0 


4. Spalte 2 der Tabelle 2 enthalt einige Werte des Defektes d[w)|. Die Abbil- 
dung veranschaulicht den Verlauf dieser Funktion. Etwa bei x =0,055 hegt ein 
Maximum, dessen Wert < 0,000 24 ist. 


Tabelle 2 
8x [uo] > 104 ¢, (x) — €2 (x) €3 (*) d [to] - 104 
0 — 0,174 2 0 | 0 — 148,19 
0,5 2,315 2 0,0049 | 0,0002 = 33,52 
1 0,783 2 163 3} 9,09 
1,5 == (O8DE | 2 802>ina| 35 14,89 
2 = 11 DO D 434 | 66 6,74 
2 = 0,530 2 | 539 99 a S| 
3 0,407 | 1,996 604 | 129 = 644 
3,5 0,932 1,946 | 624 149 — 4,86 
4 0,726 | 1,889 598 156 0,21 
4,5 0,036 —| 1,830 oot we 149 5,29 
5 — 0,605 1,774 | 53508 129 7,08 
535 — 0,688 | 12730 321 99 B98 
6 — 0,130 | 1,705 207 66 Sob LIne y ACS) 
6,5 0,510 1,709 108 35 = WL NS 
7 0,281 | 1,753 | 40 13 OG 
ioe — 0,978 | 1,846 0,0006 0,0002 27,95 
8 | = 0,074! | D ) 0 126,90 


. 5. Wir leiten hier eine méglichst genaue Majorante her, wahrend wir in I 
eine besonders einfache verwendet hatten. Die Voraussetzungen des Zusatzes 2 
sind erfiillt, so da statt (1.5) nur (2.8) zu gelten braucht. Fiir 1%, vy) =3 5-2 
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erhalt man bei beliebigem z 


hy (%, Ug + 2) = 3U9(2% — Uy) — 6 (ty — x) z — 3.2? 


O ea VSS 
= vA = ole SS 
<{yasee nag} POMEL 3h (ley sor [PETS 

mit €=$u,(é) ~4 und 
f 0 i GE 
— . | 2 3 
poet eee lenis ror | Fie | 
d [ug] 
f 
0,001. 
10 
a5 r 


Fig. 1. Defekt d[u)] 
Bei w =u, gilt (2.8) mit dieser Funktion f(a, y) also fiir alle Punkte x, w(x) +2 


aus dem Gebiet 
a Cre NOES ors —o<y<o. 


Im Gebiet Re 
OPO S09 OSS O(n) 4G (a) a 


besitzt f(x, y) auBer (1.3) auch die Eigenschaft (1.4). — Man kann f(x, y) auch 
fiir OS y< oo erklaren, muB dann aber fiir y>u)(x) — x (wegen des negativen 
Vorzeichens bei y?) einen anderen analytischen Ausdruck verwenden. 


6. Zur Losung der Vergleichsaufgabe 


[2[“]|S—n"—ft(% 7), (0) =n(1) =0 (4.1) 
machen wir den Ansatz 7 =Ax(1—.x) (4 =const). Die rechte Seite der Un- 
gleichung erhalt dann die Form 


=n!" — F (2,19) = 1%) A + e0(%) A? + g(x) - A 
mit Funktionen c; (x), von denen einige Werte in den Spalten 3 bis 5 der Tabelle 2 
angegeben sind. Die durch Vernachlassigung der in A nichtlinearen Glieder ent- 
stehende Ungleichun 
nie |4[uo]| S (2) 4 


ist fur 0,00024 


A = 0,00012 = °° > Max [4 [eo]] 


C,(*) 
erfiillt. Man stellt fest, daB auch (4.1) mit diesem Wert A gilt; die in A nicht- 
linearen Glieder haben nur wenig Einflu8. 
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Fiir die mit diesem Wert A gebildete Funktion 7 (x) ist auch (2.1) (bei @9=0) 
erfiillt. (2.4) braucht wegen y= — co, y=- oo nicht nachgepriift zu werden. 
Man erhilt also die Fehlerabschatzung 


| u*(x) — u(x)| S74 (x) = 0,00012%(4 — x) S 0,00003. 


Bei der Lésung dieser Aufgabe wurde die am Institut fiir Angewandte Mathe- 
matik der Universitat Hamburg aufgestellte Rechenanlage IBM 650 benutzt: 
4. um das System der nichtlinearen Differenzengleichungen nach einem in [2] 
beschriebenen Iterationsverfahren zu lésen (Rechenschritt 2), 2. um das lineare 
Gleichungssystem fiir die «, aufzustellen und zu lésen (Rechenschritt 3) und 
3. um die Werte der Funktionen u,, d{u | und c; zu berechnen (Rechenschritte 
3, 4 und 6). . 

Wenn man die «, so bestimmt, daB 

7 


D (Ho (%_) — 0(%2))? 

k=1 
méglichst klein wird, so erhalt man die Werte % (x,) in Spalte 4 der Tabelle 4 
und die wesentlich gréBeren Defekte d[u%,)| in Spalte 6 der Tabelle 2. 
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Note on the Expansion in Powers of Poisson's Ratio 
of Solutions in Elastostatics 


ELI STERNBERG & ROKURO MUKI 


Abstract 


In this paper we examine the power-series development with respect to 
Potsson’s ratio of the solution to the second boundary-value problem (surface 
tractions prescribed) in three-dimensional classical elastostatics for the case of 
vanishing body forces. The individual terms of this expansion, by means of 
steady-state thermoelasticity theory, are found to admit an independent physical 
interpretation. This interpretation, in turn, permits certain conclusions concerning 
the dependence upon the elastic constants of solutions to space problems. In 
particular, it is shown that for a body of arbitrary connectivity all stresses are 
independent of the elastic constants if and only if the dilatation (and the mean 
normal stress) is a linear function of the cartesian coordinates. The feasibility 
of successive approximations of solutions appropriate to sufficiently small values 
of Porsson’s ratio is also considered. 


Development of the solution in powers of Poisson’s ratio 


In indicial notation, and with reference to rectangular cartesian coordinates 
x,;, the problem to be considered is governed by the displacement equations of 
equilibrium 

1 
UST aah ie (1) 


together with the stress-displacement relations 


2 
Tig = Uy, 7 4, 5 5 ite, a - (2) 
Here u; and t,; denote the components of displacement and stress, respectively, 
the usual conventions for summation over repeated indices and space-differentia- 
tion being employed; 6;; stands for the Kronecker delta, while ~ and o designate 
the shear modulus and Porsson’s ratio. The field equations (1), (2), which must 
hold throughout the region ZY occupied by the medium, are subject to the 
boundary conditions PR eer “om ee, 3) 
where J, are the direction cosines of the outward normal to the boundary &%, 
whereas ¢; are the components of the given surface tractions and are assumed 
to be functions of position alone. 
AG 
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In view of the manner in which the parameter o enters (1) and (2), the 
solution to the foregoing boundary-value problem depends analytically on o for 
|o|<%. We may therefore expand the corresponding displacements and stresses 
as power series in o in a neighborhood of o=0. Thus 


[e,2) lo, 2) 
u,(x,0) = iw (x)o", — %;; (%, 0) =D th? (x) 0”. (4) 
n=0 n=0 
The argument x here stands for the triplet of coordinates (4%, %2, %3). Equa- 
tions (4) contain no explicit reference to the shear modulus: it is clear from 
dimensional considerations that in the second boundary-value problem the dis- 
placements are inversely proportional to « while the stresses are independent 

of w. Inserting (4) in (1), (2), (3), noting that 


$236 es (20)", (5) 


and making use of the algebra of power series, we obtain 


co n 
> [ods + Ds am ale| o” =o, (6) 
n=0 m=0 
+ (gt OC) > ['s lon 
Dt; oO" = { 2d [uy + w]e" + 6;; D | yo mah] or (7) 
n=0 n=0 n=1lm=0 
CO 
ne EO ae. ON eee (8) 


Since (6), (7), (8) must hold identically in o, we have 


us"). se ul”), se {? (SOO (9) 
a = a [al) + ul”) — 6, f] (» =0,1,2,...), (10) 
ON petty att eetOuw ot 4p Deeaoensl 4G (41)? 
provided , 
n—1 
f(x) =, {” (x) a= 8. > "a (n SSH YS oe “) A (12) 
m=0 
From (9) and (12) follows 
Mi=0  (w=0,1,2,...), (13) 


whence /”) (x) is harmonic in J. 


Interpretation. Successive approximations 


Let S(x, 6) and S™ (x) (n =0, 1, 2, ...) denote the (ordered) array of u; (x, o), 
7; ;(x,0) and u!" (x), x(x), respectively. According to (4) 


S (x) = S(x, 0), (14) 


Se that S (x) is the solution of the original problem for the special case in 
which Potsson’s ratio vanishes, as is also confirmed by comparing (9), (10), (44) 


a} . P . 
The given surface tractions were assumed to be independent of the elastic 
constants. 
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for n =0 with (1), (2), (3) foro =0. With a view toward a physical interpretation 
of S” (x) (n=1) we recall? that (1) and (2), in the presence of a temperature field 
T(x), take on the modified form 


1 _ 2(1+0)« 
Magyriitges se Uy tem ae tes (15) 
20 2(1 
Tip = MMs, 5 +g 0; j he he a 0; ; tL, (16) 


L2G 
if « is the coefficient of thermal expansion. A comparison of (9), (10) for n=1 
with (15), (16) reveals that S” (x) (7 =1, 2,...) are the solutions of a sequence 
of boundary-value problems in steady-state thermoelasticity. Bearing in mind 
(14), (13) we find specifically that u‘”) (x) and ey (i) (=A 2)... arethesthernal 
displacements and stresses produced in an elastic body occupying the region G 
by the harmonic temperature distribution T=/™)(x) in the absence of body 
forces and surface tractions, provided the shear modulus of the material is again y, 
PolIsson’s ratio o =O, and the coefficient of thermal expansion « =. 

Equations (9), (10), (11), which characterize S”(x) (7=0,1,2,...), in 
principle provide a systematic means for the determination of successive approxi- 
mations to the solution S(x,o) of the original problem in the vicinity of o =0. 
Thus, once S(x) has been established, /" (x) is known from (12) and the 
determination of S®(x) becomes possible; a knowledge of S™ (x) and S (x), in 
turn, permits the computation of /(?)(«) etc. As far as the actual process of finding 
S™ (x) for known f(x) is concerned, we note that the general solution of (9) 
admits the representation 

ug”) = G+ x, Yh) — 3p” (17) 


with 
eae eh Giese) iy Orie ae.) (18) 


Equations (17), (18), except for the obvious modification needed to accommodate 
the temperature field T =/”(x), constitute the special form assumed by the 
PAPKOVICH-NEUBER solution® of the displacement-equations of equilibrium if 
o=0. The boundary-value problems characterizing S™)(x) for n=1 may 
alternatively be attacked with the aid of GoopiER’s [3] method of integration 
of the thermoelastic field equations. 

Unfortunately, the successive-approximation scheme just described holds 
little promise of practical significance since the actual determination of each 
S™ (x), in general, is not apt to be appreciably simpler than the task of finding 
S(x, 0) explicitly. The application of the present approach to the three-dimen- 
sional treatment of the plane problem appears particularly inviting at first sight: 
in this instance S)(x) is available beforehand from the associated plane-strain 
solution’ and one may expect S®(x) to supply useful three-dimensional correc- 
tions provided o is sufficiently small. Yet, the particular difficulties peculiar to 
the determination of the three-dimensional solution S (x, o) of the plane problem 
are also encountered in connection with S\” (x) for n=1. 


2 See, for example, MeLran & Parkus [J]. Numbers in brackets refer to the list 
of references at the end of this paper. 

3 See, for example, SOKOLNIKOFF [2], p. 330. 

4 Recall that the plane-strain solution is the exact three-dimensional solution of 
the plane problem for = 0. 
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Influence of the elastic constants on the stresses 


The physical interpretation of the coefficients in the expansion (4) arrived 
at earlier permits certain conclusions regarding the dependence on the elastic 
constants of the solution S(x,o) to the second boundary-value problem for 
vanishing body forces. 

We prove first that all stresses belonging to S(x,o) are independent of the 
elastic constants if and only if the dilatation ™, , (and thus® the first invariant 
of stress 7;,) is a linear function of the cartesian coordinates (%,, %2, 3). If the 
region Zis simply connected, this theorem follows at once from the characterization 
of the stress field 1,;(%,o) through the stress equations of equilibrium and 
compatibility 

T7540; Cy err SF (19) 


in conjunction with the boundary conditions (3). Moreover, the condition 
The, 17 =O Clearly remains necessary for t;; to be independent of o even if Zs 
multiply connected. Here, however, the suffictency of this condition cannot be 
inferred from (49), (3) which then no longer determine the stresses completely 
but must be supplemented by the requirement of single-valued displacements. 
To show that t, ;;=0, 7.€. 4% ,;;—=0, assures the independence of the stresses 
of the elastic constants regardless of the connectivity of YZ, we note from (4) 
that “, ,,;=0 implies uf; =0 (n =0, 1, 2,...). Hence, because of (12), 7%) =0, 
so that the temperature fields /”) (x) are linear functions of (%,, ¥2, %3) for n=1. 
Now, according to a well known theorem in steady-state thermo-elasticity 
theory, the thermal stresses produced by a temperature distribution 7 (x), in the 
absence of loads, vanish identically if (and only if) T(x) is linear®. Consequently, 
and by virtue of (11), 7!" (x) =0 for m=1. The second of (4) therefore yields 


that 1; ;(%, 0) 79 (x), whence t;; is independent of a. 


vq 

In view of the second of (19), t;,,;;=0 is equivalent to t;; ,,—0, and we 
may conclude that the stress field of S(x, o) does not involve the elastic constants 
if and only if all cartesian components of stress are harmonic (rather than merely 
biharmonic) functions of position. In these circumstances, the maximum 
principle for harmonic functions predicts that each stress component 71,; attains 
its maximum and minimum on the boundary. According to a result of Pérya [4], 
the harmonicity of all t;; guarantees further that the maximum and minimum 
values of the principal stresses, as well as the maximum of the total strain energy 
density, along with the maxima of the energy densities of volume change and 
of distortion, all occur on the boundary’. The problem of SAINT-VENANT torsion 


probably affords the most significant illustration of the type of degeneracy under 
discussion. 


° Recall that t,,= [2u(1+6)/(4 — 20)]Uz ,, aS is apparent from (2). 

* See [1], p. 9. The proof given in [1] applies only to simply connected regions. 
T o complete the argument for multiply connected regions it suffices to observe that 
if T(*%) =cy+c;%;, with co, c; constant, the single-valued displacement field u,;= 
0 (Co %; — 30;%;%;+C;%;4%;) generates vanishing stresses in the sense of (16). 

’ Thus, in particular, if the stresses do not involve the elastic constants, plastic 


deformation must initiate at a point of the boundary provided the material obeys 
the Hencxy-MiskEs yield condition. 
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It is evident from the preceding remarks that the complete absence of the 
elastic constants from the stress distribution is characteristic of an exceedingly 
narrow Class of three-dimensional problems in elastostatics. This state of affairs 
is in marked contrast to the two-dimensional theory of generalized plane stress 
where all stresses are independent of the elastic constants for a wide variety of 
physically important problems®. We turn, finally, to examples of space problems 
in which the stresses may be predicted to be partly independent of the elastic 
constants with the aid of the present analysis. Consider an elastic body occupying 
the half-space *,=0 or an infinite elastic slab bounded by the planes x,= +h. 
In either case let an arbitrary finite portion of the boundary be subjected to 
given normal and shearing tractions, while all stresses are required to vanish at 
infinity®. Here t,,(%,o), and thus ™, ,(x,o), must tend to zero at infinity; 
consequently, by (4), the same is true of u\", (n =0, 1, 2,...) and hence also of 
the temperature fields / (x) (7 =1, 2,...) defined by (12). On the other hand, 
it was shown in [7] and [8] that a harmonic temperature distribution, vanishing 
at infinity, applied to either of the two regions under consideration, gives rise 
to a thermal-stress field that is plane and parallel to the boundary!®. Hence 
a\") (x) =O (n =1, 2, ...), and (4) insures that the transverse stresses %,, T12, T13 
are independent of the elastic constants in both problems at hand. This prediction 
is in agreement with the known general solution for the semi-infinite body under 
distributed surface loads"; it is also consistent with special results for the problem 
of the slab due to SNEDDON [//] and Mux! [12]. 

These limited conclusions concerning the influence of the elastic constants 
on the stresses in elastic solids are perhaps of some interest in three-dimensional 
experimental stress analysis, in connection with questions attached to the transi- 
tion from model to prototype. Considerations of this kind are also relevant to 
the stress analysis of viscoelastic materials in view of the correspondence prin- 
ciple? that links the linear theories of elasticity and viscoelasticity. 


The results communicated in this paper were obtained in the course of an investiga- 
tion conducted under Contract Nonr 562(25) of Brown University with the Office of 
Naval Research, Washington, D.C. 
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Berechnung der Druckverteilung an Gitterprofilen 
in ebener Potentialstromung 
mit einer Fredholmschen Integralgleichung 


ERICH MARTENSEN 


Vorgelegt von H. GORTLER 


§ 1. Einleitung 


Zur Berechnung der Geschwindigkeits- und Druckverteilung an endlich oder 
abzahlbar unendlich vielen Hindernissen in einer zweidimensionalen Potential- 
stro6mung gibt es eine Reihe von Arbeiten, in denen das Problem durch Anordnen 
von Singularitaten im Inneren der Hindernisse behandelt wird. Bei dieser Methode 
besteht jedoch die wesentliche Schwierigkeit, auf die wohl zuerst von KARMAN [7] 
aufmerksam gemacht hat, darin, da sich die wirkliche Potentialstrémung des 
AuBengebietes im allgemeinen nicht bis an den a priori festgesetzten Ort der 
Singularitaten im Inneren analytisch fortsetzen la4Bt. Mit wachsender Zahl der 
Freiheitsgrade einer solchen Singularitatenanordnung oder -belegung konvergiert 
das zugehorige Str6mungsfeld daher im allgemeinen auch nicht gegen die gesuchte 
Lésung des Problems, selbst wenn es mit Hilfe der zur Verfiigung stehenden Frei- 
heitsgrade gelingt, Stromlinien durch eine beliebige Anzahl von Punkten der 
gegebenen Berandungen hindurchzufiihren. Insbesondere bekommt man _ bei 
derartigen Verfahren keinen Anhalt fiir den Fehler, mit dem eine bestimmte 
Naherung behaftet ist. 

Die beschriebene Schwierigkeit hat wahrscheinlich erstmals PRAGER [2] im 
Falle eines umstrémten Einzelprofils dadurch vermeiden kénnen, daB er aus- 
schlieBlich die Profilkontur mit Singularitaten, und zwar nur mit Wirbeln belegt, 
da Quellen auf der Berandung die Umstrémungsbedingung verletzen wiirden. 
Als weiterer Vorteil dieser Methode erweist sich, daB diejenige Wirbelbelegung, 
die die Randbedingung erfiillt, gleichzeitig mit der gesuchten Str6mungsgeschwin- 
digkeit am Profil identisch ist!. Auf diese Weise fiihrt PRAGER das Problem auf 
die Auflésung einer Fredholmschen Integralgleichung zweiter Art fiir die (exakte) 
Geschwindigkeitsverteilung am Profil zuriick, die ihrerseits mit der gesuchten 
Druckverteilung in bekannter Weise durch die Bernoullische Gleichung ver- 
kniipft ist. 

1 Diesen Zusammenhang hat schon vorher PRANDTL [3] unter etwas allgemeineren 
Voraussetzungen hergeleitet und folgendermafen beschrieben: ,,Man denke sich das 
Innere des K6rpers durch ruhende Fliissigkeit von dem Druck p)+ 4, wie er im Stau- 
punkt herrscht, ersetzt; dadurch tritt an die Stelle der K6rperoberflache eine Wirbel- 
schicht mit einem Geschwindigkeitssprung von dem Betrage v. Diese Wirbelschicht 
ist das gesuchte gebundene Wirbelsystem.“ 
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Die Methode der ausschlieBlichen Wirbelbelegung des Randes ist auch zur 
Behandlung eines Gitterprofils (Abb. 1) in zweidimensionaler Potentialstr6mung 
herangezogen worden. GOLDSTEIN und JERISON [4] entwickelten auf dieser 
Grundlage ein Verfahren zur Bestimmung einer Gitterkontur mit vorgegebener 
Geschwindigkeitsverteilung, wahrend Isay [5] bei gegebenem Profil eine Integral- 
gleichung erster Art fiir die Wirbelbelegung aufgestellt hat. Isay kann die 
Existenz von Lésungen dieser Integralgleichung und die Konvergenz eines Itera- 
tionsverfahrens zur praktischen Bestimmung der Lésungen jedoch nur beweisen, 
wenn eine bestimmte Ungleichung, in die nur die Gittergeometrie eingeht, erfiillt 
ist. Es liegt daher nahe, in Anlehnung an die Pragersche Arbeit [2] eine Integral- 
gleichung zweiter Art fiir Gitterprofile herzuleiten und zu untersuchen. Dies 
soll in der vorliegenden Arbeit geschehen. Es zeigt sich, daB die Frage der Existenz 
von Lésungen dann ganz allgemein durch die Fredholmsche Theorie [6], [7], [8] 
beantwortet werden kann. Die Identitat von Wirbelbelegung und Geschwindig- 
keitsverteilung wird dabei nicht vorausgesetzt und itiberdies der Begriff der 
Wirbelbelegung ganz vermieden. 

Zur numerischen Auflésung der erhaltenen Integralgleichung kénnte man eines 
der bekannten und sehr genauen Verfahren von NystrOm [9], [10], [77] heran- 
ziehen, die z.B. auf der GauBschen oder Tschebyscheffschen Quadraturformel 
beruhen. Es soll jedoch fiir ein derartiges Verfahren die bequemere Rechteck- 
formel verwendet werden, da diese im Falle periodischer Funktionen die gleiche 
Genauigkeit wie die GauBsche Formel besitzt [72]. An Hand von Beispielen wird 
gezeigt, daB das Berechnungsverfahren numerische Naherungen liefert, die mit 
kleiner werdender Schrittweite sehr schnell gegen die strenge Lésung konver- 
gieren. 

§ 2. Formulierung des Problems und unmittelbare Folgerungen 

Das gegebene Gitterprofil (Fig. 1) bestehe aus kongruenten, einfach zusam- 
menhangenden Innengebieten By), 8,1, 8_,... mit doppelpunktfreien, stetig 
gekriimmten Berandungen ©), ,,,€_,..., die mit gleichem Abstand¢> 0 (,,Gitter- 
teilung’’) in y-Richtung angeordnet sind. Die ©), €,,,€_,... sollen sich gegen- 
seitig weder beriihren noch iiberschneiden. Die Abszissen der Berandungen 
Gy, ©,,, ©... gehdren simtlich dem Intervall a. <x<a, an. Die Menge aller 
Punkte der (%, y)-Ebene, die weder zu 8, 6), B=... noch zw @,;€.,;C.,... 
gehoren, bezeichnen wir als AuBengebiet &. Es gibt immer eine Parameter- 
darstellung x(y), y(p) fiir ©) mit folgenden Eigenschaften: 


x(p) und y(y) sind zweimal stetig differenzierbar fiir alle p; 
x(p) und y(q@) sind periodisch in 27; 


«(p) und y(@) durchlaufen die Berandung ©, mit wachsendem PoSPLM+22 
genau einmal im Gegenuhrzeigersinn; 


ftir alle p verschwinden die Ableitungen %(p) und ¥() nicht gleichzeitig: 
(1) [x(e) ? + Ly) P >. 


Unter diesen Voraussetzungen sind alle ,,Stromfunktionen“ W(x, y), die den 
folgenden Bedingungen geniigen, gesucht: 


V(x, y) ist eindeutig und zweimal stetig differenzierbar in Y; 
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Y (x, y) ist eine in 2% harmonische Funktion: 


(2) ALY: 
in WU besitzt W(x, y) die Periodizitatseigenschaften 
(3) L(nyt)=2(%9), Blny +t) =¥,(2,9); 


OW/dn? ist stetig auf G); 
G ist Stromlinie, d.h. auf Cy gilt 


(4) CW) "Const - 
fiir ein Intervall y — > Sysyt+ + * existieren die beiden Grenzwerte 
20H 
(5) dim_[e* FW (x,n)| =08 


im Sinne gleichmaBiger Konvergenz. 
Der Zusammenhang mit der Strémungsgeschwindigkeit (x, y) =(v,(x, y), 
vy (x, y)) ist in W und auf €),€,,,€_,... gegeben durch! 


(6) ==—TOt( Pt) = erad yt |: 
in Komponenten zerlegt : 
(7) v=, ( 


Mit (2) folgt div» =0 und rot » =0 in Y. 

Sei ferner t der in Gegenuhrzeigersinn gerichtete Tangentialeinheitsvektor von 
&,, so gibt, da © Stromlinie ist, 
(8) v = (, t) 
die Geschwindigkeitsverteilung auf ©, an und ist positiv oder negativ, je nachdem 
die Umstré6mung an der betreffenden Stelle im Gegenuhrzeiger- oder Uhrzeiger- 
sinn erfolgt®. Man bekommt aus (6) und (8) 


p—=i(lerad tit) = (erad- alt thy; 
da [{f, t] nach innen gerichteter Normaleneinheitsvektor ist, folgt 
OV 
2 ae 


auf ©); v ist auf ©) stetig. In iiblicher Weise beschreibt 


(10) (he Sitti ax © 


die Zirkulation um ©). 


1 9/dm soll stets die Differentiation in Richtung der a4uBeren Normalen einer ge- 


schlossenen Kurve bedeuten. 

2 Es geniigt, fiir y irgendeinen bestimmten Wert anzunehmen, da dann aus (3) 
sofort die Giiltigkeit von (5) fiir alle y folgt. 

3 Durch diese Bedingung werden z.B. solche (harmonischen und _ periodischen) 


Funktionen ausgeschlossen, die sich im Unendlichen wie Goj ans sin 282 verhalten. 


47,;,£ sind im folgenden die Einheitsvektoren in x-, y- und z-Richtung. 
5 Die Definition des Vorzeichens befindet sich hier in Ubereinstimmung mit 
PRAGER [2]; sie ist in der Literatur nicht einheitlich. 
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Wir stellen im folgenden einige Eigenschaften der zugrunde liegenden Funk- 
tionenklasse ¥/(x, y) zusammen. 7 
Satz 1. Der iiber irgendeinem Periodenstreifen genommene Mittelwert von » ist 
rechts bzw. links vom Gitter jeweils gleich einer vektoriellen Konstanten 04 .0= 
(v Vs 00.) B20. D_go=(V_0,41 V—oo,y) 4-H. fiir alle Z>0 und alle y gilt 
$00, #2 “+ 00,¥ = 5 


(11) 


B 


(21,41) 


ines Mea AB 


8y 


Fig. 1. Gitterprofil mit Integrationswegen € und (hy 


Beweis. Wir betrachten zunachst ein rechts vom Gitter gelegenes einfach 
zusammenhangendes Gebiet 8 mit doppelpunktfreier, stiickweise glatter Beran- 
dung (©; oberer und unterer Teil von € sollen kongruent sein und in entsprechenden 
Punkten voneinander den Abstand ¢ haben (Fig. 1). Dann ergibt sich bei Beach- 


tung der Periodizitatseigenschaften (3) die x-Komponente in (11) mit Hilfe von 
Fig. 1 aus der Eindeutigkeit von Y(x, y): 


t t 
M+ 9 Yor 5 
1 1 : 1 ov ; 
pf rele an—4 f veleon)an=4 f as=o1; 
ee Vo ; € 


. ales soll stets die in Gegenuhrzeigersinn gerichtete tangentiale Differentiation 
auf einer geschlossenen Kurve sein. 
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die y-Komponente in (41) ist eine Folge des Gau8schen Integralsatzes: 
t 


t 
nt > Yor 5 
1 1 1 ov 1 
mG [ % (0.7) dn = [FS ds=—— [ [4% axdy =o. 
Di : Ce rs € B 


Entsprechend gestaltet sich der Beweis fiir die Verhdltnisse links vom Gitter. 


Satz 2. Zwischen den zu beiden Seiten des Gitters genommenen Mittelwerten 
von » besteht die ,,Umlenkformel 


(12) Di-co Doo = ie 


Beweis. Das zweifach zusammenhangende Gebiet Sj werde von Cf und ©, 
berandet; Gf habe dieselben Eigenschaften wie ©; auBerdem sollen Cf und 6, 
keinen gemeinsamen Punkt besitzen (Fig. 1). Sodann findet man mit Hilfe von 
(10) und (11) 


nt5 Yo+ ; se 
E 1 
V40c0,4 — Vico x = = | V_(%,7) dn — [ v.(o,7) dy — 7 [= ds 0, 
V5 Uh cx 
nt > wt 5 
1 1 
U4-00;9 — © —oo,9 Vy{%1,7) dy — oh ‘ Vy (%,4) an 
t t 
ee 2 Yo= ry 
1 OV 1 r 1 Tae) If 
Pee z . ee % Ce : 
t J on 2 t ap APdxdy t f on t 
Cy B* C, 


Es ist im Hinblick auf das folgende niitzlich, einen konstanten Vektor v= 
(Peo, eoi0co,9) AurCh 


(13) Does (Port toc) 
zu definieren; dann gilt mit (12) und (43) 

| Ber 
(14) LER ta ahee, 


Weiter findet man aus (7) und den x-Komponenten von (11) und (14) zunachst 
fiir alle Z>0 und alle y 


4 |v (es 42,9 + S)—¥ (0, £Z,9-f| = 00,5 


2 


integriert man nun 0Y/@s tiber irgendein Teilstiick der Berandung Gj (Fig. 1) 
mit den Endpunkten (x, y) und (x, y+¢), so erhalt man ganz allgemein in Y% 


(15) (x, 9 +t) = P(x, 9) +t V0, 2° 


Satz 3. Der tiber irgendeinem Periodenstreifen genommene Muttelwert von 
Wx —v,, ,y ist rechts bzw. links vom Gitter jeweils gleich einer Konstanten 
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Choo D2W. Coo, ah. fiir alle Z>0 und alle y gilt 
yl 

(16) f (Plan 2.0) — (a, $2) Klas £2, 1) — Mee, 1] 49 = Cs 

le 

oder hiermit gleichbedeutend 
feu 

(17) + f [ey £Z.0) — (ag £2) Yas + Z,7)] A — 20,2) = Ceo 
mit 


y— 


2 
Beweis. Wir betrachten das Gebiet 8 mit der Berandung © in Fig. 1 und 
bekommen auf Grund des Greenschen Satzes und Gleichung (15) 


Viste ; Yo . 
tes 
Ff Pea Vola 1) dt, 28 | LY (%0, 0) 20% (0, 1)] 4 +020, Yo 
t t 
Lae ont A 2 


=+/ (3 x as al [Y(x, vy) —P(x, y—2)] Fr ds+eo,2 | ds 


€ Coben Goren 
1 ox oy = 
=+ ff (WAx — x AY) dndy— Vo. f & eee O, 
¥ Coven 
Te nes Ox ov ss, : ee 
da Ax und AY in $ und ieee — *] auf © pen, dem oberen Teil von ©, ver- 


schwindet. Entsprechend verlauft der Beweis fiir den Mittelwert links vom 
Gitter 
Im Anschlu8 an (16) bzw. (17) definieren wir noch 


(18) Cap Bl Orcn cat Oe 
Wir bezeichnen c,, und »v,, als die ,,charakteristischen Mittelwerte‘‘ von ¥. 
Lemma 1. Fir jedes Intervall y — ca SySyt+ ae existieren die beiden Grenzwerte 


— 20% 


(49) lim. je" “Tia (= 0 


x—> +00 


im Sinne gleichmadfBiger Konvergenz. 


Bewers. Fiir alle x aus a, <x <-+ oo und alle 7 aus y — BA S95 wicca sind 
W(x, y) und ae : 


é + F(x, 9) =h(x,n) 


voraussetzungsgemaB stetig. Da (5) im Sinne gleichmaBiger Konvergenz gilt, 
existiert zu vorgegebenem ¢>0 ein a’, >a, derart, daB 


[A(em)| <2! 
fiir alle xa’, und alle y — a. = y+s giiltig ist. Ferner gibt es eine Schranke 
M>O0, mit der im Intervall y — a <ynSyt+ a die Ungleichung 


|P(a',.)| <M 
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erfiillt ist. Fiir alle j 
head eo Max {as goin an} 
wt 


und alle y — + <n y+t ergibt sich dann 


22% 2H 2 “20% 
=<. ito eet fe? ud 
ay 
2m a+ 2% Qux 2a 
e= Cinite Morea t (e# —e ?t 


In gleicher Weise findet man den Grenzwert fiir x > — ov. 


§ 3. Zusammenhang zwischen der Stromfunktion im AuBengebiet 
und der Geschwindigkeit am Profil 
Satz 4. Durch thre charakteristischen Mittelwerte C.,, V9. und thre Geschwindig- 
kettsverteilung v auf Gy ist die Stromfunktion ¥ in YX vollstandig bestimmt durch 


(20) V(x, y) = Coo + Yoo, x ¥ — Yoo, y % a 3s [Feds e 
G 
wober 
(21) OE ee ee ee 


\/ [sof oe — cos AES | 


und (€, ) laufender Punkt auf Gp ist. 


Beweis. Wir schlieBen zunachst aus, daB der Aufpunkt (x, y) im abgeschlos- 
senen schraffierten Gebiet in Fig. 2 liegt. Sodann betrachten wir ein dreifach 
zusammenhangendes Gebiet ®, das von einer doppelpunktfreien, stiickweise 
glatten Kurve ©, einer den Punkt (x, y) einschlieBenden Kurve & und von ©, 
berandet wird, wobei keine dieser Berandungen mit einer anderen einen gemein- 
samen Punkt besitzt (Fig. 2). Ferner seien der obere und untere Teil von © 
kongruent mit dem Abstande ¢, wahrend © rechts und links bzw. aus den 


Geradenstiicken 
é = a. , a 


bestehen soll; wir bezeichnen mit (&, 7) stets einen laufenden Punkt in © und auf 
den drei Berandungen von . Fiir & gelte die Parameterdarstellung 


1% O COS T 


(22) §&(0;7) =x + 2% rea n (a; T) =y += are sin 


oO SINT 


1 Diese Beziehung wird von GoLpsTEIN und Jerison [4] auf anderem Wege und 
unter weniger allgemeinen Voraussetzungen, die u.a. das Verhalten im Unendlichen 


betreffen, hergeleitet. 
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bei konstantem, hinreichend kleinem o>0; mit wachsendem 0StS27 wird & 
genau einmal im Gegenuhrzeigersinn durchlaufen. Die Funktion (21) hat auf & 


den konstanten Wert 
t 


(23) F(x, y;&(o; 2), (0; 7);#) =n aE 
wihrend man fiir die Normalableitung auf 

OF m2 G2 ry 
(24) = - \ 4 -- = cos 2T 


’ 

ES, 
YU Uy 
“hy, Mf fi 


Fig. 2. Gitterprofil mit Integrationswegen SG und 


bekommt. SchlieBlich hat das Linienelement von §& fiir wachsende t die Gestalt 


2 52 
1+ a CoS 27 
(25) aS =v0 : —— ~ dt. 
4 wt 


= o2 = 2 : 
(1 = COSe r| (1 gree sin? r| 


2 

t #2 

Indem wir nun beachten, da8 fiir alle Punkte (€, 7) in G und auf G, &, Gp 
Cojf oT ee 


at ee COs 7 = 0 


ist, konnen wir auf Y(é,7) und F (x,y; &,; 2) den Greenschen Satz anwenden 
und erhalten 


[ [(P4F —F AY) ae an 


(26) © , 
© af ees flee eae fue eas 
S K C 
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da AF und AY in & verschwinden. Bevor wir jetzt die drei Randintegrale iiber 
©, K, Cy einzeln untersuchen, vermerken wir noch 


. 27(*¥—é) 
Ql a aN i 
eS Coj = cos oe 
(27) 
a eam) 
Ct ere t 
OMT aes eee 


Zunachst stellt sich dann heraus, daB die Integration iiber den oberen und unteren 
Teil von © zusammen keinen Beitrag liefert; wegen (3), (15), (21) und (27) 
bekommt man namlich 


AB aes oF aw Be oe 
{ey F | as ! [ (PRE asa trme f Fas: 


Soben Sunter 


Soben 
da f eine harmonische Funktion ist, deren sémtliche Singularitaten an den 
Stellen (x, y), (%, y+1), (¥%, y22) ... liegen (Fig. 2), kann der Integrationsweg 
des rechts stehenden Integrals in eine gerade Verbindungslinie der Endpunkte 
abgewandelt werden, und man bekommt mit (27) 
OF 7 Olf 
fi on hs | on 
Soben a 


Mit (7), (11) und (17) ergibt sich weiter fiir den rechten Teil von & 


de — 0% 
n=yt 5 


F(x, yia'.. it) + 2] — 


—|F(x, 95a.) —F (xa) 


¥,(a', .n)} dn = 


(28) ve yey 
, , , x , 
— 7% +f [Ma',.n) —a!, ¥e(a', nag += [ ¥(a',.n) an 
ys y- 5 
y+ 


et (dae a) On. ot ieh 


ao TC a +S Veo, aU V4.0, y *) ? 
Arch. Rational Mech. Anal., Vol. 3 
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wobei wegen (21) und (27) 


22 (*—a',) 
oo 22(y—n) 
; 2nx é COS = i 
Plx,ysa', mit) = 2% Li ee 
a == Tt SN § 
1—2¢ +f cos == Lae : 
: an ety { 
NCES SA Rea: ae iat) oy), sant) 
a a a(yv—yn Sa 
1=2e- * cos = +e * 
+. 
bedeuten. Es existieren die Grenzwerte 
22 
wees 2a(yv—1 

lim P(x, y;a',,n3t) = —77e t cos ay 

a‘, —> +00 vA if 
2u%x 

, r 2nKx 2 rey 

lm... (4, 952, 43 tye * tose: 2a 
a, —> +00 t 


so daB (28) mit Hilfe der Grenzwerte (5) und (19) gleichfalls den Grenziibergang 
a,—>-+ co auszufihren gestattet?: 


: . nes 
lim | (v = -F Jas= = Gude todnY itney Be 
Heo Hee on on 

rechts 


Entsprechend bestimmt sich 


lim (Ye FF) ds = — (0c +202) Y20,y7) 
6 on on : : 
Slinks 


Durch Addition der letzten beiden Gleichungen folgt bei Benutzung von (13) 
und (48) 


(29) Ww te [(¥ om cee <)as mh J 27 (Coo ae Voo, x rR Uo0, y %). 
oes Doors 


Wir wenden uns nun in (26) dem Integral tiber ® zu und erhalten wegen (22) 


bis (25) 


; On on 
K 
22 mo 
(1 + = 00s 21} W(&(o; t), n(o; T)) j aw 
a ——— : -~@t — —— ds; 
j aae Figa) 220 on 
< | (1 at sana ee r| [1 mellem Mees | 
blot t t “ 


da der zweite Summand rechts auf Grund des GauBschen Integralsatzes ver- 
schwindet, kann man den Grenziibergang o->0 ausfithren: 


2% 
. oF OV P 
(30) lim ( in ag) a8 = — | P(E(0; 2), 90; 2)) dv = — 20 P(x, 9). 
K 0 


* Die Vertauschung von Limes und Integration ist hier und im folgenden auf Grund 


gleichmaBiger Konvergenz des Integranden erlaubt, ohne da in jedem Fall besonders 
darauf hingewiesen wird. 
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Fir das in (26) noch verbleibende Integral tiber © bekommen wir auf Grund 
von (4), (9) und wiederum mit Hilfe des Gaufschen Integralsatzes 


OF ae aF e 
(61) he oi F an) OS const [ ds + [Fods = | Fods. 
©, C, ©, ©, 


Durch Einsetzen von (29), (30) und (31) in Gleichung (26), in der die (nunmehr 
erlaubten) Grenziibergange a! > -+ co, a’ — «, o > 0 ausgefiihrt zu denken sind, 
folgt die Behauptung (20) fiir alle in 2 gelegenen Punkte (x, y), die nicht dem 
abgeschlossenen schraffierten Gebiet in Fig. 2 angehéren. Jedoch kann man sich 
die Giiltigkeit auch hierfiir folgendermaBen leicht klar machen. Zunachst beachten 


wir, daB mit (3), (4) und (9) und aus Stetigkeitsgriinden ee =0 und a == 
auch auf ©,,, ©_,... gilt; damit ist ferner ¥konstant und 0W/dn stetig auf €,,, 
@_,.... Liegt nun z.B. der Punkt (x, y) im abgeschlossenen schraffierten Gebiet 


zwischen ©, und ©,, und denken wir uns jetzt das Profil €,, und ® ganz im 
Inneren von G, so fiihrt dieselbe Uberlegung wie oben auf 


(32) PUY) = Coo + U,V — Yoo,y% + see f Fods. 
( 


ast 


Tatsachlich gibt (20) fiir einen solchen Punkt (x, y) denselben Wert wie (32) an, 
da (21) in entsprechenden laufenden Punkten (§, 7) auf © und ©,, gleich ist. 


Satz 5. Die Sivémungsgeschwindigkeit konvergiert 1m Unendlichen gleichmapig 
gegen thre (konstanten) Mittelwerte im Endlichen, d.h. fiir jedes abgeschlossene 
y-Intervall existieren die beiden Grenzwerte 


(33) lim v(x, ¥) = Ukoo 


4% —> co 
am Sinne gleichmabiger Konvergenz. 


Beweis. Auf Grund von (7), (20) und (21) gilt in YU ganz allgemein 


1 oie 2 OF . 
— | 
(64) b@,N=Po tse f (Fi — Fe ijeas 
So 
mit 
Sin nu (%—) 
CL 0 Se t 
g PES : i, 
Ox t Coj eS cos ee 7) 
(35) an 
OF Tt t 
oy t Cof Brie Fay vee 
Es existieren die Grenzwerte 
. OF Se 
im= —— = + —, 
x—> 400 OX t 
? oF 
in = 
x—> +00 0Y 


17% 


246 ErtcH MARTENSEN: 


im Sinne gleichmaBiger Konvergenz, so daB (34) gleichfalls die Grenziibergange 
x—>-+ co gestattet: 


; j = De 
lim v(x, ¥) = Dap le ay feds = 45513 


x—> +00 
C, 


in Verbindung mit (14) folgt dann die Behauptung (33). 


§ 4. Beziehungen fiir den Durchgang des Aufpunktes 
durch den Integrationsweg 


Lemma 2. Sei 0 irgendeine auf ©, stetige Belegungsfunktion,; dann gelten fiir 
die Normalableitungen der in X und By eindeutig erklarten harmonischen Funktion 


a cgetl 
(36) Ulxs=s5 | Feds 
C, 
beim Durchgang des Aufpunktes durch ©, die Sprungrelationen* 


oU 
on 


oU 
fe on 


aes 1 AP 6U 
can De on 


oU 
C, on 


= eg 
i a? 


(37) 


wihrend die Funktion selbst stetig bleibt: 
(38) U|, — U|,=0. 


Beweis. Wir schlieBen an das beim Durchgang durch ©) bekannte Verhalten 
des in X& und B, eindeutigen Potentials 


(39) U(x =5 | (Inge ,)eds 


0 


einer einfachen Belegung @ auf ©) an (7 ist der Abstand zwischen Auf- und Inte- 
grationspunkt), wobei (37) und (38) sinngemaB gelten [13]. Es geniigt dann zu 
zeigen, daB sich die Funktion 


1 t? 
Usa) = Ul.) =e | (Bee Bods 
© 


mit ihren Ableitungen beim Durchgang durch 6, stetig verhalt. Da ohne weiteres 
einzusehen ist, daB 


: al 2U—2) — ong 28M) 
(40) G(x, 936,431) =In-—_, —2F =n : : 
4 ny Pa) 4 [22a | 
t t 


fiir alle (x, y) aus & und By und fiir alle (x, y) ++ (&, 7) auf ©, analytisch ist, bleibt 
lediglich der Fall zu untersuchen, daB (x, y) in einer Umgebung von (€, 7) liegt. 
Zu diesem Zweck betrachten wir 


\ Vas th 2(Coj X—cos Y 
(41) G(X, Y) =In ee ) 


1 a@ und 7 bedeuten AuBen- und Innenseite von Gy. 
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in einer Umgebung von X = Y =0. Es gilt 


Co Xe $ a yeute 
1 L 
Rel aay a als Gl ae 2, | ) (2u+2)! 
X21 V2 
4 PG ede (ae 4) yrutr2 
42 gst 
( ) 2s (2u+2)! Xe y2 


oo u 
i ee ADs! Vx se Faly 


falls die letzte Reihe konvergiert. Dies ist in der Tat fiir alle komplexwertigen 
X und Y der Fall, da man bei Vorgabe einer reellen Zahl M>0 fiir alle |X| <M, 
|Y|sM 


co 
f=.2 an Sin M 


bekommt. Da die pent eda (42) mit 


pV) LLG ! a (X2 — V2) 4 565 OS ee Vest yay: 


beginnt, konvergiert beim Ubergang zum Logarithmus 


Cea 


1 2 2 1 4 22 eae 

(x? — v8) — 1 (x4 6x8 Y2 + Va) + 

in einer Umgebung von X = Y =0. Dann ist auch G(x, y; &, 7; ¢) in einer Um- 
gebung von x=&, y=y analytisch und bleibt damit U*— U einschlieBlich der 
Ableitungen beim Durchgang durch ©, stetig. 


§5. Aquivalenz des Problems mit einer Fredholmschen Integralgleichung 
zweiter Art 
Satz 6. Notwendig und hinreichend dafiir, daB es zu gegebenen Konstanten C,,, 
Doo und zu einer auf S, erkldrten Funktion v eine Stromfunktion Y mit den charak- 
teristischen Mittelwerten C.,, V9, und der Geschwindighkettsverteilung v auf Gy gibt, 
ist die Stetigkeit von v auf ©, und das Bestehen der Fredholmschen Integralgleichung 
zweiter Art 


(43) v2 oo [ Fods=2 (0,150 + 20,9). 
Cy 


Beweis. Wir schicken die Bemerkung voraus, da8B (37) und (38) ihre Giiltigkeit 
beibehalten, wenn sich U um einen in %, By, © regularen additiven Bestandteil 
von (36) unterscheidet. Gibt es nun zu gegebenen C5, 0.9, v auf Gy eine Strom- 
funktion Y mit den charakteristischen Mittelwerten c,,, ),, und der Geschwindig- 
keitsverteilung v auf ©, so ist Y durch diese Angaben vollstandig bestimmt und 
durch (20) in & gegeben. DefinitionsgemaB ist 0W/dn und damit v auf G stetig. 
Aus (9), (20), (36) und (37) folgt notwendig 


av 
on 


= Cy as o [ Fuds=— 
PT a LL ake 9 On 2n On 
Cy 


U 
2? 


on won er 
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und es entsteht (43) unter Benutzung der auf ©, geltenden Relationen 


Ox oy oy Ox 
(44) BAB GN 
on Os on és 


Seien umgekehrt Konstanten c,,, 0. und eine auf Gy stetige Funktion v gegeben 
und bestehe zwischen diesen GréBen die Integralgleichung (43). Dann wollen wir 
zeigen, daB die mit c,,, D.. und v auf ©, gebildete Funktion (20) eine Stromfunk- 
tion Y mit den behaupteten Eigenschaften ist. Zunachst ist Y in %& zweimal 
stetig differenzierbar und sind (2) und (3) in Y& erfiillt. Mit (20), (36), (37), (43) 
und (44) bekommen wir weiter 


(45) Og GS! UO oy Ox | o | Fuds—>=—9; 


én\|a ane, 2 = Yeo, 4 Bm — "99 an " Om On 

©, 
damit sind die an der AuBenseite von ©, genommenen Normalableitungen o'P/én 
auf ©, stetig. Wir betrachten nun Y durch (20) auch im Innengebiet By erklart; 
dann folgt in Verbindung mit (36), (37) und (45) fiir den Durchgang durch ©, 


é Uy 
on 


_ aw 
a ed(oi 


+ v= 0. 


a 


(46) 


Wir wenden nun den Greenschen Satz auf By an: 


Jas: 


7 


ae ee soll 


da ¥ in 8, ih oe und auf ©, stetig ist, igh wegen (46) 


lh (WH? + WY) dxdy=0 


und damit 


(47) Y= Constusdnidy: 
Weiter bekommen wir mit Hilfe von (20), (36), (38) und (47) 
(48) At ae Peer eoney 


und ist daher (4) auf Gy erfiillt. Unmittelbar zu iibersehen ist ferner, daB (20) 
die Bedingung (5) erfiillt. Damit hat (20) alle Eigenschaften einer der Betrachtung 
zugrunde liegenden Stromfunktion ¥. Dann besagt (45), daB Y auf ©, zu der 
ee Geschwindigkeitsverteilung v fithrt. Aus (20), (33), (34) und (35) folgt 
erner 


2 (D400 + Poo) = 2[, lim v(x, y) + lim »(x, y)] =v, 


Cully ¥ besitat gemaB (13) den gegebenen charakteristischen Mittelwert 
SchlieBlich bilden wir mit (20) ib 


bs Di sl pe aa! ee GEOG Se 
©, 
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Aus (21) und (35) erhalten wir fiir alle y die Grenzwerte 


lime x F) = ts 
x%—>-Loo i t 
und folglich auch 
li = = pe 6 ge 
(49) ee pe 0.) = fae a f Feds. 
G, 


Dann kénnen wir in (16) den Grenziibergang Z-> oo ausfiihren und bekommen, 
da der durch (49) gegebene Grenzwert des Integranden nicht mehr von 7 abhangt, 


1 
(50) C400 = Coo +4 févds. 
G 
Damit besitzt Y% gemaB (18) ebenfalls den charakteristischen Mittelwert c,,. 


Aus den Satzen 4 und 6 lesen wir unmittelbar die Eindeutigkeit ab, mit der 
Y durch c,,, v,, und v auf ©, bestimmt ist, wenn v stetig ist und der Integral- 
gleichung (43) geniigt; gabe es nadmlich zwei Funktionen Y mit den gleichen 
charakteristischen Mittelwerten c,, und v,, und den gleichen Geschwindigkeits- 
verteilungen v auf ©, so waren beide Funktionen notwendig durch (20) darstell- 
bar und somit identisch. Es werden schlieBlich durch alle Kombinationen c,,, 
v,., Vv, Wobei v auf W, stetig und (43) erfiillt ist, auch alle Stromfunktionen ¥Y mit 
den eingangs genannten Bedingungen erfaBt, da ja die Stetigkeit von v auf Cy 
und die Integralgleichung (43) notwendige Eigenschaften von YW sind. Wir 
brauchen uns im folgenden daher allein mit der Integralgleichung (43) zu be- 
schaftigen. Da iiberdies nur die Geschwindigkeit auf ©, und nicht die Strom- 
funktion selbst interessiert, konnen wir die Konstante c,, ganz auBer Betracht 
lassen und suchen daher bei konstantem, aber beliebigem »,, alle auf ©, stetigen 
Lésungen von (43). 

Eine Lésung von (43) ist jedoch notwendig stetig auf ©). Fiihren wir namlich 
die Parameterdarstellung x(q), y(qy) fiir ©y ein und setzen bei Beachtung von (1) 


4 a o(y) Ly 
a vp) (GME + FOP 


so geht (43) in 


2a 


652) oly) — 2 [ Kp, vi oly) dp =2[0.0,04(9) + %,99 (9) 


PRE A 
0 


iiber, wobei, wie sich anschlieBend herausstellt, 


Ver SC ie a ane aie ATOZ, — oo, 2H17) =F) 


1 Die Einfiihrung einer Funktion (gy) durch (51) wird (bis aufs Vorzeichen) von 
Isay [5] vorgeschlagen und tragt wesentlich zur Vereinfachung der folgenden Rech- 
nungen bei; (gy) kommt jedoch keinerlei physikalische Bedeutung zu, da keine 
Invarianz gegeniiber Parametertransformationen besteht. 
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ein fiir alle @, y stetiger und in 27 periodischer Kern ist; w(y) und damit v(g) 
sind infolgedessen im Falle ihrer Existenz ebenfalls stetig und in 2% periodisch. 
Die Stetigkeit von K(q, y; ¢) ist fir p+ y+ 2ka, k=0, +1, +2..., ohne weiteres 
evident, folgt jedoch auch fiir y= q+ 2ka aus der Existenz des (nicht mehr von 
t abhangenden) Grenzwertes 


ie tim Ke. 95) = ko) (ROP + 5)? 
mit 
(55) k(o) = &() ¥(~) — (9) ¥ (9) 1 


(Vi (oP + [v(p)]2)° 


als der Kriimmung von ©, an der Stelle g. Aus (10) und (51) erhalten wir noch 
fiir die Zirkulation 


(56) T= Jfa(g) dg. 


0 


§ 6. Allgemeine Lésung der homogenen Integralgleichung 
Um zu Aussagen iiber die Existenz und Mannigfaltigkeit der Lésungen der 
Integralgleichung (52) zu gelangen, hat man nach der Fredholmschen Theorie 
[6], [7], [8] die zu (52) transponierte homogene Integralgleichung 


27% 


(67) Mo) — se | Ko) uly) ay =0 


0 
auf die Existenz von nichttrivialen Lésungen «(g) hin zu untersuchen. 


Satz 7. Die Integralgleichung (57) besitet die nichttrivialen Lésungen 


(58) u(y) = const + 0. 
Beweis. Es geniigt offensichtlich, die Identitat 
2m 
(59) = [Kyitdp=1 


2 
0 


nachzuweisen. Wir untersuchen vorerst die Funktion 


an(y—n) 6, 2a(¥—8) _ 2n(¥—€) 22(y—n) 
(60) H(x,¥;€,n;) = asad 2 : 


co EB — gs 22=F) 


ae OO} 


fiir alle (x, y) aus YU, Bo, ©, wenn (€, 7) auf Cy liegt. Wahrend die Analytizitat 
von H fiir alle (x, y) + (&, 7) ohne weiteres einleuchtet, so gilt dies doch auch fiir 
(x, y) =(&, 7), wie anschlieBend gezeigt werden soll. Es gentigt hierzu, 


(61) H(X, Y) — ¥ SinX—X siny 
Gof X —cos Y 


* Auf Grund des vorausgesetzten Umlaufssinns von gy bedeutet positives bzw. 
negatives k(~) nach au8en hin konvexe bzw. konkave Kriimmung von @). 
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in einer Umgebung von X = Y =O auf Analytizitat hin zu untersuchen: 


lee) X2utT38 co y2urs 
Ont oe pig a aa 
=o (2u+3)! = (24+ 3)! 
u=0 u=0 
Jel, a ¥ 3 X22 on ny yeute 
1 
2 Gata 2,'— " eaea 
=p (2u+3)! X?+4 Y? 
u=0 
= X Y 5 1 X2u+2 +(= 1) yeure 
(62) Ao (2u+2)! x21 Y2 
~ 1 = ( = ee year tia) 
ay — (24+ 3)! = we 
Ss ¥ 1 by ( Neste ead 
nao (24+2)! \,&5 
= XY -S (PY xy) Bs), 


Diese Reihe konvergiert in einer Umgebung von X = Y =0, da alle in (62) vor- 
kommenden Reihenentwicklungen konvergieren bzw. durch die absolut konver- 
gente Reihe (42) majorisiert werden. Damit ist (60) fiir alle (¥, y) aus U, Bo, Gy 
und alle (,7) auf ©) analytisch. Dies gilt infolgedessen auch fiir die aus (60) 
folgenden Beziehungen 


10H a = sin 22) 

in Oh gy Bale—B) — og, 2A) 
(63) ee ar 

THA. wy 2 i i 

i OS) ot Gof eae S68 geet . 


Ersetzt man in (60) und (63) x, y, &, 4 bzw. durch x(q), y(@), x(w), y(), so sind 
(60) und (63) fiir alle wy, y und alle solche ¢>0 stetig, die nicht zur Bertthrung 
oder Uberschneidung der ©), €,,, €_,... fiihren; folglich ist mit (53) auch 


i aK __ (9) OH ky) OH 1 OH 
(Ab at toy tf 0% t ap 
fiir alle gy, y und alle ,,zugelassenen“ ¢ stetig. Dann ist es aber erlaubt, 
22 
1 
(65) J(yst} => | Keys tdy 


PHL, 
0 


unter dem Integralzeichen nach ¢ zu differenzieren, und wir bekommen unter 
Benutzung der Eindeutigkeit von H auf ©) aus (64) und (65) 
2% 2% 
Ey en RRO Cy rs ea | cao 
(66) pe Dae :" ale 2nt J Op eae 
0 0 


1 Der Verfasser verdankt diese Identitat einer Anregung durch Herrn Dr. K. Jor- 
GENS. 
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fiir alle y und alle zugelassenen ¢. Wir fithren jetzt in (53) fur i= pa hon 
k=0, +1, +2..., den Grenziibergang t— oo! aus: 


(y) [x(9) —*(y)]—4 (e) vy (9) —y (y)] 


lim K(g, y;t) = K(g,y; 0) =22 


67) EES [¥(~) —#* (y)]?+ Ly (e)—y (y) ]? 
d y(~)—y(y) 2 
39 are te ee 2 
ap 8 x (g)—x(y) ? 


der Grenzwert lim K(y, p-+2ka; t) unterscheidet sich nicht von (54), da ¢ nicht 
t+ co 

explizit in (54) vorkommt. Da schlieBlich (67) fiir p>qg-+2kza gleichfalls den 

Grenzwert (54) besitzt, ist Jim K(, yp; t) fiir alle q, py stetig. Fiir alle y folgt nun 

aus (65) und (67) ae 

(68) Lr. f(g; Ct 


t— oo 


Die Gleichungen (66) und (68) gestatten noch nicht, allgemein den Schlu8 zu 
ziehen, daB J(w; ¢) identisch den Wert 1 besitzt, da es z.B. im Falle ,,verzahnter“ 
Gitterprofile (Fig. 1) solche ¢ gibt, die gréBer als der gegebene Profilabstand sind 
und dennoch zu Berithrungen oder Uberschneidungen der €,,€.,,€_,... fithren. 
La8t man jedoch im gegebenen Gitterprofil €,,, €_,, €.;, ©_,... aus, betrachtet 
also an seiner Stelle cin aus @,€..,€_,,0,,,€_,... bestehendes Profil mit 
doppeltem Abstand 27, so sind jetzt fiir alle groBeren Abstande keine Beriithrungen 
oder Uberschneidungen mehr méglich. Daher ist mit (66) und (68) allein der 
SchluB 


(69) Jw; 2t) =1 


erlaubt. Durch Integration der aus (53) folgenden Identitat 


Sin y(v)1 


a [¥(~)—*(y)] 
t 


(70) K(g, p; t) =2K(,y; 2) —2“ arctg 
ae Sin 


erhalten wir dann aber unter Benutzung von (05) und (69) die Behauptung 
1 = 
2a | K(y, yp; 1) dp = Jp; t) = 2] (py; 2t) —1=1. 


Satz 8. Die Integralgleichung (57) besitzt nur die nichttrivialen Lésungen (58). 


Beweis. Sei y(p) eine Lésung von (57). Dann ist (@) fiir alle gy stetig und 
in 2a periodisch und damit auch auf G stetig. Fiihren wir in (57) die Bogen- 
lange als Parameter ein, so erhalten wir unter Verwendung von (27) und (53) 


7 i eset] We Aid es 
(71) Sp iene (od tee 
Cy 
i Hierbei wird die Tatsache benutzt, daB sich die ©, €.,, €_,... fiir hinreichend 
groBe ¢ nicht mehr beriihren oder iiberschneiden kénnen. — 


” Ae wy; co) ist der Kern der bekannten Integralgleichung fiir ein einzelnes 
Tori | 2]. 
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wobei ¢/én auf die Integrationsvariabeln é, 7 bezogen ist. Wir definieren sodann 
mit Hilfe von mw eine in YX und %, eindeutige harmonische Funktion 


(72) Meh = ~ wds. 
P ae 
Beim Durchgang durch ©, erleidet V einen Sprung 


(73) Vie-Vig=—-#, Vie V=—, 
wahrend die Normalableitungen stetig bleiben: 

ov Val ee 
(74) On\a Oni; i 


Die Beziehungen (73) und (74) leuchten ein, wenn man das entsprechende Ver- 
halten des in Y& und YB, eindeutigen Potentials 


(75) Vales Mele a | ae B aay) ds 


© 


einer doppelten Belegung « auf ©) heranzieht [13]; dann ist namlich 


V*(x,9) — Vix) =a [| 


0 


eg (In C 2F)) ds 


on 4 772 y? 


eine in YU, By, Cy analytische Funktion, da mit (40) auch 


OG 0 t2 
pe ee ey | = IE 
on on (In 477 7 “ 


fiir alle (x, y) aus %, Bo, Cy analytisch ist, wenn (&, 7) auf ©, liegt. Nunmehr 
ergeben (71), (72) und (73) 


(76) Viz=Ve, 0 


2 


Wir wenden den Greenschen Satz auf das zweifach zusammenhangende Gebiet ¥¢ 
mit den Berandungen Gy und Gy (Fig. 1) an: 


e of 2 OV OV 
om ffwarauesunaves= fv Sar fren 


©, 


WS 


Da AV in SF verschwindet, es wegen (21) und (72) auf dem oberen und unteren 


Teil von GF entgegengesetzt Boke Werte besitzt und eV/dn auf ©, stetig ist, 
folgt aus (76) und (77) 
(78) ff (V2+V2)dxdy= f VV,ds— f VV,ds. 
C5, rechts G, links 
LaBt man die Integrationswege CF recnts UNA CF jinks IN Fig. 1 in je einem parallel 


zur x-Achse gelegenen Periodenstreifen bzw. gegen + oo gehen und beriicksichtigt 
dabei die aus (27) und (72) folgenden und fiir alle y geltenden Grenzwerte 


lim (V9) (291 =0, 
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so geht die rechte Seite von (78) gegen Null, und man bekommt V =const in BF 
und somit 


(79) on 


auf ©. Wir wenden schlieBlich den Greenschen Satz auf ¥y an: 


(80) fy (VAV+V2+V2)dudy= Lae Jas. 


Da AV in By und wv wegen (74) und (79) auf ©) verschwindet und V]|; auf 
n |i 
Grund der aus (73) und (76) folgenden Beziehung 


(31) Vie 


auf ©, stetig ist, liefert (80) 
[f+ V2) dxdy =o. 
Bo 


Also ist V in 8» und damit auch V |; auf ©, konstant, so daB (81) die Behauptung 
us =const auf &, ergibt. 
Auf Grund der Fredholmschen Satze besitzt die zu (52) gehérige homogene 
und zu (57) transponierte Integralgleichung 
2n 


ne K(y, yp; t) yy) dy =0 


0 


(82) ¥(y) — 


wie (57) eime nichttriviale Losung y(g) und geht die allgemeine Lésung aus 
dieser durch Multiplikation mit einem willkiirlichen konstanten Faktor hervor. 
Jede Lésung von (82) ist fiir alle m stetig und in 2a periodisch. 


Satz 9. Fir jede nichttriviale Lésung y(p) der homogenen Integralgleichung (82) 
galt 
27 
(83) Sy (p) dp +0. 


Beweis. Es sei y(¢) eine nichttriviale Lésung von (82). Wir nehmen an, y (¢) 
erfiille die Bedingung 


2 
(84) Jv(p) 4p =0. 
Dann ist 
(85) 0 v() 


Ve @ME+D(OE 


nach Einfithrung der Bogenlange als Parameter auf @, stetig und geniigt wegen 
(43), (52), (82) und (84) den Bedingungen 


86 hey Ds sis 

(86) etl | Fods=o, 
G 

(87) Oasis. 
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Wir erklaren nun in Y% und Wy eine eindeutige harmonische Funktion U(x, y) 
durch die Vorschrift (36). Dann folgt aus (36), (37) und (86) 


ou 
(88) |= 


Da ferner U auf ©, stetig ist, liefert der Greensche Satz fiir das Innere 
ff (UAU + U2 + U2) dady= Ane Se) as, 
SMU SU Nandy =0, 
d.h. U ist im ganzen Inneren konstant, und (38) ergibt 
(89) Ul, ==const 


auf ©). Indem wir die aus (37), (87) und (88) folgende Beziehung 


(90) {> as= [ (52 


On 
C, C, 


g)ds =0 


voranschicken, wenden wir den Greenschen Satz auf Sf (Fig. 1) an 


Lae an | (on)] 


G 


ds; 


hier verschwindet insbesondere das Integral tiber ©) wegen (89) und (90), wahrend 
sich die Integration iiber ©j, wie man leicht sieht, auf den rechten und linken 
Teil von Cf reduziert: 


(01) [f (U24+U2)dxdy= f UU,ds— f UU, ds. 
By is rechts GG, links 
Wir erhalten aus (21) und (36) fiir alle y die Grenzwerte 
%—> +00 x 


wegen (35), (36) und (87) gilt auBerdem fiir beliebige y 


lim [x2 U,(x,y)]=_ lim aq | M(B 4 Z)eds=o, 
=x0° %—> +00 ze 


folglich auch 
lim | (U(x, y) U,(x, ¥)] = 


x—>+ 


LaBt man die Berandungsstiicke C% yecnts und Cotinks bzw. gegen + oo gehen, 
so liefert (91) 
U =const in 8 
und somit 
aU 
(92) San rane 


On \q 
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Dann ergeben (37), (88) und (92) 
0U 
on 


aU 


= —o9=0 
a S 
a on 


O 


auf ©). Wegen (85) ist schlieBlich fiir alle p 


d.h. »(q) ist im Widerspruch zur Voraussetzung die triviale Losung der Integral- 
gleichung (82). 
Eine unmittelbare Folge von (83) besteht darin, daB man ohne Einschraénkung 


(93) Srl) ag=' 


setzen darf; die Gleichungen (82) und (93) haben dann genau eine Lésung. Kénnte 
man ndmlich zwei voneinander verschiedene Funktionen y,(g) und y.(g) an- 
geben, die beide (82) und (93) erfiillten, so ware y,(~) — y2(g) eine nichttriviale 
Lésung von (82) mit 


I fate) — Yo(p) 4p =0; 


dies stande im Widerspruch zu Satz 9. Im folgenden verstehen wir unter y(q) 
stets diejenige Lésung der homogenen Integralgleichung (82), die der Neben- 
bedingung (93) geniigt. Die allgemeine Lésung von (82) ist dann durch const y (¢) 
gegeben. 


§ 7. Allgemeine Lésung der inhomogenen Integralgleichung 


Die inhomogene Integralgleichung 
27 
1 : 
(94) a(—)— 52 | Klp.vit) aly) dy =2%(9) 
0 


besitzt nach den Fredholmschen Satzen sicher eine Lésung, da die inhomogene 
Seite von (94) zur allgemeinen Lésung mw(y) der transponierten homogenen 
Gleichung (57) orthogonal ist: 


22 


a x (p) w(p) dg =const fi (y) dy =0. 


Die allgemeine Lésung von (94) setzt sich dann aus einer speziellen und der 
allgemeinen Lésung der zugehérigen homogenen Gleichung const y() additiv 
zusammen. Wegen (93) laBt sich der konstante Faktor vor y() durch die Forde- 


Tung 
20 


(95) J a(p)dp=0, 

0 
der eine spezielle Lésung «(p) der Gleichung (94) zusitzlich gentigen soll, stets 
cindeutig bestimmen. Andererseits gibt es genau eine Losung « (q), die (94) und (95) 
gehorcht. Gabe es zwei voneinander verschiedene Lésungen a (gy) und «(p~) von 
(94) und (95), so ware a (~) — a (p) im Widerspruch zu Satz 9 eine nichttriviale 
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Lésung der homogenen Gleichung (82) mit 


27 


4 [a (@) — a(y)|dy=0. 


Im folgenden verstehen wir unter «(p) diejenige spezielle Lésung der inhomogenen 
Integralgleichung (94), die die Bedingung (95) erfiillt. In gleicher Weise bekommen 
wir genau eine Funktion 6(g) durch 


(96) Bi) — 55 | Kw. wi Bly) dy=29(9), 
(97) [ B(p)dp =0. 


0 
Dann stellt aber 
Voo, x % (7) on U00,y B(¢) 
eine spezielle Lésung der inhomogenen Integralgleichung (52) dar, deren allge- 
meine Losung nunmehr durch 


(98) W(P) = Veo, x %(P) + Yoo, yB(Y) + Ty (9) 


gegeben ist; der willkiirliche konstante Faktor I’ bekommt darin wegen (56), (93) 
(95) und (97) die Bedeutung der Zirkulation. Auf Grund von (51) lautet die 
gesuchte allgemeine Lésung fiir die Geschwindigkeitsverteilung auf ©) schlieBlich 


(99) U(~) =2w, AY) + 2x, B(v) +I C(¢@) 
mit ee (9) a 

A (9) Te@F+D OF’ 
(100) Pia eee 


Vert ool 
v(@) 
TE@ PD (y)]? 


Da die Funktionen «(¢), B(yv), y(v) nur von der Geometrie des gegebenen Gitter- 
profils abhangen, gilt dasselbe auch fiir die ,,Grundlésungen™ des Problems 
A (gy), B(y), C(y), mit denen man die allgemeine Lésung (99) beherrscht. Die 
GroBen y(y), A (vy), B(y) sind dimensionslos, wahrend «(¢), B(p) eine Lange und 
C (y) eine reziproke Lange darstellt. 

Von Seiten der Praxis interessiert die Kenntnis der Geschwindigkeitsverteilung 
auf ©, bei gegebener ,,Anstrémung™ v_.. und Zirkulation J’. Es empfiehlt sich 
dann, (99) mit Hilfe von (14) in 


(101) v(9) =0-co, x A (p) + 0c,» B(p) +5, BY) + C(@) 


umzuformen. Die ,,Abstrémung“’ ».,, berechnet man aus (12). Haufig sind 
jedoch anstelle von J’ und v_,, An- und Abstrémung v_,, und v,,, vorgeschrieben, 
allerdings mit der erforderlichen Einschréankung v_.~/y = Vt00,1 =Vco,x- Fiir 
diesen Fall bekommt man die gesuchte Geschwindigkeitsverteilung aus (12), (43) 
und (99) in der Form 


(102) 0(p) = V0, (p) + ¥-o,y [2 B(Y) — tC (P)] + Y+00,y Le Bp) FEC ()I. 
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Mit v(y) bekommt man die an Gy herrschende Druckverteilung aus der 
Bernoullischen Gleichung 


PAY) =P aoe ay v(@) A 
uo) ws d—oco =<? ae 
wenn v_, =|v_.| den Betrag der Anstromung, p+. den Druck fiir x -+-F © 


and fi. +¢0 Poet = v,, den zur Anstrémung v_,, und zur Dichte ot 


gehorigen ,,Staudruck“ angibt. 


Von praktischem Interesse ist ferner die auf ©, wirkende Gesamtkraft &. 
Man erhdlt sie aus dem Kutta-Joukowskischen Satz 


(104) RK = oT" [v.05 f], 


fiir dessen Komponenten man aus (12), (43), (14) und der Bernoullischen Gleichung 
die Beziehungen 


IE 
K,= 9! = 0 |V70,9 2 55) = t% co Uh 0o) = = #400 Pc) 


ES — eo Pigs ey 02 Vo0, « (Vaoaly = % oh ee) 


(105) 


erhalt. 

In Hinblick auf die Anwendung geniigt somit die Kenntnis der Grundlésungen 
A(g), B(), C(g), mit deren numerischer Bestimmung wir uns jetzt beschaftigen 
wollen. 


§ 8. Interpolation und numerische Differentiation einer an Aaquidistanten 
Stellen gegebenen periodischen Funktion 


Fiir die numerische Rechnung seien die ©, beschreibenden Funktionen x (g) 
und y(gy) an den aquidistanten Stellen 


(106) %y = FE, oH 0, £1, £2... 


als x,—x(p,) und y,=y(,) gegeben, wobei N22 eine ganze Zahl bedeutet. 
Wegen 9,42n=P,+ 2% und damit 


(107) Xu +2N aa x (Py 22) =X ’ Vut2N == V (Pu - 22) rr. Vu 


gibt 2N genau die Anzahl der Konturpunkte auf ©, an, deren Koordinaten 
gegeben sind. Fiir die numerische Berechnung der Druckverteilung kommt es dann 
darauf an, allein mit den vorhandenen Informationen 1%»... %)y—1, Vo-»-Von—p t 
die Grundlésungen A (gy), B(), C(m) an den Stellen ,, also die Zahlenwerte 
Ay... Agy-1, By... Ban_1, Cy... Cgy_y zu ermitteln?. Bei der Berechnung des 
Kernes K (9, y; t) der gegebenen Integralgleichung taucht nun das Problem auf, 
die Ableitungen x(q) und y(g~) numerisch zu bestimmen. Wir befassen uns daher 


zunachst mit der Differentiation einer an aquidistanten Stellen bekannten periodi- 
schen Funktion. 


* Nicht zu verwechseln mit der durch (36) und (85) voriibergehend eingefiihrten 
Belegungsfunktion @, 

* Im folgenden ist mit /,, stets der Wert einer periodischen Funktion f(y) an der 
Stelle y, gemeint. Ist eine Funktion f(y, y) in beiden Verdnderlichen periodisch, so 
wird entsprechend T(Pu» Pr) =fuy Besetzt. 
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Eine in 22 periodische Funktion /(g), deren 2N Funktionswerte /, bekannt 
seien, werde durch eine Fouriersche Summe 


: N-1 
(108) = a, COS Mp + 2 b, sin up 
derart interpoliert, daB 
(109) ha = tn 


wird. Dann berechnen sich die Koeffizienten a,,, b,, eindeutig aus den Besselschen 
Formeln [14] zu 


, Nat 
age ie 


(110) a, +15, -+> fc rr pee Pe (NGS 


Wir interessieren uns fiir die Ableitung 1 (9), die in der numerischen Rechnung 


an die Stelle der nicht bekannten Ableitungen i (9) treten soll. Es folgt aus (108) 
und (110) 


i(y) = He {a0 + + (a, +ib,) eH e4 ay eine 


N— 
se) By. el! (pv—9¢) =f: | 


1 
u=l1 


u=0 


4 + ef (Pr—)) 


4—e'N (%—9) 
ees war} 


1 
A i! 
; N-1 — 
/, Re \( ae yet ae 
1 
7M {( 
1 


j, ctg —* sin N(p, — 9), 


2 | 


A 2N—1 : 

¥ me 41 SI 1 sin N(g,—¢) N cto Po” Seine a i 

/(@) are » if F ta Py—P g§ 2 (P ) ’ 
D 


strebt p>,“ =0...(2N—1), so erhalten wir unter Beachtung von (106) 


zunachst 
; 2N-1 
pe es ee 
v=0 
yeu 


denken wir uns nun die /,, durch /,,,2y=/, fiir beliebige w erklart, so kommen 


wir zu der bequemeren Formel 
2N-1 


(141) i x iO ee 


Arch. Rational Mech. Anal., Vol. 3 18 
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mit 
(112) p, =~ t=! cig BE SP ay 


fiir die wir nachtraglich die Beschrankung 4 =0...(2N—41) wieder aufheben 
kénnen. Wir vermerken noch 


2N-1 2 
(113) x, 1y=0 
Ke 
als Folge von (411) und (112). 


§ 9. Numerisches Berechnungsverfahren 


Will man die Integralgleichung (52) mit. Hilfe irgendeiner Quadraturformel 
in ein lineares Gleichungssystem iiberfiihren, so macht sich die Spur K(q, @; ?¢) 
des Kernes K(g, y; ¢) besonders stérend bemerkbar, weil in ihr nicht nur die 
ersten, sondern auch die zweiten Ableitungen von x(y) und y(gm) vorkommen. 
Abgesehen vom praktischen Mehraufwand, den die Bestimmung der %(y) und 
y() erfordert, stellt doch schon die numerische Berechnung von *(g) und y(q) 
aus (114) und (112) eine Fehlerquelle dar und gilt dies in erhdhtem MaBe fiir die 
zweiten Ableitungen. Aus diesem Grunde wollen wir die Integralgleichung (52) 
mit Hilfe der Identitat (59) so umformen, daB bei anschlieBender Anwendung 
der Rechteckformel alle Kernelemente, die die zweiten Ableitungen enthalten, 
herausfallen : 


(114) [ (KW. 9: 9 o() — K(p, v2) o(p)] dp = 2[P.., 4 (9) + P20, 5 (9): 


Denken wir uns die Integralgleichungen (82), (94) und (96) entsprechend umge- 
formt, so bekommen wir im Zusammenhang mit den Nebenbedingungen (93), 


~ 


(95), (97) fiir die Naherungen «(¢), B(~), y(p) das lineare Gleichungssystem 


21 ah ay 
Dy Kyo — Koi eC — Ko en-i 
u=0 
u+0 ~ a ~ 
ee raha R [ Bo Yo 
; 10 >; le ae eI oNenT ~ 
u=0 x )) 
Beas Xy By 1 
Dia ee ea : 
ke a 2N-1 oo B. 5 
QE _pP 2N-1 P2N-1 Y2N-1 
Ken-1,0 Ken-11 » Ea 2N-1 
u=0 
wM+e2N-1 
22 27 rst 27 
(115) is = 
2% 2Vo 0) 


2%en-1 2Ven-1 0 
O 0) 4\ 


Dieses System ist jedoch nur scheinbar tberbestimmt ; da namlich die Summe der 
ersten 2N Gleichungen des Systems wegen (113) identisch verschwindet, kann 
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man ohne Einschrankung der Allgemeinheit irgendeine der ersten 2N Gleichungen 
fortlassen. In (115) bedeutet K(9, y;t) diejenige Kernfunktion, die aus (53) 
hervorgeht, wenn man %(y) und j() bzw. durch ¥(p) und ¥(q) ersetzt. 

Die Berechnung der Kernwerte K uy tur w=» stellt einen wesentlichen Anteil 
des Rechenverfahrens dar. Den entstehenden Aufwand kann man aber noch 
weitgehend einschranken, indem man die numerische Bestimmung der K uy durch 
Anwendung der Additionstheoreme auf algebraische Operationen zuriickfiihrt. 
Zu diesem Zweck stellt man zu Beginn auBer den Ableitungen Sn Vu» =0... 
(2N —1), die Zahlenwerte 


20% 


Ehe=e aalts 0 (QIN), 
2a xy 
a ad tea p=Orw 2N —1), 
(116) cate 


25 
R,=cos="2#, w=0...(2N—1), 


S,=sin =, w=0...(2N—1), 


bereit. Mit Hilfe von (53) wird dann fiir uw,» =0...(2N —1), w=—», 


(117) K _ 22 Vn (ih, E,)"|—2%, F, E,(S, Ry—R, S,) 
wa of 1—2F, E, (R, Ry+ S, Sy) + (Fy Ey)? 

oder auch 

118 K Tess Vi Eto) = tt — 2%, By tela Ry S) . 

( ) ake t (2, F)* 24) F, (ty iy S55) 4 3 


wegen E,f,=1 sind (117) und (118) miteinander identisch. 

Bei Benutzung elektronischer Rechenmaschinen mit Festkomma-Arbeitsweise 
machen die Bereiche der E eC i Schwierigkeiten. In diesem Fall sollte man 
wenigstens die E,,, /, und die anfallenden Produkte £,,F, mit Hilfe eines Gleit- 
komma-Unterprogramms bestimmen; je nachdem, ob dann FE, bzw. E,, F, nicht 
gréBer als 1 ausfallt, kehre man mit FE, bzw. E,,f, zur Festkomma-Darstellung 


zuriick und benutze zweckmaBig (117) bzw. (118). 
Nach Auflésung des Gleichungssystems (115) findet man schlieBlich die 
Ndherungen der Grundlésungen A(p), B(), C() entsprechend (100) aus 


A, a0 (ZN 1), 
© Vat Th 
(119) By =e (je ON A); 
a SA 
ai ~ 
C,=— =0...(2N —1 
Gi commen ay ( ) 
HM be ue 


§ 10. Auswahl einer geeigneten Parameterdarstellung fiir die Kontur 


Wir haben bisher vorausgesetzt, daB fiir ©) eine Parameterdarstellung ~ (9), 
y(p) gegeben ist. In praktisch interessanten Fallen wird man jedoch eine solche 
Darstellung im allgemeinen nicht vorfinden; ©, wird vielmehr graphisch oder an 


18% 
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diskreten Punkten numerisch gegeben sein, die wiederum eine graphische Auf- 
tragung von G, gestatten. Wir beschreiben im folgenden eine Methode, nach der 
man Funktionen x (gy) und y(g) mit den gewiinschten Eigenschaften fiir bestimmte 
praktisch vorkommende Gitterprofile konstruieren kann?. 

Auf ©, gebe es zwei Punkte %9=(%o, Yo) und ty=(*y, Yn) mit den (in Gegen- 
uhrzeigersinn gerichteten) Tangentialeinheitsvektoren f)= (%o.1Vo).2 wae y= 
(xy, Yy), den (nach auBen gerichteten) Normaleneinheitsvektoren 1y= (Yo, — Xp) 
bzw. tyw=(yy, —*y) und den (zweckmaBig aber nicht notwendig maximalen) 
Kriimmungen k,>0 und ky >0, die sich unter den Voraussetzungen 


(120) by = (tp, Yo — tw) > 0, On = (tw, tn — Lp) > 0 


als ,,Hinterkante‘ bzw. ,,Vorderkante“ auszeichnen lassen. Nun stellt die Glei- 
chung 


(121) g(%,5 p) = (Mo, t — zo) cos* S — (ty, — Ey) sin? F =0 


sicher fiir jedes OX ma eine Gerade dar; es gilt namlich 


Cg \2 0g 2 4 1 
(FE) + (s8) => [1 — (mo, tty)] +> [1 + (Mo, ny)] 205? g, 
wobei die rechte Seite sicher nirgends verschwinden kann, da sonst tg= My ware 
und aus (120) dann b)+ by =(ty— ty, Yy—Ln) =O folgen wiirde. Speziell stellen 


(122) 2(4, 95.0) = (ng, FS) 0; ey) = "ek Se) = 0 


die Tangenten an die Hinter- bzw. Vorderkante dar, wahrend die Geradenschar 
(121) fiir 0<@m<a das von (122) eingeschlossene Gebiet schlicht tiberdeckt. Von 
@, wird nun die weitere Eigenschaft vorausgesetzt, daB die Tangenten (122) auBer 
den Beriihrungspunkten yy bzw. ry und alle dazwischen hegenden Geraden 
auBer je einem Schnittpunkt (Schnittwinkel = 0) mit der ,,Ober-“ bzw. ,,Unter- 
seite’’ von © keine gemeinsamen Punkte mit ©, besitzen diirfen. Ersetzt man 
dann in jedem Schnittpunkt mit der Unterseite y durch 2” —g, so erhalt man 
auf ©, eine eindeutige Zuordnung zwischen der Bogenlange und einem die Kontur 
im Gegenuhrzeigersinn durchlaufenden Parameter 0S gS 2a. Praktisch bedeutet 
dies, daB man N +1 Geraden g,= g(x, y; p,) =0, w=0...N, in die gegebene 
Zeichnung eintragt und die Koordinaten der Schnitt- bzw. Beriithrungspunkte 
Xu» Yur #=0...(2N —41), abliest (Fig. 3). 

Da ©, stetig gekriimmt ist, folgt die zweimalige stetige Differenzierbarkeit 
von «(y) und y() ohne weiteres fiir 0<m<a und a<g<2za, wahrend die 
Stellen »=0 und wy=2x bzw. xy und ry besonders untersucht werden miissen. 
Wir nehmen dazu einfachheitshalber an, daB die Konturkriimmung k(s) in einer 


* Grundsatzlich kénnte man natiirlich die Bogenlange als Parameter wahlen; dies 
diirfte jedoch deshalb unzweckmafig sein, weil man dann fiir aquidistante Parameter- 
werte an stark gekriimmten, wie an weniger stark gekriimmten Stellen der Berandung 
die gleiche ,,Punktdichte“ bekame, wabrend eine dichtere Besetzung mit Kontur- 
punkten in der Umgebung starkerer Kriimmungen wiinschenswert erscheint. 


* Akzente bedeuten Differentiation nach der im Gegenuhrzeigersinn wachsenden 
Bogenlange s. 
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Umgebung von yy bzw. s =0 analytisch ist. Dann gilt 


k= ky + ks ++ 


und folgt hieraus der Ortsvektor in kanonischer Darstellung [15] 


(123) tp plas tog ltgs? ae eet hs tgs? toes 


Durch Einsetzen in die Geradengleichung (121) erhdlt man unter Benutzung 
von (120) fiir die Abhangigkeit zwischen s und in einer Umgebung von s =0 
bzw. py =0 


1 jG ¢ 
(5 fos? + & kos? + 2) cos* 2 = [by — (fp, ty) s +---] Sue 


TY ge0\ge0 \gu0 


Fig. 3. Konstruktion einer Parameterdarstellung «(¢), y(q) flir eine gegebene Kontur ©) im Falle 2N = 12 


Beiderseitiges Radizieren liefert, wenn man beachtet, daB s und @ gleiche Vor- 
zeichen besitzen, 


Yias(t Ro § f+-}cos © = YOy[1— Su tad s +... sin’. 
1 2D 2 a8 


| 6 Ro 2by 2 
dann ist aber 
ki 
aga Heit aa Gass 
@ = 2arctg os 
2bn ques (to, ny) s ate =e 
(124) is 2by 


ee ae | i eokelevene) 
ap be "16% |  2on 


analytisch in einer Umgebung von s =0, und ebenso ist s eine analytische Funk- 
tion von @ in einer Umgebung von y=0. Eliminiert man s in (123) mit Hilfe 
von (124), so hangt auch x in einer Umgebung von » =0 analytisch von @ ab. 
In gleicher Weise findet man regulires Verhalten von y an der Stelle p=a. 


a4.) 
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§ 11. Beispiele 
Vorgegeben seien reelle Zahlen a>0, >0, Uo0,41 Voo, y» I’. Dann definiert 


(125) vf — 10% = 0_0, 03 (0) + Yo, 0h) + Ped) 
ee t 1 1 1 
* 
wt(¢) = =a ) 
ie t a 8 PE 
t 1 1 1 
26 wht) = : ) 
ae ee eee ee 


* 1 1 : 1 1 
wt =—1-(—- + PA 
smi( C+ altg as pays Soe = 


in der ¢-Ebene auf dem Kreis 
(127) t=ae'%, OY S28, 
und iiberall auBerhalb desselben mit Ausnahme der Stelle € = —a Ctg one eine 


Potentialstrémung v*=(v¥,v}). Der Kreis (127) ist auf Grund des bekannten 
Spiegelungsprinzipes Stromlinie bzw. Teil einer solchen, denn es gilt 


(128) we (ad) = 16" WW, (yr ee, 3, 
wobei die 
; Sin — sin ¥ 
Mh) 2a 4ma 4na y 
ee + Gin 7 C08 x 
t 1 
@ = i ee. 
Cee Wal) 2na [ Ana = 4nG ) 
Cp} ; + Gin 7 08 7 


ie! 1 F 
3 (Z) 42a (14 4za ares AiI0. 
oj rs + Gin COS 
reelle GréBen sind. 
Denkt man sich nun die ¢-Ebene entlang der reellen Achse von — co bis 


—a tg aa aufgeschnitten, so tiberfiihrt die konforme Abbildung 


. t 22a © ve CRUE: 
4 (6) = —— AY Sa Sy : = woeA 
(130) zg=x+tty 3 in (Goi : + — Sin 
den Kreis (127) der ¢-Ebene in ein Gitterprofil der z-Ebene mit der Teilung ¢, 
den Begrenzungen a..=--a und der Parameterdarstellung fiir C, 


«= shin |/Coj 424 4 Gin 44 cos y, 
(131) 


‘ Fiir den arc tg sind die Hauptwerte zu nehmen. 
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wobei © mit wachsendem 0< ¥< 2 genau einmal im Gegenuhrzeigersinn durch- 
laufen wird. Die Geschwindigkeitsverteilung (125) transformiert sich in die 
Potentialstrémung » =(v,, v,) der z-Ebene gemaB 

(132) Vz — 1 Vy = Ugo, 7 Wy (2) Tt Yo0, y We (2) + Iw, (z) 

mit 


(133) w, (2) = SE wh(e) = 22 (¢ + aGtg 274) w¥ (0), v=1, 2,3. 


24 
t 
ubergangen x—> + oo in der z-Ebene Aquivalent sind, folgt aus (126) und (133) 


Da die Grenziiberginge €-—> co und (> — a Ctg in der ¢-Ebene den Grenz- 


(134) lim W,(z) = —1, 


hi a Was 
pete (2) 2s 


Damit enthalt (132) einerseits die Berandungen ©), ©,,,€_,... als Stromlinien 
bzw. als Bestandteile von Stromlinien und erfiillt zum anderen fiir «> ++ o die 
Bedingungen 


. a Va 
(135) lim Uy = Yoo, x lim By = Voo,y S75 


e+ x—> +00 


folglich ist (132) die allgemeine Losung des Problems und ergeben sich die Grund- 
l6sungen aus der Definition (99) als die Werte der w,(z) auf ©). Beachtet man, 
daB jede Geschwindigkeitsangabe auf ©, bei Umstr6mung im Gegenuhrzeigersinn 
einen positiven Wert haben soll, so folgt fiir die Grundlésungen aus (128) und (130) 


(136) es W,(x) = ane a |/soi a8 | Sin aes cosy, v=1,2,3. 
Sin 4 


t 


Die Parameterdarstellung (131) ist nun insofern fiir die numerische Rechnung 
ungeeignet, als sie die Symmetrie des Gitterprofils beztiglich der y-Achse nicht 
ausnutzt. Wir wenden daher das in § 10 beschriebene Verfahren zur Ermittlung 
eines Kurvenparameters 0< gS 2a an und erhalten mit ry=at, ty=1, Iy= — At, 
Ny = —t aus (121) 


(137) x(p) =acos @. 


Eliminiert man x aus (131), so folgt 


(138) sini hae 3s arc COS 


1 Fiir den arc cos sind die Hauptwerte zu nehmen. 
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wobei das positive Vorzeichen fiir 0S psa und das negative fiir ti pS2z gilt. 
Entsprechend ergeben sich aus (129), (131) und (136) die Grundlésungen 


DI oe 2IO (P) 
A(y) = — 2 Cig a: a recsaee 
ae 
(139) B= 2g 
Sin 
t 
, Si a 
OG) ee 
Sin —— 


t 
Insbesondere ist (138) in einer Umgebung von g =0 analytisch und besitzt dort 
die Potenzreihenentwicklung 


1 62a 2 8 SATE 2 EE iy 224, 
(140) y(9) =a} a(t ’ Ctg i a t Jett] Bee ree 


im Zusammenhang mit (55) und (137) folgt weiter fiir die Kriimmung in einer 
Umgebung von gy =0 


2 
(141) k(p) = 77 Cty =4% — ome ght, 


d.h. an der Stelle gy =0 und damit aus Symmetriegriinden auch an der Stelle 
2% 27a 
; Stg ; 


t 
, si =—(0,4 und —— =0,2 
2a 2a 


ist das Gitterprofil (137) und (138) in Fig. 4 bis 6 dargestellt. Das numerische 
Verfahren wurde im Falle 


g =a nimmt die Kriimmung von ©, die maximalen Werte 


an. 


os : : . Fe é Gone 
Fiir drei verschiedene ,, Teilungsverhaltnisse me =| 


—=1 mit 2N=12 und 2N = 24 Konturpunkten, 


2a 
— =0,4 mit 2N=20 und 2N = 40 Konturpunkten, 
= =0,2 mit 2N=28 und 2N =56 Konturpunkten 


durchgefiihrt. Die Koordinaten x,, y,, die exakten Werte (139) fiir die Grund- 


losungen A,, B,,C, und deren Naherungswerte 4,, By C, als Resultat der 
numerischen Rechnung sind in den Tabellen 1 bis 3 aufgefiihrt ; dabei geniigt es aus 


Symmetriegriinden, diese GréBen fiir y =O... S anzugeben. Es zeigt sich, daB 


man fiir kleiner werdende Teilungsverhiltnisse ¢/2a und damit fiir wachsende, 
mit a dimensionslos gemachte Maximalkriimmungen 


ak(0) = ae ci 


eine gréBere Anzahl von Konturpunkten benétigt, um eine gleichbleibende 
Genauigkeit zu erhalten; zum anderen ergibt sich bei Halbierung der Schritt- 


weite bzw. Verdoppelung der Konturpunkte eine sehr rasche Konvergenz gegen 
die strenge Lésung. 
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Die numerischen Rechnungen wurden mit der mittelschnellen elektronischen 
Rechenmaschine G2 (ca. 20 Operationen pro Sekunde) des Max-Planck-Instituts 
yA 

Cy, 


SI 


Fig. 6. Gitterprofil mit einem Teilungsverhdltnis ¢/2 a=0,2 


fiir Physik und Astrophysik in Miinchen durchgefiihrt. Ohne Ausnutzung von 
speziell gegebenen Symmetrieeigenschaften dauert die vollstandige Bestimmung 
der Grundlésungen mit der G2 z.B. fiir 2N =24 Konturpunkte ca. 80 Minuten. 
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On the Uniqueness of Compressible Fluid Motions 


JAMES SERRIN 


One of the fundamental questions which should be answered concerning any 
problem of applied mathematics is whether it is well set, that is, whether solutions 
actually exist and whether they are unique. We shall be concerned here with 
the initial value problem for compressible fluid flow, and we shall study in 
particular the uniqueness of its solutions. An exact statement of the problem 
will be found in §1, while the main theorems are given in § 3 (viscous fluids) 
and § 4 (non-viscous fluids). Stated roughly, these theorems assert that fluid 
motion in a bounded region “ = ¥Y (t) is uniquely determined by its initial velo- 
city, temperature, and density distribution together with certain boundary 
conditions. The form of the boundary conditions is especially interesting for 
non-viscous, non-heat conducting fluids, it being found, in fact, that such con- 
ditions are superfluous at any point where the fluid is leaving Y at a relatively 
supersonic speed. This result is intimately associated with the existence of 
characteristic manifolds in the flow region, and allows us to prove for spatial 
flows a theorem of the type given by CourRANT & FRIEDRICHS for one-dimensional 
motions. 

The uniqueness of the initial value problem for incompressible viscous fluids 
was proved as long ago as 1929 by E. FoA. His proof, which applies equally 
well to the non-viscous case, was rediscovered and presented more elegantly by 
Do.ivzeE (see § 72 of reference [15] for a resumé of this work). More recently, 
GrRaAFFI showed that the initial value problem for compressible fluids is unique, 
provided that the pressure and density satisfy a piezotropic relation p =/(@) 
with /’(@) >0. Our work is an extension of GRAFFI’S, in that we make no assump- 
tion concerning the equation of state of the fluid, beyond requiring it to satisfy 
certain universal thermodynamic conditions. This generalization is non-trivial 
from both the physical and mathematical point of view, for it is well known that 
a relation between pressure and density is usually incompatible with the motion 
of a tri-variate fluid, whether or not it is viscous. Moreover, treatment of the 
general case requires that we use the full set of equations governing the fluid 
motion, including the complete energy equation. 

The proof makes use of the classical energy method, developed originally 
by HapAMARD, ZAREMBA, and FRiEDRICHS & LEwy in connection with purely 
mathematical problems, and by REYNOLDS, ORR and FoA for the study of 
incompressible fluid motions. Perhaps the most novel feature of our treatment, 
in comparison with these, is the systematic utilization of the transport equation 
for the derivation of energy integral identities, thus taking into account in a 
natural way the special form of the convection terms in the flow equations. 
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Though a more abstract treatment does seem possible, this has been avoided on 
the grounds that it would cloud the fundamental issues and make the whole 
paper more difficult. 

In view of the rather strong dependence of the viscosity and heat conduction 
coefficients on temperature in actual liquids and gases, we have included in the 
final section of the paper a brief discussion of this case. Also, although we do 
not treat heterogeneous fluids here, our results remain valid even in this case. 

The reader should observe that the uniqueness theorem for non-viscous fluids 
specifically assumes the flows in question to be continuously differentiable. The 
proof can be extended without difficulty to cover flows with first order singular 
surfaces (Mach waves), but in its present form it does not apply to motions 
admitting shock waves. Such an extension would be very important, since there 
are numerous cases where the only possible solution involves a shock wave, or 
in which shock waves occur in the initial data, and in none of these is there 
reason to expect anomalous behavior so far as uniqueness is concerned. Further 
discussion of this point will be found in § 4. 


1. Statement of the problem 


We consider a fluid flow occupying a finite region “~ =¥Y (¢), and subject to 
certain prescribed conditions on the boundary Y of ¥. Our purpose is to answer 
the following question: given the velocity, temperature and density distribution 
at some initial instant, does this uniquely determine the subsequent motion, or 
might it happen that two different motions could be compatible with the same 
initial conditions? —To make the problem precise, we shall now set down the dif- 
ferential equations governing the fluid motion, and the exact boundary conditions 
which we shall consider. These take a somewhat different form depending on 
whether the fluid is viscous or not. 


A. Viscous compressible fluid 


The motion of the fluid is governed by the following standard equations of 
fluid mechanics, 


do : 
a7 tt edive= 0, (1) 
d ‘ : 
0 =ef+divT, Ti (pA die) Bead), (2) 
E ; 
Die q=—xegradT. (3) 


In these equations the symbols have their usual meanings, except perhaps for 
the stress tensor J’, the deformation tensor D, the specific internal energy £, 
and the heat flux vector gq. The notation is that of ordinary Cartesian or Gibbsian 
vector analysis, so that, for example, 


(div T);=7,,;, I =identity matrix, 7T:D= T;;,D;;- 
The material derivative d/dt is of course a differential operator and is defined by 


d 0 
ap st ppae’ Veted, 
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The reader unfamiliar with equations (1)—(3) will find a complete derivation 
and explanation in the memoir of C. TRUESDELL, The mechanical foundations of 
elasticity and fluid dynamics, J. Rational Mech. & Analysis, vol. 1, 1952. 

To equations (1)—(3) must be added certain equations of state relating the 
thermodynamic variables o,f, E, and T. For this purpose we shall take 


pb=h(e,T) and E=E(g, T) (4) 


as the defining relations; here #(o, T) and E(o, T) are given twice differentiable 
functions of the variables g and T, whose explicit form depends of course on the 
particular fluid in question. In order to simplify the proof somewhat, we shall 
actually treat in detail only the case where the second equation has the special 
form 

Ser nec Consts= 0: (5) 


the proof can be carried out equally well in the general case provided that the 
specific heat c,=(0E/0T), is positivel. It may be added that (5) is quite well 
satisfied by most real gases and liquids at ordinary conditions. 

On the boundary of W we suppose the velocity distribution to be given 
(this will be the case in particular if YW is bounded by material walls to which 
the fluid must adhere). In addition to prescribing the velocity on -Y, it is 
necessary to assign certain further conditions in order to determine the fluid 
motion. Various possibilities present themselves, depending on whether at a 
particular boundary point the fluid is entering ¥, or leaving ¥, or is constrained 
by a material wall. For convenience in tabulating these cases, let m denote the 
outer unit normal to ¥, and G the outward normal velocity of . Also, let 


U=v-n—G (6) 


be the relative normal speed of particles at the boundary. Then we have the 
following table of 
Boundary Conditions. At all points of the boundary Y, and at all instants 
of time, the velocity vector 1s given. Furthermore, 
a) At all points where U<O the absolute temperature and the density are 
prescribed. 
b) At all points where U>O the absolute temperature 1s given. 
c) At all points where U=0 one of the following conditions 1s assumed : 
cl) The temperature is prescribed. 
c2) The heat flux q-n ts prescribed. 


1 The condition (@E/éT),>0 is a general attribute of all real materials, as shown 
in textbooks on thermodynamics. 

It is also possible to treat heterogeneous fluids which behave as inert mixtures; 
in this case the defining equations (4) become 
dn 
=ai? 
where ” is a parameter specifying the constitution of a given particle. The additional 
details necessary to handle cases (4) and (4’) can readily be supplied by the reader. 


p=p(e,T,n), E=E(e,T,n), = 05 (4’) 


274 JAMES SERRIN: 


c3) The heat flux is assumed to be proportional to the difference between the 
fluid temperature T and a given wall temperature Ih, that 1s, 


q:n=k(T— 1), RO. 


Cases a), 6), c) above correspond respectively to fluid entering ic : fluid 
leaving Y, and fluid constrained to move tangentially to 7. Condition cd) is 
the familiar Newtonian law of cooling. We remark that any one of conditions 
c1)—c8) could be used at a boundary point where U>0, though only the first 
one seems to be of any physical interest. 

This completes the specification of the differential equations and the boundary 
conditions for a viscous compressible fluid. 


B. Non-viscous compressible fluid 


Assuming that the fluid is devoid of heat conduction as well as viscosity, 
the differential equations (1)—(3) reduce to the familiar system 


22 + edivv=0, (7) 
o = of —gradp, (8) 
oo 

a (9) 


where S denotes the specific entropy. The thermodynamic equation of state is 
conveniently assumed given in the form =#(o, S), where the function f(g, S) 
is continuously differentiable and 
(42) ty (10) 
doa/s 


Ss 
corresponding to a veal speed of sound in the fluid. 


At each point of the boundary we suppose the normal velocity component 
v-n to be assigned (this is a weaker condition than was necessary for viscous 
fluids). The case of most interest is naturally that where the fluid is bounded 
by material walls and the normal velocity component is determined by the motion 
of these walls. In this case, in fact, no further boundary conditions should be 
imposed. The general situation is given in the following table, where U=v -n—G. 


Boundary Conditions. At all points of the boundary Y, and at all instants 
of time, the normal velocity component v-n is given. Furthermore at all points 
where U<O the complete fluid state—velocity, density, and entropy?—is assigned. 

It is understood in this statement that no further conditions are to be imposed 
at points where U=O. On the other hand, the fact that the entire fluid state 
must be prescribed at points where fluid is entering Y will appear reasonable 
upon a moment’s reflection. 

If heat conduction is present, the differential equations and boundary condi- 


tions take a somewhat different form. It is convenient to defer this case until 
a later section. 


* In place of the entropy one could of course assign the temperature or the pressure 
and thus indirectly determine the entropy. 


Uniqueness Theorems for Compressible Fluids DYS 


C. The imtial value problem 


For both problems A and B above, a solution of the initial value problem is 
defined to be a set of continuously differentiable functions satisfying the given dif- 
ferential equations and the boundary conditions, and reducing when t =0 to assigned 
initial distributions*. These solutions need not be continuable for all times f>0: 
it is not our intention, nor is it necessary in the proofs, to restrict consideration 
to solutions which exist for all time. Concerning the boundary of “, we assume 
it smooth enough so that the divergence theorem is applicable at each instant, 
and also so that the developments of the following section are valid. It would 
be pedantic to state these agreements with any greater precision. 


The solutions which we consider will be understood without further comment 
to have positive density at every point of the closure of YW. 


2. The transport equation 


Let F(x, y,z,t) be a continuously differentiable function defined over a 
region ¥ = ¥Y (¢) occupied by a moving fluid. In the sequel we shall have frequent 
occasion to use the transport formula 


4; | eFav= es dv heUFaa, (14) 
t Sf 


where @ and w are the density and velocity of the fluid motion, and U denotes 
as usual the quantity v-n—G. To prove (11), we begin with the rather obvious 


identity * 
d é 
aif eF =/3 (ef) + Herc, (12) 


3 We do not consider “weak solutions” or solutions involving shock waves. On 
the other hand, with slight modifications it would be possible to deal with flows 
having first order singular surfaces. 


4 Formula (12) is well known and physically obvious: it says that the rate of 
change of the integral foF arises as the sum of internal changes of oF together with 
changes due to the fact that the boundary of Y is moving with normal speed G. 
A proof may be given as follows: 


Suppose the moving region Y consists of those points (7, v,z) such that 
aT (AR OMGME) ie V9 SA (HOG; 2) 2 WGN) 5 


where a, b,c range over some closed region # in R%, and f/f, g,h are continuously 
differentiable functions of all their variables; suppose also that the mapping from ¥ 
to Y is one-to-one with non-vanishing Jacobian (Y is a differentiable homotope 
of #), and that the divergence theorem is applicable to Y at each instant. If the 
“velocity field’ (%,7,z) is denoted by w, then it can be shown that 


si fe-f 2 Down, p= 9%. y,2,%), 


(cf. [15], § 4, or [16], p. 53). Since w-n=G is the outward normal speed of the 
boundary of ¥, this completes the proof of (12). 

It is evident that the conditions laid down on Y in the preceding proof embrace 
all reasonable configurations. In any case, the regions which we shall consider in 
the rest of the paper are supposed to belong to the class for which (12) holds. 
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in which (as in the following work) the conventional volume and surface infini- 
tesimals are omitted, and the integrals are understood to be extended over either 
W or F as the context indicates. With the help of the divergence theorem, one 
sees that the right hand side of (12) is equal to 


fle (oF) + div(eF »)| +P eF(G—v-n). 


On the other hand, 


) : OF 00 : - dF 

© (ef) + divieF v) =¢ \ar +0. grad F} ips Fie + div (o v)| es 
using the equation of continuity (1) or (6),and the definition of the operator 
d/dt. Combining the results of the last three displayed lines immediately establishes 


the validity of (11). 


3. The main uniqueness theorem for viscous fluids 


With the conditions and agreements of §1 understood, we can now state 
the fundamental uniqueness theorem for compressible viscous fluids. 


Theorem 1. Let the coefficients of viscosity and heat conduction be constants 
satisfying the inequalities 
20, 23) = Oe ae 


Then there can be at most one solution of the initial value problem described in 
§ 1A. Moreover, if x =0 no boundary conditions on the temperature (or heat flux) 
need be wmposed at points where U=O. 

The important case where 34+ 2u=0 is not included in this theorem, but 
seems to require a more delicate treatment, see § 7. It may be remarked that 
Theorem 1 applies also to flows occupying unbounded regions, provided that 
certain boundary integrals at co can be assumed to vanish. 


~ 


Proof. Let (v,o, T) and (@, 6, T) be two solutions of the problem under 
discussion. We shall show that they are in fact identical. For convenience in 
notation, a tilde over a flow quantity will mean that the quantity is to be evaluated 


~ 


for the flow (v,@, T). In addition, for an arbitrary flow quantity F we define 
f= PF; 


J2 


thus in particular, v’ = — v, o'=6 —@, etc. The following identities will be used 
frequently throughout the paper: 
AB—AB=AB'+A'B, (13) 


dF dF dF’ ; S 
a ae ae -grad F. (14) 


The first of these is obvious, while the second follows from an easy calculation 
using the definition of the operator d/dt. Note that dF’/dt means 


- +o grad?” 
and not (dF/dt)’. 
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Since both flows satisfy the equation of continuity, we obtain by subtraction 
and use of the identities (13) and (14) 


£0 +! gradé + g dive’ + o' div® =0. (15) 
Multiplication by 0’ yields 
a (S 0’) = — (o’v'- grado + oo’ divv’+ o'2dive) = — TI, 


where J is defined by the last equality. Now using the transport equation (11) 


there results 
d 1 oath, i 12 
ai} ee fel heVe (16) 


This formula will be important in the sequel. 


Again, since both flows satisfy the equations of motion, we have by sub- 
traction 


dv 1 ~ Caen” , . / 
(+e - grad @) + 0 =a 04 1 diy T’. 


Forming the scalar product of this equation with the vector v’, and rearranging 
terms, yields 


EG Peer [o’v’- (a — f) +ov'-D-v'+T':D')]+div(T’-v’), 


where for brevity we have put v’?=v’-v’ anda=dv/dt. Finally, with the help 
of the transport equation and the fact that v’=0 on the boundary of W, one 
obtains 


ats ov t= — [| [o'e". f)+evr': D-v'+T': D’). (17) 


Turning next to the energy equation, we observe that in virtue of (5) it can be 


written in the form 
aT 


Oa = W—divg (W=T:D), 


where without loss of generality we have set c,=1. Applying the method already 
described in connection with the previous equations, the reader will have no 
trouble in proving the identity 


af eT t= f(T’ +a' -grad T’ — oie iy Te -grad T] — 


(18) 
—p[la nts eur. 

The next step in the proof is to estimate the various terms which appear 
on the right hand sides of equations (16), (17) and (18). For this purpose it is 
convenient to restrict attention to some frxed time interval OS¢S<t, where t 
is an arbitrary positive number (if the solutions are known to exist only for a 
finite time 0<¢X7), then t is naturally chosen less than t)). We shall let N 
be the generic notation for an upper bound; the number N will then not be the 
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same in each estimate, though it will of course be possible to determine its size 
at any stage. The expression oJ appearing in (16) has the following estimate, 
valid for OSS: 


|oI| =|ee'v’- gradé + oo’ dive’ + go'*divd| 
<N(o?+0'2) +20" (0 = dive). 


Here the basic tool is CAUCHY’s inequality 2ab<a?+b?, and « is an arbitrary 
positive number to be fixed later; N depends on ¢ and on bounds for the magnitude 


~ 


of (v, 0, T), (®,6, T) and their derivatives during the time interval 0OX¢St. 
In addition, because of the boundary conditions the term @ Ug’? in (16) is 20. 
Inserting these estimates into (16) yields the inequality 


af _ 12 2 2 Spee 

<< <[Nle a ea er O<t<t, 
where it is convenient to use the notation 

S=Sso0", H=Jfzov%, L=JzoT”. 


In considering equation (17) it will be advantageous to write T in the form 
T=— p1+YV, whereV is the viscosity tensor, V= AOI+ 2uD. Thensince I:D = 
div v =O, we have 

TDS OLY Ds 


It is now easily verified that 


£* S[N(@2+02+5") + [(cO%—V':D), 0StSr. 


Turning finally to an estimate for the right hand side of (18), we see that 
w= T:D—T:D=T:D'4+T :D 
=T:D'+(—?'+10) 0+ 2uD’:D, 
whence for 0X<#<r1, 
ie a al =|): D'—pT OLA TOL 2nT D:D 
<N(T'?+ p’2) + 602+ D':D’). 


(19) 


Furthermore q’- grad T’= —x grad T’?<0, and on the boundary of Y, in all 
cases, 9QU T’*=0 and T’q’-n=0. When the above inequalities are inserted into 
(18) there results 


Fes {Net + +9247) + fe(O2+D':D’). 
Adding the various estimates for the derivatives of f%, % and Y, and using the 


obvious inequality p’?< N(0’?+ T’2), gives easily 


£(g+ H+ LYSN $+ H+ L) + | [e(30%+ D': D') —V': D'). 


This completes the main estimate. It is necessary now to consider the term in 
braces in the last integral. 
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The following identity becomes obvious if the matrix D is written in diagonal 
form, 


3D: D=GO67+ A? (A =/(4, dy)® + (dy — dz)? + (dg d,)’) ; 


here 0 =div v =Trace D, and A is the octahedral invariant defined in terms of 
the characteristic values d,, d,,d, of D. From this identity follows 


VD =107%+ 2u D:D’ =4{(34 + 2u) O24 2A}, 
and also 
6302+ D’:D') = < (100/24 A’2), 
Thus if ¢ is chosen small enough, we have 


£(30'2+ D’: D’)<V':D’, 
and hence 
F(PLHA+ AEN P+ LLL), 05t&r. 


This differential inequality has the integrating factor e %’, and the solution 


eee a Pee ee ee StS r. 


Since the integrals J, % and ¥ are zero when ¢t =0, it follows that they remain 
zero throughout the whole time interval in question. But this forces v’= 9’ = 
T’=0, and so the two flows are identical at least until =r. Finally, t being 
arbitrary, the flows are identical so long as they exist, and the theorem is com- 
pletely proved. 

Another consequence of the final inequality is that solutions of the initial 
value problem, so long as they remain uniformly bounded, depend continuously 
on the initial data. 


4. The main uniqueness theorem for inviscid fluids 


We consider here the most important case A = =x =O, and refer the reader 
to §1 for a description of the problem. The principal result is the following 


Theorem 2. Let the equation of the state p=p(o, S) be of class C* for all values 
of the variables, and suppose that 


op 
| 00 \, 4 (20) 
Then there can be at most one solution of the initial value problem described in 
§ 1B. Moreover, if for some solution of this problem we have 


U=w-n—GEZc, = (Ze) 


on a subset S* of the boundary, then the boundary conditions on S* are redundant 
with respect to the given data; that is, the conditions assigned on the remainder of 
the boundary, together with the given initial values, uniquely deternune the solution. 


(To the author’s knowledge, the only previous uniqueness theorem for com- 
pressible non-viscous fluids with equation of state =#(o, S) is that given by 
Courant & FriEpDRICHS ([Z], pp. 82—84 and pp. 48—53, 70—75). There is a 
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large body of literature concerning hyberbolic systems of equations in more 
than two variables [5], but this work does not seem to have immediate application 
to our problem.) 

Before commencing the proof of Theorem 2, there are several points of related 
interest which we wish to discuss. 

To begin with, the second part of the theorem includes as a special case the 
result that initial values over a region # uniquely determine the flow (assumed 
continuously differentiable) throughout a certain “‘hypercone”’ in (x, y, Z, t) space 
with its base on &. If for geometric clarity we consider the motion to take place 
in two spatial dimensions, then this conical domain of determination will appear 
as in Fig. 1. In order to determine the exact shape of the cone it is necessary 
to know one solution of the initial value problem; the lateral boundary, or 


ey 


x 


Fig. 1. Domain of determinacy for the region & (two spatial dimensions) 


mantle, of the cone is then characterized by the condition that U=c at each 
of its regular points®. In view of the independent importance of this result, we 
shall restate it in the following precise form. 


Corollary. Consider a solution (with continuous first derivatives) of equations 
(6) —(8) in a region B* of (x, y, 2, t) space bounded by the plane t=0 and a mantle 
S* on which Uc. Suppose another solution with continuous first derivatives is 
given in R* which assumes for t =0 the same values as the first one. Then the two 
solutions are identical in B*. 


Another point which should be made in connection with Theorem 2 is that 
the given initial value problem may easily have no solutions at all, or only 
solutions which break down after a finite length of time. This occurs most 
obviously if we simply allow the boundary to penetrate the initial value cone 
A*, but it may also occur even when the boundary conditions appear quite 
innocuous, as has been shown by Luprorp in a very interesting paper. 


* The equation U =c specifies a Monge cone at each point of (x, y, z, t) space where 
the solution is known. The determination of the lateral boundary thus involves 
solving a first order partial differential equation; the mantle is then the inner envelope 
of the family of Monge conoids whose vertices lie on the boundary of R. Since the 
relation U =c can be shown to be a characteristic condition of equations (6)—(8), 
it follows that the lateral boundary is a characteristic manifold of the solution (ef. (21, 
pp. 374—376 (plane flow), [12], pp. 112—114, or [15], § 51; in the last named reference 
the characteristic condition is derived in exactly the form used here). 

The Monge cone noted above has obvious affinities with the Mach cone in (C53) 
space. Indeed the latter is just the projection of the Monge cone; this projection 
is real if and only if the flow is supersonic at the point in question. 
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The fact that initial values uniquely determine the subsequent flow throughout 
a portion of the flow field, and the fact that other solutions break down after 
a finite time, arises, of course, because we are considering continuous solutions 
of a hyperbolic system. If we allow solutions with shock waves, then the initial 
values no longer dictate the subsequent motion, and the breakdown of a continuous 
solution no longer disqualifies the future motion from consideration. For this 
reason, it may be conjectured that the first statement of Theorem 2 remains 
true for flows with shock waves. The work of LAx, OLEJNIK, and LADYSHENSKAYA 
(see reference [/7] and the accompanying bibliography) may provide the clue 
to a deeper understanding of the situation. 


5. Proof of Theorem 2 


~ 


Let (v, @, S) and (, 6, S) be two solutions of the initial value problem under 
discussion. Also, let Y* denote the (possibly empty) set of boundary points 
where the basic solution (v, 9, S) satisfies Uc. In order to prove Theorem 2, 
it must be shown that the two solutions are identical, and this must be done 
without recourse to the assigned boundary conditions on the set 4*. 

Retaining the notation introduced in the proof of Theorem 1, and applying 
the methods developed there, it is clear that from equation (7) follows 


af sour =— flee @—f+er'- D-v'+v'-gradp']—Gt evo’? (21) 


(this should be compared with the analogous equation (17) ; note that the boundary 
term in (21) cannot be assumed to vanish). The energy equation dS/dt=0 
likewise yields 
eS" +v'- grad S=0 (22) 
and 
d 1 ¥, 2 ete / , E CS £e ze y 4 2 
aifses 2 fese grad S h5 evs (23) 


Now let 0StS1 be an arbitrary time interval throughout which both solutions 


(wv, 0, S) and (®, @, S) exist, and let N be the generic notation for an upper bound. 
Then we have the following estimates valid for OS¢#StT: 


leo’ @—P|SNio'? +0), [ev D-v' 
|e S’v’- grad S| < N(v’? + S"); 


= Nv, 


and (24) 


~ 


here N depends on bounds for the magnitude of (v, e, S), (6, @, S) and their 
derivatives during the time interval OS/St. 

The remaining term in (21), namely wv’ - grad f’, requires a separate treatment. 
The algebraic effort in this estimate can be greatly reduced by agreeing that an 
equality sign with a dot placed over it means equality modulo terms which contain 
products and squares of the quantities v', e', and S’. Since such terms can be 
majorized by inequalities of the form (24), it is clear that this procedure retains 
at each stage the essential quantities. In any case, setting 


4-8) 
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we have v’- grad p' + v'- (c? grad g’ + A grad S’) 
— (c29' + AS’) SHE BAR PSD (25) 
she aera gifs ee o + div(ctg’ + AS) e" 
rH 


where the final step is a consequence of equation (15). Now 


1 , ao Be 
Ga hare nie *); 20) 


and using equation (22), 
4 POO ae Oe ny 1 > 


When (26) and (27) are inserted into (25) there results 


c2 


/ Jie d 
v'- grad p ~enlog 


g?+ 40'S) +div(o’ + 4S')e 
Q 


Substituting this formula into (21) and then applying the divergence theorem, 
the transport equation and inequality (24) yields the inequality 


d 1 oven Oe ppg ‘\s [Net +0 
af sel ea 2 Se ase (or 2 


~ Feu (v L _ 0’? 4 s 0’ S') 2(c2g' +48’) vn], 


S’2) 
(28) 


this being valid for O<¢SXt. In addition, from (23) and (24), 
d fi 2 12 Lo\ ga b ail 12 
af pests [Nwt+s')— hl ous’ (29) 


Now let # be a large positive number, whose value will be fixed in a moment, 
and define quantities J and K by 


12, OP. 4 ror / 
jav4 el ra igo as ghee bie 
Ke oU J -2(co -+-AS’) om, Of 24. 
We observe that # can be chosen so Nici that 
1 OS = (t+ SEND (30) 
for OStS<T; indeed, recalling that c>0, we may take 
= Max (225 +5). 
ese OG |) Dy 


With this choice of #, it follows from the boundary conditions (p. 274) that 
K20 at all points of Y. Even more, however, K=O on the set S* irrespective 
of the boundary conditions. For on L* we have U>c, so that either 


KZ gc] +20'(2e +AS)=Q 


or 


K20¢] —20'(%o’ -AS)\—0,, 
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Since both the quadratic forms Q are easily shown to be non-negative definite, 
the required inequality follows at once. 


With the results of the preceding paragraph in mind, let us multiply inequality 
(29) by # and add it to (28). Since K=O, the resulting boundary integral is non- 
negative, and we obtain 


d / / , 
Fil OFS} Met o'*+ $4); 
therefore, with the help of (30) and the fact that 9 and c are positive, 


afersNfes. OStsSt. 


It follows by integration and use of the initial conditions that { @ J=0, and then, 
using (30) again, that v’=o0’=S’=0. This completes the proof of Theorem 2. 

In the following section we shall deal with the previously excluded case of 
an inviscid fluid with non-zero heat conduction. 


6. The excepted case A=yu= 0, x>0 
The type of fluid considered in this section, having zero viscosity but positive 
heat conduction, is perhaps rather unlikely to occur in a physical problem. 
Nevertheless, the results of this section are quite interesting from a mathematical 
viewpoint, and will also be useful in a later part of the paper. 


The differential equations for the problem are obtained by specializing (1) —(3) 
to the forms 


“2 4 edivy=0, (31) 
oa” = ef —gradp, (32) 
oe =— pdive —divg. (33) 


As in the treatment of viscous fluids, we shall consider only the equations of 
state A 
p=p(9,%), E=«,T, (34) 


where c, is a positive constant. It will also be assumed that the pressure-density- 
temperature relation is twice continuously differentiable and satisfies the ine- 
quality 
op 

eee > 0, (35) 
corresponding to a veal isothermal speed of sound. Inequality (35), like the 
earlier conditions (GE/0T),>0 and (0f/0@);>0, 1s a general attribute of all 
real materials. Moreover, it is always the case that 


Gar > ae (36) 


\ 


see [15], §30. The boundary conditions appropriate to equations (31) —(33) are 
as follows: 
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At all points of the boundary S, and at all instants of time, the normal 
velocity component v -n is given. Furthermore 


a) At all points where U<O the complete fluid state, velocity, density, and 
temperature, 1s assigned. 


b) At all points where UZO the temperature 1s given. 


c) At all points where U=0, one of the conditions c1)—c3) listed on p. 273 1s 
assumed. 


These conditions being understood, we have the following uniqueness theorem 
for non-viscous, heat-conducting fluids. 


Theorem 3. There can be at most one solution of the initial value problem 
corresponding to the equations and boundary conditions described above. Moreover, 
if for some solution of this problem we have 


WHO) = 6 2G; a=(ZP) 


on a subset S* of the boundary, then the boundary conditions on S* are partly 
superfluous with respect to the given data: that is, the conditions assigned on the 
vemainder of the boundary, together with the temperature on S* and the given 
initial values, uniquely determine the solution. 


Exactly as in the previous section, this theorem implies the following corollary: 
Consider a solution (with continuous derivatives) of equations (31)—(33) in a 
region &* of (x, y, z, ¢) space bounded by the plane ¢=0 and a mantle /* on 
which Ua. Suppose another solution with continuous derivatives is given in 
&A*, which assumes for t=O the same values as the first one, and whose tem- 
perature agrees on Y* with the temperature of the first solution. Then the two 
solutions are identical. 


Since this situation is very much like that of the preceding section, we confine 
our discussion to the remark that the Monge cone for the present problem lies 
inside the Monge cone U=c, in virtue of (36). Thus the region &* here is “‘larger”’ 
than the region of determination for a non-viscous, non-heat-conducting fluid, 
assuming that the two can be compared at all. 


Proof of Theorem 3. The difference between the present problem and the 
problem 4 =A =x =0 lies entirely in the energy equation. Equation (21) then 
holds exactly as before, but in place of (23) we now must deal with 


a faer =f |a-grad T’— 0’ 7’  _ 9T'v'. grad F — r' "| — 
(37) 
—O|T'q'-n +i our), 


where ®= # div v (this identity is proved in the same way as (18)). A perusal 
of the various terms on the right hand side of equation (21) shows that, as before, 
the only difficult quantity is wv’. grad p’. Using the “dot” notation to denote 
equality modulo squares and products of v’, 9’, and 7’, and setting B =(0p/0T),, 
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we have (cf. formulas (25) and (26)) 


v’- grad p’ =v’. (a* grad 9’ + B grad T’) 
=—a@o'divv'+ Br’. grad T’ + div (a? 0’ v’) 


. d a , , / = A 2 
= Ale 0 ) + Bv’- grad T’ + div(a?0’v’). 


Substituting this formula into (21), and then applying the divergence theorem, 
the transport equation, and the usual kind of estimates, we obtain 


af fad (+4 2) [N( (24 y'24-772) + [el grad T’|[2— 
1 / 2 , Ces? 
= hg Alaa ask ea 


this holding as always for OStXt. 


(38) 


> 


Turning next to equation (37), one has the following identity for the term 
Oe: 
T'®'=T'(pdive)’ =T'(pdive’ + p’ div’) 
=div(pT'v’) —pv’-gradT’— T’v’- gradp + T’p' dive. 


If this formula is inserted into (37), then from the divergence theorem, the 
boundary conditions, and the usual kind of estimates, we get the inequality 


4 fhores[Net+e2+79+fe—*|gradTP (39) 


Now define 


K= U(o? 0"? + a? 0?) + 20a%0'v'-n, on S. 
Clearly K=0 in virtue of the given boundary conditions. Moreover 
|2eao’v'-n| S 2a] ag'|-|Qv"| 
Sa(g?v'? + a’ "), 


so that K =>0 on S* irrespective of the boundary conditions. Adding inequalities 
(38) and (39), and using the fact that K2O, gives without difficulty 


dM 
dt 


SNA+ [(2¢ = »)\erad I"), 021=r, 
where 


Aa fore Serer), 40 


Setting «4x and solving the resulting differential inequality leads easily to 
the identical vanishing of v’, 9’, and 7”. 


7. The excepted case 3A+2yn=0 
As we have already remarked, Theorem 1 does not apply in the important 
case when 3A+2y=0, that is, when the Stokes relation holds. However, this 
case is covered by the following 
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Theorem 4. Let the coefficients of viscosity and heat conduction be constants 
satisfying u>o, x>O, 3A+2u=0, 
and suppose also that (0p/8e)p>0. Then there can be at most one solution of the 
initial value problem described in § 1A. 

Proof. Because of the similarity between Theorem 4 and Theorems 1 and 3, 
we can make considerable use of formulas already developed in the preceding 
sections. In fact, since the only difference between equation (2) and equation (32) 
is that the former has an extra term div V on the right, we have, by applying 
to equation (2) the same development which led from (32) to (38), 

cs i 
a fselvt+ a o'*) < [Nl + 98+ Te) + felgradT |? 
in this inequality, the boundary conditions and the fact that V’:D’=20 have 
been employed. Besides this estimate, which takes the place of the inequality 
for d%#]dt in the proof of Theorem 1, it is necessary to find an estimate for 
f4oT'? corresponding to the earlier inequality for d #/dt. To this end, we observe 
that from (19), 


PLATT DP POLAG TOL201 DD 
=v'- T.grad T’+ 10 v'- grad T’+2uv’- D. grad JT” 
+ {squares and products of v’, p’, and T’} + {divergence terms}. 


Inserting this expression for 7’ ¥”’ into (18), and applying the divergence theorem, 
the boundary conditions, and the usual sort of estimates, yields 


a {ser ts [Not +v24p2+7%)+ [(e—x |gradT’ 


Now adding the two preceding inequalities, there results 


2 


44 NA+ [(2e—*|gradT’2, 0StSt, 


where -@ is defined by (40). The rest of the proof is standard. 


Remarks. In addition to the various combinations of inequalities restricting 
A, 4, and x which have already been treated, there are two others which deserve 


notice, namely EON On ee (41) 


and 
0; 324-20 =057) eo (42) 

Case (41) can be discussed exactly as in Theorem 1; when x =O the appropriate 
boundary conditions are those in § 1B, and when x>0 the appropriate boundary 
conditions are those of §6. The remaining case (42) does not seem to be acces- 
sible by our methods. 

Cases (41) and (42) exhaust all the situations for which the constitutive 
equations (2), and (3), satisfy the thermodynamic restrictions 


V:D20, q-gradTSo. (43) 


It would be interesting to consider the initial value problem for more general 
constitutive equations obeying (43) or some modification of (43) 


Uniqueness Theorems for Compressible Fluids 287 


8. Temperature dependent viscosity and heat conduction 


As mentioned in the introduction, the coefficients A, “, and x are rather 
strongly dependent on the temperature in actual liquids and gases. For example, 
in the case of air 


RT? ba 
w= 7 Conse: const = 200, 


while for water ~(100°)/u(0°) =6.314% The following theorem shows that no 
significant qualitative changes are brought about by treating non-constant 
viscosity and heat conduction. 


Theorem 5. If the coefficients of viscosity and heat conduction in Theorem 1 
ave assumed to be differentiable functions of 0 and T, satisfying 


BOF RRA SOP or S00 


then the conclusion of Theorem 1 remains valid. 


Proof. If one repeats the steps in the proof of Theorem 1, he is led with 
only a few modifications (specifically, in the estimation of T’ Y”’) to the inequality 


x (4+ H+L)SN P+ H4L) + [ [6302+ D': D')—V’': D'+q'-gradT’]. 


On the other hand, 


~ ~ 


VD =(AV'+H0)0' + (2u D’ + 2u’' D): D' 
=10'? + 2uD':D'+R,, 
where | R,|<N(0’?+ 7T"?)+ (024 D’:D’). Similarly, 
gq’: grad T’ = (— x grad T’ — x’ grad T) - grad T”’ 
= —x|grad 7’ |?-+ Rp, 
where R,< N(0’?+ T’?) + |grad7’|?. It is evident then (by the argument used 
near the end of the proof of Theorem 1) that if ¢ is chosen small enough, 
6302+ D’: DD’) —V':D'+q'- grad T’S N(o’? + T’?). 
Thus we have 
I (F+HLL)SNI+H+), 
and the theorem follows at once. 


Note: The above research was supported in part by the United States Air Force 
Office of Scientific Research under Contract No. AF 49 (638)-262. 

6 Water, being an incompressible fluid insofar as its usual hydrodynamic properties 
are concerned, does not properly come within the scope of this paper. Nevertheless, 
by slight modifications Theorem 5 can be shown to hold also for incompressible fluids. 
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